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Abstract 



In this article we review the present understanding of azimuthal and single spin asymmetries 
for inclusive and semi-inclusive particle production in unpolarized and polarized hadronic collisions 
at high energy and moderately large transverse momentum. After summarizing the experimental 
information available, we discuss and compare the main theoretical approaches formulated in the 
framework of perturbative QCD. We then present in some detail a generalization of the parton 
model with inclusion of spin and intrinsic transverse momentum effects. In this context, we exten- 
sively discuss the phenomenology of azimuthal and single spin asymmetries for several processes in 
different kinematical configurations. A comparison with the predictions of other approaches, when 
available, is also given. We finally emphasize some relevant open points and challenges for future 
theoretical and experimental investigation. 



1 



Contents 



1 Introduction 

2 Experimental results on azimuthal and transverse single spin asymmetries 

2.1 Transverse single spin asymmetries for p^p — > h + X processes 

2.1.1 E704-Fermilab results 

2.1.2 BNL-RHIC results 

2.2 Azimuthal asymmetries in polarized semi-inclusive deeply inelastic scattering . 

2.2.1 HERMES results 

2.2.2 COMPASS results 

2.2.3 JLAB-CLAS results 

2.3 Azimuthal asymmetries in unpolarized processes 

2.3.1 SIDIS processes 

2.3.2 Drell-Yan processes 

2.3.3 e + e" — > nn + X processes 

3 Theoretical approaches to single spin asymmetries 

3.1 Parton model approaches including intrinsic transverse motion 

3.1.1 QCD developments 

3.1.2 The impact parameter picture 

3.1.3 Models and constraints for TMD distributions 

3.2 Twist-three effects in collinear pQCD 

3.3 A semiclassical model 

4 Spin asymmetries in the generalized parton model 

4.1 Polarized cross sections and helicity formalism 

4.2 Spin and TMD parton distribution functions at leading twist 

4.3 Spin and TMD fragmentation functions at leading twist 

4.3.1 Quark and gluon fragmentation into an unpolarized hadron 

4.3.2 Quark fragmentation into a spin 1/2 hadron 

4.4 Helicity amplitudes for the elementary process ab —> cd 

4.5 Kernels for doubly polarized processes 

5 Phenomenology of single spin asymmetries 

5.1 Unpolarized cross sections and SSA's in pp — > C + X 

5.1.1 Pion, kaon and heavy meson production 

5.1.2 Prompt photon production 

5.1.3 Transverse A polarization in unpolarized pp collisions 

5.2 SSA's for double inclusive production in pp collisions 

5.3 Azimuthal asymmetries and SSA's in Drell-Yan processes 

5.4 e + e~ annihilation in two nearly back-to-back hadrons 

5.5 Semi-inclusive deeply inelastic scattering 

5.5.1 Azimuthal dependence in unpolarized SIDIS 

5.5.2 Transverse SSA's in SIDIS 

5.5.3 Phenomenological extractions of Sivers and Collins functions 

6 Conclusions and outlook 



1 Introduction 



Polarization phenomena in hadronic processes deeply challenge high-energy physics theories and per- 
turbative Quantum Chromodynamics (pQCD). It has been known since a long time that experimental 
results for several polarized processes can hardly be explained by means of available theoretical ap- 
proaches. Some relevant examples are: i) The large spin-spin correlation measured in pp elastic scatter- 
ing and its unusual behaviour as a function of the relevant kinematical variables; ii) The so-called "spin 
crisis" , and the fact that quarks only account for a small fraction of the total proton spin; Hi) The large 
transverse polarization of hyperons produced in unpolarized hadronic collisions; iv) The large transverse 
single spin asymmetries (SSA) measured in inclusive pion and kaon production in polarized hadronic 
collisions and in semi-inclusive deeply inelastic scattering (SIDIS). Due to space limitation, in this work 
we will not cover at all the topics i) and ii). The rich, low-energy aspects of topics Hi) and iv), which 
may be confronted only with nonperturbative approaches, will also remain uncovered. A comprehensive 
overview on these subjects and on related topics can be found in several review papers, addressing both 
their theoretical and experimental aspects [II [31 IS IS 13 El El El El IIOl IIIl [T^ 1131 1131 113 IISl IIZI - Here we 
will focus on the high-energy regime and the treatment in the context of pQCD approaches of transverse 
SSA's and, marginally, of the transverse polarization of spin 1/2 particles in unpolarized processes. 
The SSA for the inclusive process A^B — >• C + X is defined as 



where der^ stands for the invariant differential cross section, Ecd 3 a^^ /d 3 pc, for the production of 
hadron C, with 4-momentum p 1 ^ = (Ec,Pt,Pl), m the scattering of a hadron A off an unpolarized 
hadron B, with A upwards (downwards) transversely polarized w.r.t. the production plane. More 
precisely, in the AB center of mass frame, with A moving along the +Z direction and C produced in 
the XZ plane, the upwards (downwards) transverse spin directions are respectively along ±Y. The 
SSA, or analyzing power, A^, is also denoted as "left-right" asymmetry since, by rotational invariance, 
the number of left going hadrons C when A is downwards polarized, da^, is the same as the number 
of right going hadrons C with A upwards polarized. Notice that left (right) is defined w.r.t. the plane 
spanned by the momentum and spin directions of the polarized hadron A looking downstream its 
direction of motion. 

The transverse polarization measured in the AB — > + X process, where C is a spin 1/2 hadron, 
is defined, analogously to An in Eq. flT|), as Pt — (dcr^ — da^) / dcr unp , where now the arrows in 
stand for the transverse spin directions for the final hadron C. 

Most of these processes show peculiar and unexpected dependences on the relevant kinematical 
variables: the rapidity, y = (l/2)\n[(E c +PL)/(E c —pL)], or the pseudorapidity 77 = — (1/2) lntan(#/2), 
the Feynman variable, xp = Pl/Pl'^i an< ^ Pt, where 9 and pl {pt) are respectively the production angle 
and the longitudinal (transverse) momentum of the observed particle w.r.t. the direction of the colliding 
beams. The magnitude and behaviour of these polarized observables are very difficult to account for 
in pQCD approaches, as well as in nonperturbative models, and have motivated a huge amount of 
theoretical work. Although significant progress has been done, as it will be illustrated in the rest of 
this review, much remains to be clarified. Moreover, a new family of second-generation experiments 
are starting to produce new sets of data with larger statistics, at higher energies and with an enlarged 
coverage of the kinematical variables. For these reasons, we will try to collect the main ideas on 
SSA's and give an account of their present status, both from the experimental and, in more depth, the 
theoretical side. 

The study of SSA's is strongly related to our understanding of the structure of hadrons and their 
spin and orbital angular momentum content in terms of partons. This leads naturally to a special 
attention to the role played by parton distribution and fragmentation functions in pQCD approaches. 
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In this review we will then present, in some detail, a new class of parton distributions, beyond the 
usual collinear ones, resulting from the inclusion of spin and intrinsic transverse momentum effects. We 
will also show how they can help in describing the observed SSA results. In particular, two spin and 
transverse momentum dependent (TMD) mechanisms, introduced in the early 90s, will play a special 
role: the Sivers effect, describing the asymmetry in the azimuthal distribution of an unpolarized parton 
inside a transversely polarized hadron [HI [19]; the Collins effect, describing the azimuthal asymmetry 
in the fragmentation of a transversely polarized quark into an unpolarized hadron [20, [21] . 

Among the topics not covered in this review, that are of some relevance in the context of SSA's, we 
mention the generalized parton distributions (GPD) and the dihadron fragmentation functions (DiFF). 
GPD's represent another generalization of ordinary parton distributions and indeed could provide a 
unified parameterization of many different aspects of hadron physics. We refer the interested reader to 
the following papers [22j [231 12H [251 E6J [271 12E1 [29], which emphasize the different aspects of the physics 
encoded in these quantities. We also mention their possible connection to the TMD distributions, 
discussed in Refs. [501 EH E21 E21 EH ESj- Concerning DiFF's (describing the probability that a parton 
hadronizes into two hadrons plus anything else) we recall their role in spin studies, in particular for 
the extraction of the transversity distribution function in double inclusive processes in pp collisions 
and SIDIS [361 E7J, an d the azimuthal correlations in e + e~ annihilations. Some seminal papers on the 
DiFF's, with their definitions and properties, are [381 (211 EHl SOI SH S2l S3l S3] . 

The outline of this review is the following: in section [2] we present an overview of experimental 
results on SSA's, including also the azimuthal asymmetries observed in unpolarized collisions, that 
are intimately connected to our main subject. In section [3] we discuss QCD theoretical approaches to 
SSA's, which represent nowadays the most reliable attempts towards the understanding of the observed 
results. Their most recent developments are also presented. We will discuss in particular two pQCD 
approaches: one based on the extension of the collinear factorization theorems, at leading twist, with 
the inclusion of intrinsic transverse momentum effects; another one based on a twist-three collinear 
formalism. As a representative example, a short account of a semiclassical model will be also given. In 
section [H we will give a detailed description of the first pQCD approach within the helicity formalism. 
Its phenomenological implications and a comparison with the results of other theoretical approaches 
are collected in section [51 Finally, perspectives and conclusions are gathered in section [61 

This work is to a great extent based on our longstanding and fruitful collaboration with M. Anselmino, 
M. Boglione, and E. Leader; without their contribution it would have been impossible. We are also grate- 
ful to A. Kotzinian, S. Melis, and A. Prokudin for valuable collaboration in the last years. Many people 
have contributed with enlightening discussions and clarifying remarks to this work, before and during 
its completion; in particular it is a pleasure to thank A. Bacchetta, D. Boer, C. Bomhof, L. Gamberg, 
Y. Koike, A. Metz, P. Mulders, F. Pijlman, P. Schweitzer, W. Vogelsang, F. Yuan, for their theoreti- 
cal insights; H. Avakian, L. Bland, F. Bradamante, M. Grosse Perdekamp, D. Hasch, S. Heppelmann, 
J. Lee, A. Martin, L. Nogach, R. Seidl, F. Videbaek, for their kind support on the experimental side 
and for providing us with preliminary results of their experimental collaborations. 

2 Experimental results on azimuthal and transverse single 
spin asymmetries 

A comprehensive account of all experimental information available today on polarization phenomena 
in high-energy hadronic physics is well beyond the scope of this review. Forced by space limitation, 
in the following pages we will only summarize the most significant results, limiting ourselves to those 
concerning the main subjects of this paper: transverse single-spin and azimuthal asymmetries. The 
selection presented somehow reflects our personal perception. We apologize for those papers and topics 
not included and discussed here. The interested reader may fortunately find complete and comprehensive 
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accounts on this subject in Refs. [JJ, [21 [91 [T3] . 

Whenever possible we will refer to published work. For ongoing experiments we will also quote 
available preliminary results. We will not discuss details concerning experimental setups, detector 
performances, etc., referring to the quoted papers for all technical aspects. In order to avoid duplicates, 
we will not show plots of the results in this section, deferring instead to section [51 where a selection of 
relevant results, and their comparison with theoretical predictions and fits, will be presented. 

2.1 Transverse single spin asymmetries for p*p — > h + X processes 

Starting from the late 70s and during the 80s a number of experimental collaborations presented early 
results on SSA's for the pp — > h + X process. For its impact on the later developments we mention, for 
instance, the work on the transverse A polarization by Bunce et al. |45j . SSA's significantly different 
from zero were found for the p^p — > it + X process in the forward production region (large positive 
xf)-, at relatively small transverse momentum, pp. These results were in accord with previous extensive 
measurements performed for similar exclusive channels (e.g., irp — > irp) and were interpreted as purely 
nonperturbative effects. Additional measurements on the pp dependence of these results showed inter- 
esting behaviours, difficult to understand and explain. All these experiments were fixed-target ones, 
at relatively low cm. energies and low pp (below 1.0-1.5 GeV/c). Therefore, no serious attempt was 
made to understand these effects in terms of pQCD involving hard elementary scattering. Moreover, a 
well-known paper by Kane, Pumplin and Repko [IS], reinforced the common wisdom that spin effects 
should become negligible at high-energy scales. 

2.1.1 E704-Fermilab results 

During the 90s, the E581 /E704 Collaborations at Fermilab (henceforth indicated as the E704 Collabo- 
ration) set the stage at a different level, because of the larger cm. energy available (y/s ~ 20 GeV), the 
improved separation between the xp and pt kinematical dependence of the SSA's and the small but 
significant extension of the pp range covered. The E704 Collaboration collected SSA data for pp and 
pp collisions using (for the first time) polarized secondary p, p beams and for several produced particles 
(mainly pions, but also photons, A and rj particles). The experiment covered essentially two different 
kinematical regions: 0.2 < xp < 0.6 and pp in the range 0.2 2.0 GeV/c (the forward production, or 
beam fragmentation, region); \xp\ < 0.15 and pp values up to 4.0 GeV/c (the central rapidity region). 
The main results gathered by the E704 Collaboration can be summarized as follows: 
1) Large SSA's (up to 30-40% in magnitude) have been measured for inclusive pion production in the 
large positive Xp region, for polarized proton and antiproton beams colliding with an unpolarized hy- 
drogen target (33, HHl HH1 [50]. For both neutral and charged pions the SSA was found to be almost 
zero at low xf, up to xp ~ 0.3, where it starts raising with increasing xp, reaching values of ~ 15% 
and 30 — 40% in size for 7r°'s and 7r s respectively [see Fig. [3j p. ETJ. The measured asymmetries 
are roughly similar in magnitude but opposite in sign for positive and negative pions (and positive in 
the p^p — > ir + + X case); for charged pions, they change sign (being somehow reduced in size) when 
going from polarized proton to polarized antiproton beams. The asymmetry for neutral pions is al- 
ways positive and almost half of that of the nearest charged pion case, in good agreement with isospin 
conservation in the fragmentation process. Information gathered on the pp dependence of the SSA's 
was less detailed. In Ref. [50J the pp dependence of the SSA for charged pions in the p^p — > 7r ± + X 
process was measured at moderate pp (0.2 < pp < 1.5 GeV/c) over a wide xf range, 0.2 < xf < 0.9. 
A threshold for the onset of the asymmetry was observed around pp ~ 0.5 GeV/c, below which the 
SSA was essentially zero and above which it increases in magnitude with pp. In general, in all data sets 
considered, the magnitude of the SSA at large positive xf was found to be more pronounced for the 
largest pp values reached. No data on the corresponding unpolarized cross sections were presented. 
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2) In the central rapidity region, \xf\ < 0.15 and 1.0 < pr < 4.5 GeV/c, the inclusive and semi-inclusive 
SSA's for neutral pions, using transversely polarized 200 GeV/c proton and antiproton beams off a hy- 
drogen target were measured |51j. These asymmetries were found to be consistent with zero within 
the experimental uncertainties over the full pp range covered, both for polarized proton and antiproton 
beams [see, e.g., Fig. [7] ( left panel), p. [45] . In this case, measurements of the corresponding invariant 
unpolarized cross sections were also provided. 

3) The SSA for inclusive direct-photon production, p^p — » 7 + X, was measured in the central ra- 
pidity region at 2.5 < pp < 3.1 GeV/c [52] • Data on the corresponding unpolarized cross section at 
2.5 < pt < 3.8 GeV/c were also provided. Only two data points were measured for the SSA, that was 
found to be consistent with zero within (large) experimental uncertainties. 

4) Some interesting results are also available for SSA's in inclusive A and 77-meson production [53| [54]. 
showing non- vanishing, sizable asymmetries in kinematical regions similar to those considered above. 

The E704 results triggered renewed interest on SSA's, given the larger cm. energy available as 
compared to that of previous experiments and the large, unambiguous effects measured. However, a 
clear separation between the xp and pt dependence of the results was not reached and the covered pt 
range was still limited, although the largest SSA's were measured at the highest pp. For all the 90s these 
results were the only ones available at sufficiently high cm. energies. The advent of new experimental 
setups started a new stage in the investigation of transverse spin effects in hadronic processes. 

Transverse SSA's were also measured, although again at lower total cm. energies, by experimental 
collaborations at IHEP (Protvino). In Ref. [55] they have shown the transverse SSA for 7r°'s inclusively 
produced in the collisions of a 40 GeV/c ir~ beam off transversely polarized protons and deuterons, in 
the region \xp\ < 0.2 and 1.6 < pp < 3.2 GeV/c. The SSA's were found to be approximately equal 
for the two reactions, compatible with zero or slightly positive for 1.0 < pp < 2.0 GeV/c, large in 
magnitude and negative for 2.4 < pp < 3.2 GeV/c. In Ref. [51] SSA's were measured for inclusive 
hadron production in collisions of a transversely polarized 40 GeV/c proton beam with an unpolarized 
liquid hydrogen target. The asymmetries were measured for ir^, K ± , protons and antiprotons, produced 
in the central region, 0.02 < xp < 0.10 and 0.7 < pp < 3.4 GeV/c. Within experimental errors, the 
measured asymmetries for tt^, and p show a linear dependence on xt = Ippj ^fs and change sign 
near xp = 0.37. For protons, a negative SSA, independent of xt, was found. 

2.1.2 BNL-RHIC results 

Starting from the end of the 90s, very interesting results on transverse SSA's and unpolarized cross 
sections for inclusive single particle production in polarized hadronic collisions came up and are currently 
collected by the experimental collaborations (STAR, PHENIX, BRAHMS) using the Relativistic Heavy 
Ion Collider (RHIC) at the Brookhaven National Laboratories (BNL). 

For the first time proton-proton collisions, with one (both) beam(s) transversely (longitudinally) 
polarized, can be studied at cm. energies as large as 200 GeV, with the possibility to reach -y/s = 500 
GeV in the near future. This represents a significant advancement with respect to all previous (fixed 
target) experiments. The setups of the three RHIC experiments allow to cover a large kinematical 
range, including particularly the previously almost unexplored negative xp region. 

Let us also remind a precursor experiment at the BNL Alternate Gradient Syncrotron (AGS). The 
E925 Collaboration [571 EE] measured transverse SSA's for charged pions using an incident 22-GeV/c 
transversely polarized proton beam on hydrogen and carbon targets, in the range 0.55 < xp < 0.75 and 
0.6 < Pt < 1-2 GeV/c. Although still at low energies and pt, this experiment gave a first, independent 
confirmation of the large SSA's measured by the E704 Collaboration. 

We now give a summary of the main high-energy RHIC results on SSA's: 
1) The STAR collaboration measured the transverse SSA for forward neutral pion production in p'p 
collisions at s/s = 200 GeV, average pseudorapidity (77) = 3.8 and 30 < E n ~ xp^fsjl < 55 GeV, 
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corresponding to 1.0 < (pp) < 2.4 GeV/c [59] . The unpolarized cross section was also measured 
at (r/) = 3.3, 3.8, 4.0 [SU ED] and found to be consistent with next-to-leading order (NLO) pQCD 
calculations [T7J EI]- In agreement with the E704 results, the SSA is small for xp values below about 
0.3 and becomes large and positive at higher xp. The spin-dependent azimuthal asymmetry for almost 
back-to-back dijets in p^p collisions at y/s = 200 GeV was measured over a wide pseudorapidity range 
for the two jets and found to be consistent with zero within experimental uncertainties [62] . 
Preliminary results for 7r° SSA's, based on a larger data sample, for both positive and negative xf 
in the range 0.2 < \xp\ < 0.6, collected at two average pseudorapidities, (77) = 3.3 (corresponding 
to 1.6 < (pr) < 3.7 GeV/c) and (77) = 3.7 (corresponding to 1.3 < (pp) < 2.8 GeV/c) have been 
presented [63]. They essentially confirm the previous data published in the positive xp region and show 
a SSA consistent with zero for negative xp [see Fig. 0] (right panel), p. 03J. A preliminary study of the pp 
dependence of the neutral pion SSA for xp > 0.4 and in xp-bins for 0.25 < xp < 0.56, using combined 
data at (77) = 3.3,3.7,4.0, in the range 1 < pp < 4 GeV/c was also performed [63]. Data show a non 
trivial dependence of the SSA on pp, which seems not to behave as a smooth decreasing function of pp 
[see Fig. |6j p. 04]. Other interesting STAR preliminary results concern the differential cross section for 
inclusive charged pion production at mid-rapidity and for inclusive neutral pion production at forward 
and mid-rapidity in pp collisions at y/s = 200 GeV [64] . 

2) The PHENIX Collaboration has measured the unpolarized cross section and the transverse SSA for 
neutral pions and charged hadrons produced at mid-rapidity (\r]\ < 0.35) in p^p collisions at y/s = 200 
GeV [651 EH1 E7]. For neutral pions (charged hadrons) the unpolarized cross section in the range 
0.5 < pp < 20 GeV/c (0.5 < pp < 7.0 GeV/c) was found to be in good agreement, above ~ 2 
GeV/c, with pQCD calculations at NLO accuracy [T71ET]. The SSA for both neutral pions and charged 
hadrons was found to be consistent with zero within errors in the range 1.0 < pp < 5.0 GeV/c [66J. 
These results are in agreement with those found by the E704 Collaboration at lower energy and in 
similar kinematical conditions [51] [see also Fig. [7J (right panel), p. 05], extending them to larger pp 
values, previously unexplored. The invariant unpolarized cross section for direct-photon production at 
mid-rapidity was also measured up to pp = 16 GeV/c in y/s = 200 GeV pp collisions [68j [69]. Data 
are well described by NLO theoretical predictions for pp > 5 GeV/c; in addition, the ratio of isolated 
photons to all nonhadronic decay photons is well described by pQCD for pp > 7 GeV/c. Other recent 
PHENIX preliminary results concern [70]: measurements of the SSA for forward neutrons and mid- 
rapidity J/ip's at y/s = 200 GeV; for the neutron case, large flat asymmetries for positive xp and SSA's 
consistent with zero in the negative xp range were observed. 

3) The BRAHMS Collaboration has recently presented preliminary results for invariant unpolarized 
cross sections for charged hadrons, ir^, K ± , proton, antiproton production in pp collisions at y/s = 200 
GeV measured at forward rapidities (y = 2.95 and 3.3) and pp up to 4 GeV/c [71]. Other interesting 
preliminary results have also been presented [72] concerning: the invariant cross section for — > n~ +X 
at y/s = 62.4 GeV and rapidity y ~ 2.7, 3.3; the SSA for 7r ± , K ± at y/s = 62.4 GeV, for both negative 
and positive xp in the range 0.25 < \xp\ < 0.6, and pp up to ~ 1.5 GeV/c; the SSA for tt , K , p and 
p at y/s = 200 GeV, 0.15 < xf < 0.3 and 0.2 < pp < 3.5 GeV/c. Among the main features of these 
results we mention: the almost flat x^-behaviour of An for at 200 GeV and 9 = 2.3° [see Fig. [5] (left 
panel), p. 01]; the comparable and sizable K + and K~ SSA's observed at 200 and 62 GeV [see Fig. [HI 
p. 06]; a large for antiproton production and an almost vanishing SSA for the proton case. 

2.2 Azimuthal asymmetries in polarized semi-inclusive deeply inelastic 
scattering 

Starting from the end of the 90s the HERMES Collaboration at DESY, and later the COMPASS Col- 
laboration at CERN and the CLAS Collaboration at Jlab, have performed a series of measurements of 
beam and target azimuthal asymmetries in semi-inclusive particle production in the deeply inelastic col- 
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lisions of longitudinally polarized leptons off both longitudinally and transversely polarized proton and 
deuteron targets. As it will be discussed in detail in the following sections, these azimuthal asymmetries 
with respect to the virtual photon direction, in the 7*-target cm. frame, can give a lot of information. 
The azimuthal angles under consideration refer e.g. to the transverse component of the observed hadron 
momentum (0/J and of the initial hadron polarization vector (0s) w.r.t. the virtual photon direction, 
measured from the leptonic plane [see Fig. [HJ p. [55]. By integrating the polarized cross section or the 
asymmetries, opportunely multiplied by some circular function of the azimuthal angles, one can dis- 
entangle several interesting contributions to single spin asymmetries and other polarized observables. 
These moments of the azimuthal asymmetries are usually indicated by A^g' , where Sb, St = U, L,T 
refer to unpolarized (U), longitudinally (L) and (for the target) transversely (T) polarized beam, target 
respectively, while W(4>) refers to some appropriate 0-dependent circular function, W — 1, sin0, cos0, 
sin 20, etc. As a phenomenologically relevant example, A s ™^ h±< ^ s ^ is given as [73] 

AS in{4> h ±4> s ) -o/ a - t , , - nv S d <Phd<p s sin(0 fe ± S ) [dajfe) - da{<f)s + tt)] , . 

a ut = A sir H0/i ± <Ps ) Fir^TTTTT^~ri 72 — i — vi • ^ z 

J d<p h d<p s [da{(f)s) + da((ps + 

Notice that in the literature sometimes the factor two in the above equation is not included. 

Throughout this section and in the sequel, we will make use of the well-known DIS kinematical 
variables: P, £, £', q = £ — £' are respectively the four-momenta of the target nucleon, the initial and 
final leptons and the virtual photon exchanged; Q 2 = —q 2 ; Xb = Q 2 /(2P ■ q) is the Bjorken variable; 
v = q ■ P/M is the lepton energy loss in the target nucleon rest frame; y = (P ■ q)/{P ■ £) is the so- 
called inelasticity; W is the cm. energy of the virtual photon-target nucleon system. Additionally, we 
also adopt the usual SIDIS variables: Ph is the four-momentum of the final observed hadron, and Ph± 
its transverse component w.r.t. the virtual photon direction of motion, with magnitude Ph± = |-P/u_|; 
z h = (P ■ Pf l )/(P ■ q) is the usual hadronic variable. 



2.2.1 HERMES results 

In its first stage, HERMES was equipped with a longitudinally polarized hydrogen target. This gave 
access to measurements of the quark helicity distributions (in DIS) and to the leading-twist azimuthal 
asymmetry A s ^ n ^ h , associated with the transverse momentum dependent (TMD) quark distribution 
in a longitudinally polarized hadron convoluted with the Collins function. In the virtual photon- 
target center of mass frame, even a longitudinally polarized (in the Lab. frame) proton target acquires 
a (kinematically) small (~ 1/Q) transverse spin component. Therefore, early attempts to measure 
transverse SSA's and some related effects (namely the Sivers and Collins effects) through the azimuthal 
asymmetry A s ^]f h were performed. Even if these transverse-spin contributions mix up with concurrent 
effects due to subleading-twist quark distribution and fragmentation functions, first indications for 
interesting physics were found. Concerning SSA's, however, the most interesting HERMES results 
came out with the advent of a transversely polarized hydrogen target. This allowed to extract the first 
unambiguous indications of non-vanishing leading-twist transverse SSA's in SIDIS, and to disentangle 
the possible contributions due to the Sivers and Collins effects. 

The kinematical regime covered by the HERMES experiment is approximately the following (some 
differences are possible for specific measurements; we refer to the original papers for more detail): 
1 < Q 2 < 15 (GeV/c) 2 ; 0.023 < x B < 0.4; y < 0.85; W 2 > 10 GeV 2 ; 2 < E h < 15 GeV, 0.2 < z h < 0.7; 
Ph± > 50 MeV/c. Notice that in some cases more selective cuts on the last two variables are used. 

We now briefly summarize the main results of the HERMES experiment, quoting only those results 
that are relevant in the context of this review paper (a full account of all HERMES results may be 
found in the papers quoted in the bibliography): 

1) First evidence of azimuthal SSA's for charged [71] and neutral [75] pion production in deeply inelastic 
scattering of polarized positrons off a longitudinally polarized hydrogen target; A s ^ h was found to be 
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significant and comparable for tt + and tt° production, and consistent with zero for ir~ . Indications of 
an increasing of this analyzing power with increasing of xb, and with Ph± up to ~ 0.8 GeV/c, were also 
found. A s ^ l2 ^ h was found to be consistent with zero for both ir + and 7r~. 

2) Measurement of azimuthal SSA's in SIDIS production of pions and kaons on a longitudinally polar- 
ized deuterium target [76]; the dependences of these asymmetries on i^, Ph±, and Zh were investigated; 
positive asymmetries for tt + and tc° and an indication of a positive asymmetry for 7r~ mesons were 
found; the asymmetry for K + was measured to be comparable with that for 7r + mesons; 

3) Measurement of the azimuthal SSA's A s ^p^ h± ^ s ^ in SIDIS production of charged pions and kaons on 
a transversely polarized hydrogen target [771 178]. For the first time two different phenomena, possibly 
associated with the Sivers and Collins effects, indistinguishable in previous data were disentangled. The 
(sin(0/j — (ps)) moment was measured to be positive and nonzero for 7r + and consistent with zero for 7r~ 
mesons [see Fig. [16] (left panel), p. [58]. The (sin(0^ + 4>s)) moment was found to be positive for ir + and 
negative and comparable in magnitude for ir~ mesons [see Fig. [T7J (left panel), p. EH]. For charged kaons 
the measured A^^ h ~^ 3 ' asymmetry came out to be positive, showing large values for K + ; A s J^^ h+( ^ s ^ 
for K~, even if with large error bars, presented large positive values [see Fig. [TS], p. I5T?] . 

4) Extraction of subleading-twist effects in SSA's in SIDIS production of charged pions on a longitudi- 
nally polarized hydrogen target [79] . By combining measurements of azimuthal asymmetries with both 
longitudinally and transversely polarized hydrogen targets (see previous point) this contribution was 
found to be significantly positive for ir + mesons, therefore dominating the asymmetries on a longitudi- 
nally polarized target (see points 1) and 2) above). The subleading-twist contribution for ir~ mesons 
was found to be small. 

5) Measurement of the subleading-twist beam spin asymmetry A s ^j in the azimuthal distribution of 
pions produced in SIDIS [80]. A s ^f h for 7r + 's was found to be consistent with zero for low Zh and rising 
with Zh up to about 0.02; for 7r~'s, A s ^f' h was found to be consistent with zero within the statistical 
accuracy; for neutral pions, the asymmetry was found to be roughly constant and of the order of 0.03, 
decreasing to zero only in the lowest /highest Zh regions. 

6) Interesting preliminary results on azimuthal SSA's in SIDIS production of charged kaons on a trans- 
versely polarized hydrogen target have been also presented (see, e.g., Ref. [HT]). 



2.2.2 COMPASS results 

The COMPASS Collaboration at CERN has performed measurements of transverse azimuthal asym- 
metries for charged hadrons produced in deeply inelastic scattering of the CERN-SPS 160 GeV/c muon 
beam off a transversely polarized deuteron target. The high-energy lepton beam allows to reach large 
values of Q 2 and a very large coverage of W. This can be very helpful in disentangling genuine TMD 
effects from perturbatively generated contributions. The large xb region covered, well below the valence 
region, allows to better estimate the first moments of the TMD distributions and test related sum rules. 
Moreover the use of a transversely polarized (isospin scalar) deuteron target well complements with 
HERMES results, in particular concerning flavour separation of the TMD distributions and fragmen- 
tation functions. COMPASS also plans to perform measurements with a transversely polarized proton 
target in the near future. 

The asymmetries A S j^p^ h±( ^ s ^ for charged hadrons, and their dependence on and, for leading 
hadrons, z^, Ph±, were measured and found to be compatible with zero within statistical errors [82~j 
[see Fig. [TH] (right panel), p. EH] and Fig. [T7J (right panel), p. EH]. The kinematical region covered was: 
Q 2 > 1 (GeV/c) 2 , W > 5 GeV, 0.003 < x B < 0.3, 0.1 < y < 0.9. Concerning the observed hadron, the 
cuts Ph± > 0.1 GeV/c, Zh > 0.2 (zh > 0.25) have been imposed for all (leading) hadrons. 

Recently, new high-precision measurements for the same observables and in the same kinemati- 
cal region but with larger statistical significance have been published [83]. These results extend and 
essentially confirm those of Ref. [82] , with statistical errors a factor of two smaller. 
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Interesting preliminary results have also been presented on A^- T for charged pions and kaons, in 
the same kinematical configuration and using data collected in 2003-2004 [HJE5]. Again, all transverse 
spin effects investigated are found to be compatible with zero. 

Very recently COMPASS has also presented preliminary results [86], in the same kinematical region, 
on the first measurement of the other six possible transverse azimuthal asymmetries (see, e.g., Ref. |87j). 

A sin(3^ s)) A *»{*K-*s) ( leading twist ) ) and A ffi S) A M2*h-*s)^ A cos0 S) A cos(2^-0 S ) (subleading 

All these asymmetries are compatible with zero within statistical errors. 
2.2.3 JLAB-CLAS results 

The JLAB-CLAS Collaboration has performed measurements of both longitudinal beam and target 
spin asymmetries for pion production in SIDIS off protons using an electron beam with energy 4.3 GeV 
and 5.7 GeV. More precisely, they have measured: 

1) The sin 0^ moment of the subleading beam spin asymmetry, A L i/ >h , for the process ep — > e'ii + + X 
above the baryon resonance region, in the range 0.15 < xb < 0.4, 0.5 < zy t < 0.8, using a polarized 
electron beam of 4.3 GeV [88] and of 5.7 GeV [89] (in this last case only preliminary results are 
available). The measured beam SSA is positive for a positive electron helicity, varying roughly in the 
range 0.02 -j- 0.08; it increases with and shows no significant dependence on the x^-range; 

2) The sin 0ft, sin 20ft moments of the longitudinal target spin asymmetry, A^f h and A S ^ L , for both 
charged and neutral pions, in the range 0.15 < xb < 0.4, 0.5 < Zh < 0.8, using an electron beam of 5.7 
GeV (preliminary results, see Refs. [891 [90]). Data for 7r +, s show clear sin 0ft and sin 20ft modulations, 
leading to (sin 0ft) = 0.058 ± 0.011 and (sin 20ft) = —0.041 ± 0.011 (errors are statistical). 

2.3 Azimuthal asymmetries in unpolarized processes 

Several azimuthal asymmetries in unpolarized cross sections for SIDIS, Drell-Yan and e + e~ processes 
have been measured. Among other explanations of their origin, intrinsic parton motion is one of the most 
accredited. Therefore, these asymmetries are intimately related to transverse single spin asymmetries 
and in particular to TMD parton distribution and fragmentation functions. For this reason, we find 
useful here to summarize some of the most relevant results. 

2.3.1 SIDIS processes 

The most general azimuthal dependence of the cross section for SIDIS off an unpolarized target can be 
written as (see, e.g., Ref. [91] and references therein): 

— oc A + B cos fa + C cos 20/j + D sin <fi h + E sin 2 fa , (3) 

dfa 

where the angle 0ft, as mentioned above, is the azimuthal angle of the observed hadron momentum 
around the virtual photon direction, measured starting from the leptonic plane [see Fig. [TH p. 155] . 

In collinear pQCD, A is of order 0(a®); the first corrections to A and the B-E contributions come 
from higher order terms, in particular B, C are of order 0(a s ), D, E are of order 0(a 2 s ). Further- 
more, the T-odd contributions D and E are present only for parity violating weak interactions or for 
electromagnetic interactions with longitudinally polarized lepton beams. Nonperturbative transverse 
momentum effects and TMD functions may give leading-twist contributions to B and D terms which 
may be non- negligible even at very large Q 2 . 

The EMC Collaboration at CERN measured the azimuthal distribution [see Fig. [151 P- ES] and its 
(cos 0^), (cos 20ft), (sin 0ft) moments for charged hadron production in SIDIS with a muon beam of 
energy 280 GeV off a liquid hydrogen target [92]. Both the xf and (for the moments) Ph± dependences 
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were measured. The main kinematical cuts considered were: Q 2 > 4 (GeV/c) 2 , 40 < W 2 < 450 GeV 2 , 
v > 20 GeV, y < 0.8, P^± > 0.2 GeV/c. The main results found were: A strong dependence of the 
(cos 0^) moment on xp, which has a small and positive value for negative Xp, becomes negative at 
xf ~ and falls below —0.1 for xf > 0.5. A small dependence of the (cos^), (cos 20^) moments on 
Ph±. The values of (cos20/ l ), (sin </>/,,) as functions of a?j? were found to be close to zero. Previous results 
limited to the forward hemisphere were also available [93] . 

The Fermilab E665 Collaboration measured the azimuthal asymmetry of forward charged hadrons 
produced in SIDIS with a 490 GeV muon beam off proton and deuteron targets [94]. The main kine- 
matical cuts imposed were: Q 2 > 3.0 (GeV/c) 2 , 100 < W 2 < 900 GeV 2 , 60 < v < 500 GeV, x B > 0.003, 
0.1 < y < 0.85, Ph > 8 GeV/c. They measured, among other observables, the (cos0h) moment as a 
function of the minimum hadron transverse momentum up to 2.5 GeV/c, imposing W 2 > 300 GeV 2 
and Xp > 0.2. Moment values are negative, starting from about —0.04 at the lowest cutoff in P h ±, 
decreasing up to about —0.12 for a P^ cutoff of ~ 1.5 GeV/c and then increasing again (but with 
larger errors) to ~ —0.03 for the largest considered cutoff in P^ of 2.5 GeV/c. 

EMC and E665 results are consistent with theoretical models including parton transverse momenta 
arising both from intrinsic nonperturbative motions and NLO pQCD radiative corrections. Apparently, 
the nonperturbative contributions dominate in the kinematical regimes covered by these experiments. 

The ZEUS Collaboration at HERA-DESY measured the azimuthal angle distribution for charged 
hadrons produced in neutral current deeply inelastic positron-proton scattering [55] . The dependence of 
the moments of the distribution, (cos <frh), (cos 2(f) h) on (the minimum value of) the transverse momentum 
of the charged hadrons observed, P™ n , was also studied. At low P/^" 1 a clear (cos<fth) contribution was 
observed, while as the value of P^j" 1 increased, a (cos 2(f>h) term became evident; (sin <ph) was found to be 
consistent with zero independently of the P^ n value chosen; ii) For hadrons produced at large transverse 
momenta the (cos<frh) moment was measured to be negative, in agreement with pQCD predictions; the 
(cos20/j) moment was measured for the first time and found to be non-zero and positive, with increasing 
magnitude as a function of P™ n . These results were interpreted as giving clear evidence for a pQCD 
contribution to the azimuthal asymmetry, in contrast with the results of EMC and E665 Collaborations 
in different kinematical regimes. 

The ZEUS Collaboration has also performed a study of the azimuthal asymmetry for inclusive jet 
production in neutral current deeply inelastic e + p scattering. Again, NLO pQCD calculations seem to 
give a good description of the observed results [96] . 

More recently, the ZEUS Collaboration presented new results on the dependence of azimuthal mo- 
ments on the pseudorapidity and the minimum transverse energy of the final-state hadrons [97] . Both 
neutral and charged hadrons were considered over an extended phase space as compared to Ref. [95J. 
The value of (cos^) was found to be negative for t] < —2, becoming positive for larger rj. NLO pQCD 
predictions reproduce reasonably well this behaviour but fail in describing the magnitude of the asym- 
metries. However, the predicted values of (cos 20^) agree with the data. A deviation of (sin0/j) from 
zero at the level of three standard deviations was observed, while (sm2(p h ) was consistent with zero. 

2.3.2 Drell-Yan processes 

Concerning the Drell-Yan (DY) process, it can be shown that after averaging over the initial hadron 
polarizations, summing over the final lepton spins, and applying invariance principles (permutation 
symmetry, gauge invariance, parity conservation and unitarity), the most general form of the lepton 
angular distribution, in the lepton-pair cm. frame reads 

— : — — oc 1 + A cos 9 + a sin 20 cos 0+ (z//2)sin 2 # COS/0, 4 

a a cos 0d<p 

where 9 and <fi are respectively the polar and azimuthal angle identifying the lepton-pair direction of 
motion [see also Fig. [TUJ p. [50]; the coefficients A, /i, v are in general functions of the other relevant 
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kinematical variables: the total cm. energy y/s, the lepton-pair invariant mass M, its transverse mo- 
mentum qp, and its fraction of longitudinal momentum Xp. The values of these coefficients depend on 
the choice of the reference frame, the so-called Collins-Soper (CS) frame being often adopted (other 
popular choices are the Gottfried- Jackson (GJ) frame and the w-channel frame). When the angular 
distributions in both cos# and <fi are measured in a given frame, one can easily compute the coefficients 
in any other one. Notice however that this is not the case when only the cos# distribution is measured. 

In the naive parton model for the Drell-Yan process, the assumption of massless quarks implies that 
the virtual photon is transversely polarized, so that A = 1, n — v — 0, and one has da/dcos9 ~ 1+cos 2 9. 
A general relation by Lam and Tung [981 EH1 HOP] - analogous to the Callan-Gross relation in DIS, is 
expected to hold in any reference frame: A = 1 — 2u. This relation is not modified by first-order QCD 
corrections [lOOj . but can be influenced by parton intrinsic motions and other nonperturbative effects. 

Concerning experimental results, we limit ourselves to mention some measurements that, as will be 
discussed in section [5J are relevant in the context of this review. 

The NA3 pm] and NA10 [T02l[T03] Collaborations at CERN have measured the angular distributions 
of high-mass muon pairs produced by 140 GeV/c and 194 GeV/c ir~ beams impinging on a tungsten 
target, and by a 286 GeV/c n~ beam on deuterium and tungsten targets. The dependence of the 
parameters A, fi, v on the cm. energy y/s, the invariant mass M, transverse momentum qp and rapidity 
y of the lepton pair, and the fractional momentum of the parton in the pion, x\, was measured. The 
kinematical region covered was indicatively (see Ref. [103] for more details): 4.0 < M < 8.5 GeV/c 2 
and M > 11 GeV/c 2 , 0.32 < q T < 2.60 GeV/c, -0.14 < y < 0.71, 0.28 < x x < 0.80. No evidence for a 
cm. energy dependence or a nuclear dependence of the angular distribution parameters was found. The 
value of A was found close to 1, as predicted by the naive parton model, and substantially independent 
of any kinematical variable, apart from the high X\ region, where indications of a decrease of A with x\, 
as suggested by higher-twist effects [1041 1105j . were found. The value of fi was found close to zero in 
the CS frame, as expected if both annihilating partons contribute equally to the transverse momentum 
of the lepton pair. The most important result of this analysis was the strong dependence of v on qp, in 
clear disagreement with pQCD expectations. In particular, starting from v ~ at the lowest qp values, 
v was found to grow up to 0.15 -j- 0.30 (depending on y/s) at larger qp values [see Fig. CD] (left panel), 
p. [51] . Since A ~ 1, this result also implied a clear violation of the Lam- Tung relation. 

The E615 Collaboration at Fermilab has measured the cross section, the transverse momentum 
distribution and the full angular distribution in cos# and as a function of M, xp, and qp, for muon 
pairs produced by a 252 GeV n~ beam interacting with a tungsten target [106] . The dependence of 
the parameters A, fi, v, on the fractional momentum of the parton in the pion, Xi, was also measured. 
The kinematical region covered was: 4.05 < M < 8.55 GeV/c 2 , < q T < 5.0 GeV/c, 0.2 < x x < 1.0. 
Again, a clear violation of the Lam- Tung relation, and a strong dependence of v on qp were found. In 
this case also, starting from v ~ at the lowest qp values, v was found to grow up to ~ 0.54(0.73) in 
the GJ(CS) frame at the largest qp values. 

Recently, the FNAL E866/NuSea Collaboration has presented results on the angular distributions of 
DY dimuons produced using an 800 GeV/c proton beam on a deuterium target [107] . in the kinematical 
range 4.5 < M < 15.0 GeV/c 2 , qp < 4 GeV/c, and < xf < 0.8. No significant cos 20 dependence 
was found, in contrast with pion-induced DY data discussed above. These results put constraints on 
theoretical models that predict a large cos 20 dependence originating from QCD vacuum effects and 
suggest that the sea-quark Boer-Mulders functions [108] are much smaller than those for valence quarks. 

2.3.3 e + e~ — > ixix + X processes 

Recently the Belle Collaboration at the KEKB asymmetric-energy e + e~ storage rings has published [109] 
very interesting results on the inclusive production of hadron pairs (namely charged pion pairs) in e + e~ 
annihilation. Using two different reconstruction methods, they have found evidence of statistically 
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significant azimuthal asymmetries. The first method is based on the reconstruction of the thrust axis 
in the e + e~ cm. frame: combining two hadrons from different hemispheres in jet- like events, with 
azimuthal angles <pi and 02 defined w.r.t. the plane spanned by the lepton momenta and the thrust 
axis [see Fig. [12] (left), p. [52], a cos(0i + <p2) modulation of the dihadron yield has been observed [see 
Fig. [13] (left panel), p. [53]. An alternative method does not require the knowledge of the thrust axis. 
In this case the dihadron yields are measured as a function of one angle, (po, that is the angle between 
the plane spanned by the momentum vector of the first hadron and the lepton momenta and the plane 
defined by the two hadron momenta [see Fig. [12] (right), p. [52]. In this case a cos(2^o) modulation 
appears [see Fig. [13] (right panel), p. [53]. As it will be discussed in section 15.41 these results are related 
and give access to the product of two Collins functions. 

The kinematical cuts imposed in both methods are the following: —0.6 < cos(6 l i a b) < 0.9, where #i a b 
is the polar angle in the laboratory frame; z\p > 0.2, where z\ = 2EhJQ, and Q is the cm. energy; 
Qt < 3.5 GeV/c, where Qt is the transverse momentum of the virtual photon in the rest frame of the 
hadron pair. Moreover, to cancel acceptance effects double ratios of normalized rates for unlike-sign (U) 
over like-sign (L) pion pairs (UL double ratios) have been considered. At large Z\,i values the measured 
asymmetries are significantly different from zero [see Fig. [13] p. [53]. In Ref. [110] preliminary results 
on the double ratios of unlike-sign, U, over all charged (C) pion pairs have been presented. These UC 
double ratios can give useful additional information on the Collins function. Belle results have been 
confirmed by recent preliminary data with enlarged statistics and much smaller error bars [11 lj . 

3 Theoretical approaches to single spin asymmetries 

We present here some of the approaches developed to explain the interesting features of the SSA data 
discussed so far. We will try to follow a historical perspective, starting with the former attempts to 
describe the SSA's observed in hadronic collisions and then entering into the latest developments. 

The first approach we address is based on a pQCD factorization scheme with inclusion of spin and 
intrinsic transverse momentum dependent (TMD) effects: we will generally refer to it as the TMD 
approach. In particular we start with its application to inclusive hadron production in hadron-hadron 
collisions and its formulation as a phenomenological generalization of the parton model. This form of 
the TMD approach will be also referred to as the Generalized Parton Model (GPM) and will find a 
more detailed description within the helicity formalism in section [4] We then discuss more recent, and 
significant, QCD developments of the TMD approach with inclusion of initial and final state colour 
interactions which modify, to some extent, the partonic interpretation of the GPM approach. Some 
details on the impact parameter picture for SSA's, which gives a useful physical interpretation of the 
TMD approach, are also presented. 

The second approach we consider extends the QCD collinear factorization theorems to higher-twist 
contributions: the so-called twist-three approach. We recall its basic ingredients, giving some details of 
its application to SSA's in pp — > it + X and its later extension to other processes. 

The issues of validity, applicability and possible overlapping of these two approaches, for different 
kinematical regions and processes, as well as the status of factorization in pQCD, are also discussed. In 
particular, we will see how the TMD approach applies naturally to leading twist asymmetries, like SSA's 
in Drell-Yan processes and SIDIS, or azimuthal asymmetries in e + e~ — > h\hi + X, at low transverse 
momentum, where also a second large scale (the virtuality of the exchanged boson) is present. With 
some caution, it can be also adopted in double inclusive production of hadrons (jets, or photon-jet) in 
hadronic collisions. 

On the other hand, the collinear higher-twist approach applies in principle to subleading asymme- 
tries, like SSA's in pp — > h + X and pp — > 7 + X, where a single hard scale (the transverse momentum 
of the observed particle) appears. For these cases the GPM can be viewed as an effective and, as we 
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will see, phenomenologically successful, TMD description. 

For completeness we will give some details of an alternative description of SSA's based on a semi- 
classical model: the so-called orbiting valence quark model. 

Among other attempts to describe SSA's, even if on a more qualitative ground, we mention: i) 
a model discussed in Ref. |112j . for pp — > ir + X, and in Ref. [113j . for SIDIS, where the quark 
chirality flip and the required T-odd phase (see below) are associated with instanton fluctuations in 
the QCD vacuum; ii) a soft rescattering mechanism [114] with a Pauli coupling that leads in SIDIS to 
a sin((j)h + <ps) dependence (like that usually ascribed to the Collins effect in the TMD approach); Hi) 
a coherence mechanism [115j . where the entire hadron contributes to the scattering process, aimed to 
describe the sizable values of observed in pp collisions at large xf- 

For technicalities and a comprehensive discussion of the topics addressed here we refer the reader to 
the original bibliography. 

3.1 Parton model approaches including intrinsic transverse motion 

It has been known for a long time, especially in the fixed target regime, that unpolarized cross sections 
for inclusive particle production in high-energy hadron-hadron collisions at moderate pt (few GeV/c) 
cannot be properly described in the collinear factorized pQCD approach, neither at LO nor at NLO 
accuracy [116] . Only recently, improved NLO calculations including threshold resummation effects have 
been developed, see e.g. Ref. [117] . These results show how for inclusive cross sections integrated over 
all hadron rapidity range the agreement with data is much less problematic. On the other hand the 
rapidity-dependent case is still under study. 

Even more dramatic is the failure of the description of the Drell-Yan low qr spectrum, where qr is 
the lepton-pair transverse momentum. Indeed, in this case no transverse momentum can be generated 
in the collinear LO approximation. All these facts have raised, starting from former work of Feynman 
and collaborators [1181 1119] , the interest in considering the role of intrinsic transverse momentum. 

In the context of transverse SSA's the intrinsic transverse motion of partons has played an even more 
essential role. Indeed, as shown in Ref. [16], in collinear pQCD An appears only as the immaginary 
part of interference terms between spin-flip and no-spin-flip partonic scattering amplitudes. Since at LO 
these are real and helicity is conserved for massless partons, it was natural to expect A^ ~ a s mj ' \fs. 

Contrary to these expectations several experimental observations, as seen in section [21 show sizable 
SSA's in the high-energy regime. 

The cross section for a generic inclusive process A B — > C + X, in a phenomenological approach 
based on the generalization of the factorization theorem with the inclusion of intrinsic motions, k±, of 
partons inside a nucleon and of final hadrons relatively to the fragmenting parton, reads, schematically: 

da oc fa/A(x a , k± a ) <g> f b / B (x b , k ±b ) <g> da ab ~* cd (x a , x b , k ±a , k ±b ) g> D c / C (z, k ±c ) , (5) 

a,b,c 

where (g) stands for appropriate convolutions both in the light-cone momentum fractions, x a ^, z, and 
in the fej_'s. The /'s and the -D's are the TMD parton distribution (PDF) and fragmentation functions 
(FF), respectively; k± a)b is the transverse momentum of parton a, b w.r.t. the hadron A, B momentum 
and k±c the transverse momentum of hadron C w.r.t. the fragmenting parton direction of motion. 

The above QCD factorization scheme - with unintegrated k± dependent distribution and fragmen- 
tation functions - has never been formally proven for the process under consideration, but only for the 
Drell-Yan process [120j and for two-particle inclusive production in e + e~ annihilation [121] (somehow 
a time-reversed Drell-Yan process). Factorization at low transverse momentum for SIDIS has been 
recently addressed [1221 1123] and proved at the same level as for the two above mentioned processes. 
Recent developments on the status of the more problematic fc^-factorization in AB — ► CD + X will be 
discussed in section 13.1.11 
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The study of spin asymmetries requires the extension of Eq. ([5]) to the polarized case. The inclusion 
of k± effects and TMD functions can in principle lead to sizable SSA's; however, if the TMD functions are 
symmetric in the azimuthal direction, would be still negligible. It was then originally suggested by 
Sivers [T8l [19] that there could exist a correlation between the azimuthal distribution of an unpolarized 
parton and the spin of its parent hadron. That is, we could have the following nonvanishing asymmetry: 

A/ a /At (X, k±) = f a /A\ (X, k±) - f a/A l (x, fej.) = f a/A i (x, k±) - f a/A \ (x, -k±) , (6) 

where the last relation comes from rotational invariance and T (!) stands for the upwards (downwards) 
transverse spin direction of the hadron A. It is easy to see that Af a / A t must be, by parity conservation, 
proportional to (p A x k±) ■ S A , where p A , S A are the momentum and the spin of the hadron A, and 
k± the transverse momentum of the parton. This function is fc^-odd and chiral-even, as will be clear 
in section HI It is also naively time reversal odd (T-odd). Notice that the above asymmetry can be 
rephrased in terms of the number density of unpolarized partons inside a transversely polarized hadron 
(Pa, k± being unit vectors and k± = \k±\): 

/a/At (X, k±) = f a /A(x, fcj.) + - A N f a/A t (x, k±) [p A X fcjj • S A , (7) 

where A N f a/Al (related to [108] . see Ref. [73]) is referred to as the Sivers function. 

In this first formulation the numerator oi Am, see Eq. (JTJ), assumed to be generated by the Sivers 
effect, was computed according to the following expression: 

dAa S[vCTS oc J2 A A/at (x a , fcia) ® fb/ B {x h ) ® dtf*-**^ x b: k ±a ) <g> D c/c (z) ; (8) 

a,b,c 

that is, by keeping TMD effects only where the collinear approximation would otherwise give zero 
(remember that the Sivers function vanishes for k± = and is fcj_-odd). 

In 1993 Collins proposed a proof of the vanishing of the Sivers function [20]. By imposing, in 
the light-cone gauge (this point would have become crucial), parity and time reversal invariance, he 
showed that strong interactions forbid such an asymmetry. However, on the basis of possible initial 
state interactions between the two colliding hadrons, in Ref. [124] it was argued that this asymmetry 
might be allowed. By a reasonable parameterization of this function a good description of E704 data 
was then obtained [1241 1125] . 

In Refs. [20] |21~] another possible asymmetric azimuthal angular dependence was proposed: a left- 
right asymmetry in the fragmentation of a transversely polarized quark into a spinless (or, more gener- 
ally, unpolarized) hadron. This implies a new TMD and spin dependent twist-two function (p q and s q 
being the momentum and the spin of the quark and pc — zp q + k± the hadron momentum), 

D c / q i (z, k±) = D c /g{z, k±) + - A N D c/g i (x, k±) {pg x fcjj ■ Sg , (9) 
where A N D c / g ^ (or in Refs. |126[ I1U8] , see also Ref. [73]) is the Collins function. Again we have 

AD c /qi (z, k ±) = D c / q i (z, k±) - D c/q i (z, k±) = D c/q i (z, k±) - D cjq[ (z, -k±) , (10) 

where T (!) stands for the upwards (downwards) transverse spin direction of the fragmenting quark. In 
this case, due to the rescattering of the hadron with the remnants of the jet (or, in more formal words, 
since the final hadron, contrary to the incoming one, cannot be represented as a plane wave state), 
time-reversal invariance does not imply any constraint, allowing the existence of such an asymmetry. 

This new effect in the fragmentation mechanism was then proposed as a transverse quark polarimeter 
with the aim of extracting information on the transverse polarization of a quark, 1.6. , clS 8b tool to measure 
the transversity distribution. This function, denoted by hi q (x) (or A T q, Sq), and defined as 

hi q (x) = f q y p i(x) - f q i/ p t(x) , (11) 
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is the difference between the probabilities to find a quark polarized along the transverse proton po- 
larization and against it. This chiral-odd function, in contrast with the unpolarized and the helicity 
distribution functions (both chiral-even) , is much more difficult to measure. Indeed, it cannot be ac- 
cessed in DIS, since it implies a helicity flip at the parton level whereas QED and QCD for massless 
quarks conserve helicity. It has therefore to be accompanied by another chiral-odd quantity. As orig- 
inally proposed by Ralston and Soper in Ref. |127j . the most promising source of information on h\ q , 
is the double transverse spin asymmetry, Att, in Drell-Yan processes, where it appears quadratically, 
but unfortunately no data are still available. 

Collins idea was to consider a new process where another chiral-odd function could appear: this is 
the SIDIS process with a transversely polarized proton target, ep^ — > e'h + X, where the azimuthal 
dependence of the final hadron distribution, w.r.t. the lepton scattering plane, would probe at the same 
time the transversity distribution function, the transverse spin dependence of the hard scattering and, 
eventually, the spin dependence in the fragmentation. This process, as we will discuss in section has 
become, indeed, the first source of phenomenological information on the transversity distribution. 

A calculation of the Collins function, based on a sigma model for pions with an effective Lagrangian 
incorporating chiral symmetry breaking, showed that AD c / q t(z, k±) could be different from zero [20J. 

This idea was then extended to the study of in p^p — > ir+ X collisions. In the factorized approach 
discussed above, the numerator of the SSA generated by the Collins effect would read 

dAa Colhas cx ]T h la (x a ) ® f b/p (x b ) ® dAa ab ^ cd {x a , x b , k ± ) ® A^/d (z, k ± ) , (12) 

a,b,c 

where dAa ab ^ cd is the partonic double transverse spin asymmetry related to the partonic spin transfer, 
Dnn = dAa/da. Notice once again that only the leading fcj_-dependence is kept. 

A first phenomenological attempt to explain the SSA in terms of such effect was presented in 
Ref. [128] . The azimuthal asymmetry in the fragmentation process was generated via the Lund string 
mechanism |129j . According to this model, in the breaking of the string spanned between the scattered 
quark and the target, the quark and the antiquark of every pair created acquire a transverse momen- 
tum (w.r.t. the direction of the string) and a polarization which are strongly correlated. Namely, the 
antiquark polarization goes like P q ~ —z x k± where z is the unit vector along the string direction. 

In order to form a pion, the polarized scattered quark (coming from the hard scattering) has to 
combine with the antiquark (belonging to the string pair) into a spin singlet. This implies a cor- 
relation between the q and q polarizations, resulting in an asymmetric azimuthal distribution of the 
pion w.r.t. the quark polarization. In this model the Collins effect is introduced only for leading, or 
"first-rank", pions (those containing the original quark spanning the string): the azimuthal distribu- 
tion of higher-rank (subleading) pions is assumed to be symmetric. It is worth to mention that even 
starting from Eq. (fT2"j) . that is a pQCD inspired formula, the authors of Ref. |128j adopted a partonic 
spin-transfer D^n equal to one (as it is for very small scattering angles) and a simple power-like pt- 
dependent partonic cross-section based on the observed inclusive pion spectra. By assuming that the 
transverse polarization of the u and d quarks in the proton are +1 and — 1 at x = 1, respectively, they 
obtained qualitative agreement with the E704 data, although underestimating their size. 

Among the peculiar features of this model we recall that for spin-1 meson (e.g. the p) production the 
spin-fcj_ correlation gives a very different Collins effect; as a consequence, A^{p) ~ — (1/3)^4^(71") [130] . 
No SSA is expected for K~ production (which has no valence quark in common with the polarized 
proton) since the subleading fragmentation mechanism is assumed to be symmetric. 

A more detailed phenomenological study of the Collins effect in pp collisions, based on the factorized 
scheme of Eq. ffT2"j) . was performed in Ref. [131| . where the authors adopted the helicity formalism with 
inclusion of TMD effects. We postpone the details of this approach to section HI Here we only recall 
that by a reasonable parameterization of the Collins function and the transversity distribution a good 
description of the E704 data was obtained. As we will discuss in section [5] a complete treatment of the 
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kinematics, taking into account the full k± dependences in the factorized formula of Eq. ffl2l) . plays a 
big role and leads to a strong revision of this conclusion. More generally we will see how the correct k± 
kinematics is crucial in the TMD approach to inclusive particle production. 

Almost in the same years, a first systematic classification of the twist-two and twist-three trans- 
verse momentum and spin dependent parton distribution and fragmentation functions was presented 
in Refs. [1261 1108j . For the twist-three case these functions have no partonic interpretation since they 
are related to quark-quark-gluon matrix elements. Further improvements beyond leading twist in the 
above classification were reached in Ref. |132j . This work has been recently reassessed and completed 
in Ref. [87J where some new developments (see also Refs. [1331 1134j ) have been taken into account; 
in particular, a systematic presentation of the relevant observables for all combinations of lepton and 
hadron polarizations at small transverse momentum and twist-three accuracy has been given. For space 
reasons, and given its relevance, in the following we will restrict to the leading-twist TMD sector. We 
refer the reader to the above mentioned papers for a more detailed discussion. 

The starting point in this classification is the correlator of quark fields (i/j) for a polarized nucleon 
(with spin S) entering the DIS processes in the diagrammatic expansion of the hard scattering amplitude, 

^(p,s,k) = J SLe ihy (P, s\^(o)Uy)\P, S), (13) 

and the analogous correlator for the fragmentation sector 

A ij (P h ,S h ,k) = J2j ^e ik - y (0\^ i (y)\P h ,S h ,X)(P h ,S h ,X\^ j (0)\0). (14) 

By integrating over k~ = (k° — k 3 )/^/2, where k^ = xP^ + k± and kj_ = (0, k±, 0), and keeping the 
leading twist terms, one gets, for the distribution sector, 



$(x,k ± ,S) 



f\V>+ + /it— — + \Sl9il + M 9it I 7 fi+ 



+ hvr i(T^rn%S T + I S L h 1L + — h 1T I — + h x ^ M 



, (15) 



where n± are auxiliary light-like vectors and M is the nucleon mass. All functions above depend on x 
and Few words on the above notation are helpful: /, g and h stand for unpolarized, longitudinally 
polarized, and transversely polarized quarks, respectively: the subscript "1" stands for leading twist; 
the apex _L indicates the explicit presence of transverse momenta with a noncontracted index [126J; 
the subscripts L, T stand for the longitudinal and transverse polarization of the hadron. Therefore, for 
instance, f\ is the TMD unpolarized PDF, is the Sivers function and hiT, h^ T are related to the 
TMD transversity distribution. For relations among different notations see also Refs. [73l [121 1135] . 

By appropriate Dirac projections, $' r ' = Tr(r$), one can single out the various sectors of distri- 
bution functions. In particular, V = {n-) a ^ a /2 projects out the fi sector, and T = |io"^(n_) /i ^^-7 5 
gives the hi sector. Among the others we only mention, in particular, the function hj^, also known 
as the Boer-Mulders function, which has been proposed in Ref. [136j as a possible tool to explain the 
azimuthal asymmetry observed in unpolarized Drell-Yan processes. 

The partonic meaning of the eight functions appearing in Eq. (Tl5l) and their relations with those 
defined in a complementary formalism, the helicity approach, will be discussed in more detail in section HI 

The same structure can be obtained for the fragmentation sector, starting from the corresponding 
correlator, Eq. ( Tl4l) . In this case one gets, among the others, the Collins function, A N D c / q i , and a 
chiral-even TMD fragmentation function of an unpolarized quark into a transversely polarized hadron 
(denoted consistently as A N D c i / q ). This function, also known as the polarizing FF, could play a role in 
describing the observed transverse polarization of A hyperons produced in unpolarized pp collisions [137j . 



17 



For the sake of clarity we give here the relation between two standard notations for the most 
relevant TMD functions appearing in the sequel (the Sivers and Boer-Mulders distributions; the Collins 
and polarizing fragmentation functions): 



A N f q/p i(x, \k ± \) = -^ff^fr |fcx|) A^ T/p (x, \k ± \) = --j^hi\x, \k ± \) 
A N D h/q , (z, \k ±h \) = ?£±*l H^(z, \k ±h \) A N D h y q (z, \k ±h \) = %l D#(z, \k ±h \) . (16) 



It is also customary to consider /cj_-moments of the TMD distributions. The most relevant ones in the 
PDF sector are defined as (similar relations hold also for the FF case, see e.g. Eq. (1761) ) 



f$\x) = Jd 2 k ± f ± {Xi fHW( x) = J # k JkA f ± {Xi | fe± |) . (17) 

3.1.1 QCD developments 

An important and clear improvement of the TMD approach, in particular concerning the distribution 
sector, came with the discovery of a mechanism to generate a transverse spin asymmetry at leading 
twist [138] . Here, by considering the deeply inelastic lepton-proton scattering process, the authors 
showed that final-state interactions due to gluon exchange between the outgoing quark and the target 
spectator system lead to single spin asymmetries at leading twist in pQCD; i.e., the rescattering correc- 
tions are not power-law suppressed. The existence of such SSA's requires a phase difference between two 
amplitudes coupling the proton target with J* = ±| (where J p is the total proton angular momentum) 
to the same final state. These same amplitudes are necessary to produce a nonzero proton anomalous 
magnetic moment. The calculation of Ref. |138] was performed in a field-theoretic model in which QCD, 
with massive quarks, is supplemented by a coloured scalar diquark field and an elementary proton field. 
Soon afterwards this mechanism was reanalyzed by Collins |139j who proved that it is compatible with 
factorization and due to a spin asymmetry in the k^ distribution of quarks in a transversely polarized 
hadron: the "Sivers asymmetry". The earlier statement by Collins, that the Sivers asymmetry has 
to vanish because of time-reversal invariance, gets invalidated by the path-ordered exponential of the 
gluon field in the operator definition of TMD parton densities. What the time-reversal argument shows 
is instead that the Sivers distribution is reversed in sign in Drell-Yan processes w.r.t. the SIDIS case. 
This result, which would imply a violation of naive universality of parton densities, offers a clear test 
of our understanding of the TMD approach in pQCD. 

The origin of this sign change can be understood by looking at the gauge link entering the correlator. 
Therefore it is crucial to start with the properly defined gauge-invariant parton density for a hadron 
with momentum P and transverse spin S (f, j), that has the following operator definition 

f q/A ,(x,k ± ) = J d ^ ^fr- e -fap + r-Hfcx-ei. (P; s$(0, r , { ± )w^W 0oo m\P, S) ■ (18) 

Here, light-front coordinates are used: ^ = (£° ±£*)/\/2 and ^ = (0, £_l,0). The symbol W^oo 
indicates a Wilson-line operator, also called gauge-link, going out from the point £ to future infinity. 
These links are path-ordered exponentials, defined as [1401 1141] 



Wfoo = V exp ( - ig s dz A + (z • Vexp ( - ig s dz± ■ A ± (z =oo,z±f), (19) 



where = J2 a t a with t a = X a /2 being the Gell-Mann colour matrices. 

In particular, in the light-cone gauge (A + = 0) the first piece in Eq. (Tl9|) reduces to one and we are 
left with the pure transverse gauge link, depending on the nonvanishing vector potential at infinity. As 
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extensively discussed in Ref . |141] , this link in the transverse direction is crucial to maintain the gauge 
invariance of TMD parton distributions under residual gauge transformations; at the same time, it is 
responsible for the final(initial) state interactions in SIDIS(DY) \138\ 1139] . In nonsingular gauges, like 
covariant gauges, where the vector potential at infinity vanishes, in Eq. ( fl9j) only the gauge link along 
the light-like direction, z^ ~ (0,0, z~), survives. In this case, as discussed in Ref. |121j and reassessed 
in Ref. |142] severe light-cone divergences arise from contributions of virtual gluons with vanishing 
light-cone plus momentum. A standard way to cut them off is the use of a Wilson line slightly out of 
the light-like direction (z + ~ 0) |143j . Notice that this complication does not affect model calculations 
of TMD distributions at lowest order. 

The original proof on the vanishing of the Sivers function [20] , obtained applying space- and time- 
reversal to the quark fields in the operator definition of the parton densities, was then incomplete, 
because it ignored the presence of Wilson lines. Under time-reversal the future-pointing Wilson lines 
are replaced by past-pointing Wilson lines so that the correct version of the proof gives |139j 



Note the change in direction of the transverse spin- vector (f goes into |). Since the past-pointing 
Wilson lines are appropriate for factorization in the Drell-Yan process |139l 1144] . the correct result is 
not that the Sivers asymmetry, Eq. ([6]), vanishes, but that it has opposite signs in DIS and in Drell-Yan: 



The TMD spin-independent parton distribution keeps its universality. 

A complementary study for the fragmentation sector has been developed in Refs. [1451 1146] , where 
the authors compare SIDIS processes and e + e~ annihilation into two hadrons. The relevant result 
is that spin and k± dependent fragmentation functions are universal, also in sign. This finding has 
played a crucial role in the first extraction of the transversity distribution from a global analysis of the 
azimuthal asymmetries observed in e + e~ — > tttt + X and SIDIS (see sections 15.41 15.5.2p . 

The role of the gauge links in TMD distributions has then been further and systematically inves- 
tigated in a series of papers [1401 1141[ 11221 11231 I147[ 1148] , which focus mainly on the gauge invariant 
definition of TMD distributions and the proper factorization in QCD at leading twist for SIDIS and 
Drell-Yan processes at low transverse momentum. We only mention here that according to the k± 
factorization for SIDIS, in the convolution of kj_ dependent PDF's and FF's with the hard scattering 
parts another term appears: a soft factor coming from soft gluon radiations and defined by a matrix 
element of Wilson lines in QCD vacuum |122[ 1146] . One of its effects is that the transverse momentum 
of the observed hadron in SIDIS could be generated not only by the intrinsic motion of the quarks in 
the nucleon and in the quark fragmentation process but also via soft gluon radiation. Notice that, at 
the present stage, all phenomenological studies have neglected this extra unknown factor. 

In Ref. |147] the energy evolution of the spin and TMD distributions is also discussed. The authors, 
following Ref. |121j . show how the Collins-Soper equation in the impact parameter space can be applied 
also to these distributions. This equation, intimately related to the regularization of the light-cone 
singularities of TMD distributions in covariant gauges, allows to resum the large logarithms arising in 
perturbative calculations for SIDIS (or DY) process at low transverse momentum. 

Another interesting issue related to the TMD azimuthal asymmetries is the potential suppression 
due to Sudakov factors coming from soft gluon radiation. This problem has been analyzed in Ref. |149] 
for a class of processes, like the Drell-Yan process or two-hadron production in e + e~ annihilation, where 
two energy scales (for DY: the transverse momentum of the lepton pair, q?, and its invariant mass, Q) 
are present. In order to extend the TMD factorization picture beyond the region of very small qx (where 
tree-level formulas are sufficient) to the region where the transverse momentum becomes moderate (still 
Qt *C Q) one has to include the Sudakov factors by proper resummation. This implies a broadening of 
the qx distribution and a suppression effect, that becomes more important with rising energy |149] . 



fq/A^V^i ^J-) (future— pointing W — fq/At(.-Ey ^J_) I past— pointing W ■ 
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The scale dependence, at LO, of the first moments of the TMD distribution and fragmentation 
functions has been studied in Ref. [150] . by using Lorentz invariance and the QCD equations of motion 
in the large N c limit. A phenomeno logical study on the Q 2 -evolution (via resummation of soft and 
collinear parton emissions at LO) of unpolarized TMD distributions in SIDIS has been presented in 
Refs. [T5TirT52] . 

Before completing this short overview of the TMD approach we want to mention the recent study 
of the gauge links in more complicated processes, like pp — > h\h2 + X, where the elementary scattering 
involves many subprocesses and hadrons are present both in the initial and final state. 

Indeed, in the simplest hadronic scattering processes, such as SIDIS, DY and e + e~-annihilation, 
only a limited number of different gauge-link structures appear. These are the future and past-pointing 
Wilson lines mentioned above. When going beyond these processes one may encounter more complex 
gauge-link structures. In Refs. [1531 I154j the paths that can appear in general scattering processes 
have been analyzed (some technical aspects are presented in Ref. j!55j ). It was then shown that the 
presence of gauge links can lead to modified hard parts, referred to as gluonic pole cross sections, for 
SSA's in p^p^mr + X [1541 1156] . By taking suitable fcj_-moments of the SSA's (with proper weights) 
the net effect coming from these highly non-trivial structures is to modify the usual sum of Feynman 
diagrams representing the standard partonic cross sections into gauge invariant weighted sums, with 
SU(3) colour factors as weights (analogous to the ±1 factors discussed above for SIDIS and DY). 

Very recently the double inclusive hadron (or jet) production in hadron-hadron collisions with a 
small qr imbalance, has received a special attention, in particular concerning the validity of the k± 
factorization. In Ref. [157] by adopting a simplified Abelian model, Collins and Qiu argue that fac- 
torization is violated, even in the extended version of the TMD approach including gluonic pole cross 
sections. On the other hand, another approach |158[ I159j . based on a careful analysis of one-gluon 
radiation contributions, led apparently to the opposite result. Some light on this apparent discrepancy 
has been shed in Refs. [16011161] and further clarified in Ref. [162j . showing a mutual consistency among 
the different approaches. 

Summarizing, the main findings are: i) by expanding the gauge link in the correlators entering the 
colour gauge invariant approach [154] 1156] at first order in the coupling constant one recovers the results 
of Refs. [158] 1159] ; ii) by extending, and properly adapting, the calculations of Refs. [1581 1159] to two- 
gluon exchange, a non factorizable term in the TMD unpolarized cross section appears, in agreement 
with the results found in Ref. [157] 1160] . This extra piece can be taken into account by redefining the 
gauge link entering the TMD parton distribution [161] . A similar complication (i.e. a non factorizable 
hard part, or equivalently, a redefinition of the gauge link entering the TMD distribution) for the SSA 
(Sivers effect) is expected at the next order of perturbation theory. In Ref. [162] . by assuming factoriza- 
tion at the diagrammatic level, the authors obtain a manifest gauge invariant expression for the TMD 
cross sections in terms of gluonic pole cross sections and TMD distributions with pro cess- dependent 
Wilson lines. In this sense the TMD distributions are in general non universal. In this approach, beside 
confirming the results of Refs. [1581 1159] , gluonic pole cross sections in TMD unpolarized cross sections 
as well as ordinary partonic cross sections in unweighted SSA's arise together with universality-breaking 
parts, in agreement with Refs. [160] 1161] . However, these pieces are shown to vanish for the integrated 
and weighted observables considered in Refs. [154] 1156] 1163] 1164] . The identification of these process de- 
pendent pieces could be an important step towards a unified picture of TMD factorization for hadronic 
processes. Further work is still needed. 

3.1.2 The impact parameter picture 

This physical picture, mostly applied to the Sivers effect, has been discussed in a series of papers by 
Burkardt and collaborators [3~Tl [32] [33l 1165] . It is based on two main ingredients: the distortion in the 
transverse plane of the distribution of partons inside a polarized target and the attractive (repulsive) 
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final (initial) state interactions in SIDIS(DY). The first aspect is strongly related to the information one 
can extract from the generalized parton distributions (GPD). 

GPD's are generalizations of ordinary parton distributions to non-forward matrix elements of a 
lightlike correlation function. They can be accessed for instance in deeply virtual Compton scattering, 
7*p — > 7p', where a finite momentum is transferred to the proton. GPD's depend on two longitudinal 
momentum fractions (x and x + £, where £ is the so-called skewness parameter) and on the invariant 
momentum transfer to the proton, t = A 2 . In the simpler case where £ = the Fourier transform of the 
GPD w.r.t. t gives the parton distribution in position space, q(x, b), where b is the impact parameter 
giving the transverse distance of the struck quark from the center of mass of the proton. 

When the proton is transversely polarized (for instance along the +Y axis) two GPD's describe the 
distribution of unpolarized quarks, namely H(x,0,t) (related to the unpolarized PDF) and E(x, 0,t), 
whose Fourier transform gives the distortion of the distribution. 

By relating the integral of E q (x, 0,0) to the g-flavour contribution (n q ) to the proton anomalous 
magnetic moment one can get an estimate of such distortions. On this basis one can expect an opposite 
sign of the distortion for u and d quarks, with u quarks mostly displaced along —X. Notice that this 
left-right asymmetry in fe-space is T-even, in contrast with the left-right asymmetry in momentum space. 

Let us come to the second ingredient in this picture, by considering SSA's for the production of 
7r + , 7T° in SIDIS. Due to charge factors and the fact that u — > ir + , 7r° fragmentation is favoured, most of 
7r + , 7T° mesons come from an initial u quark in the polarized proton. We now recall that any left-right 
asymmetry in the transverse momentum space at the quark level requires final state interactions (FSI). 
Intuitively one then expects that the FSI's are on average attractive, since it costs energy to build up the 
string of gauge fields that connects the struck quark with the spectators. More precisely, if the photon, 
moving along the — Z axis, collides with a nucleon polarized along the +Y axis, from the previous 
results we expect that a u quark tends to be displaced in the —X direction in impact parameter space. 
If the FSI's are attractive, this quark experiences a force along the +X direction. This translates into 
the preferred direction of 7r + , it mesons along +X, in agreement with model calculations of the Sivers 
function as well its phenomenological extractions. In the Drell-Yan process the initial state interactions, 
expected to be repulsive, imply an opposite sign for the Sivers function. 

3.1.3 Models and constraints for TMD distributions 

The Sivers function, being a naively T-odd entity, requires the interference between two complex am- 
plitudes with different phases. Spectator models at tree level cannot provide these nontrivial phases. 
However they can arise as soon as a gluon is exchanged between the struck quark and the target spec- 
tator |138j . In other words the presence of the gauge link, which insures the colour gauge invariance of 
parton distributions, at the same time could provide the nontrivial, relevant, phases. 

Starting from the work of Brodsky, Hwang and Schmidt (BHS) [138| . different models for the Sivers 
function as well as for the Boer-Mulders function have been developed. In Ref. [166] the authors 
extended the BHS calculation, based on a scalar spectator diquark model of the nucleon with a point- 
like form factor, to DY processes. They also computed the Boer-Mulders function for u quarks, finding 
h>i = /xr ■ The same relation was found in Refs. [167[ I168j . still with a scalar spectator but adopting 
a Gaussian form factor. In both cases a negative f^p for u quarks came out. The MIT bag model 
was adopted by Yuan in Ref. |169j . where he still found a negative Sivers function for u quarks and 
f£? = -(1/4) f$ (model [A] in Fig. GQ left panel). Moreover he obtained h\ u = (3/2) and 
= — = (3/4) f^ L . In Ref. |170j the Sivers function has been calculated in a spectator model 
with scalar and axial- vector diquarks (model [B] in Fig. [H left panel), with yet a different choice for the 
form factor w.r.t. Ref. [168] . The d quark Sivers function turns out to have the opposite sign compared 
to the u quark one and to be much smaller in size; again ~ fyp. 

Analogous results have been obtained in Ref. |171j both with scalar and vector diquarks in a light- 
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cone SU (6) quark- diquark model. The instanton liquid model for QCD vacuum with MIT bag model 
wave functions for quarks was adopted in Ref. [172] (model [C] in Fig. [TJ left panel). In this model, the 
two independent perturbative and instanton terms sum up to give the total contribution to the Sivers 
function. For d quarks the two terms almost cancel each other, leading to a small, negative Sivers 
function. On the contrary, in the case of the w-quark Sivers function, the instanton contribution and 
the perturbative one add together in fair agreement with the other models. 

From these models the emerging common result is the sign of the Sivers function for u quarks. It is 
worth to stress that while this is in agreement with the phenomenological extractions, a clear discrepancy 
is found concerning the overall size. In this respect, the study in Ref. |173j . based on the QCD counting 
rules, could help in putting useful constraints on the Sivers function, which should be one power of 
(1 — x) suppressed relative to the unpolarized PDF. Concerning the Boer-Mulders function new results 
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Figure 1: Left panel: Comparison of models calculations of x/ 1T for up (thick lines) and 
down (thin lines) quark Sivers function according to Refs. [169] [A], [170] [B], and [172] [C]. 
Right panel: Model calculations of H^ 1 ^ 2 ^ for the favoured Collins function normalized to 
the unpolarized FF [174] . Results for pseudoscalar (solid line) and pseudovector (dashed 
line) pion-quark coupling with gluon loops are shown. 

have appeared [175[ 1176] . With the exception of a previous calculation [170] . all models, including the 
new findings in the diquark spectator model [175] predict a negative Boer-Mulders function for both 
u and d quarks. The same conclusion has been reached in Ref. [176] through the calculation of the 
chirally odd GPD Et (which has opposite sign w.r.t. [33]) in different models both for nucleon and 
pion targets. 

Let us now consider the Collins function. Four model calculations for the fragmentation of a trans- 
versely quark into a pion have been presented so far in the literature [1771 11781 1179[ 1180] (besides the 
former attempts by Collins [2D] and Artru [128] , discussed above). All of them produce the necessary 
imaginary parts by adopting a simple model for the fragmentation process at tree level and inserting 
one-loop corrections. Two possibilities for the tree-level amplitudes and two possible kinds of one-loop 
corrections have been investigated, for a total of four different models: pseudoscalar pion-quark cou- 
pling with pion loops [177] and with gluon loops [17911174] : pseudovector pion-quark coupling with pion 
loops [178] and with gluon loops [180] . For gluon loop cases see Fig. [T] (right panel). We refer the reader 
to Ref. [174] for further details. 

In the above-mentioned approaches disfavoured fragmentation functions (e.g. d — > tt + ) vanish, some- 
thing apparently not supported by data, see section [5j In principle, these functions could be calculated 
in these models by considering diagrams with the emission of two pions, one of which goes unobserved. 
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The different results obtained for the Collins function do not allow to draw any conclusion even 
about its sign. The fact that more than one diagram with a different structure contributes to the 
Collins function makes it also difficult, if not impossible, to interpret the resulting effect, as we do for 
the Sivers function, in terms of simple attractive/repulsive interactions. 

Some useful constraints on TMD distribution and fragmentation functions in terms of bounds and 
sum rules are available. In Ref. |181j by considering the combined Dirac ® target spin space and requir- 
ing the positivity of the spin-matrix eigenvalues, they found, beside the well-known Soffer bound |182j 
(2 1 < (/i + 9il))j new bounds on the eight TMD distributions entering the correlator in Eq. (115j) . 

The first sum rule for a T-odd TMD function was discussed in Ref. |183j : the authors considered 
the Collins function and, by imposing conservation of the intrinsic transverse momentum in parton 
fragmentation, obtained 

£ / dzzHt {l)q {z) = Y,l dz f d 2 k ± ^A N D h/q ,(z,\k ± \) = 0, (22) 

h J h J J M h 

where H^~^ 9 (z) is defined analogously as in Eq. (j!7p and the sum is over all hadrons. 

More recently a somewhat corresponding sum rule for the Sivers function has been derived by 
Burkardt [184] . This rule states that the net transverse momentum due to the Sivers mechanism must 
vanish; more formally 

(k±)= J2_ (Ma= fdxfd 2 k ± k ± J2_ A/ o/pT (a:,fe±) = 0. (23) 

a=q,q,g a=q,q,g 

Even if intuitively expected, the nontrivial fact is its validity in presence of final state interactions, which 
might spoil the simple partonic interpretation. The Burkardt sum rule (BSR) has been also checked 
explicitely in the diquark model calculations of the Sivers function |185j . Its role in SSA phenomenology 
will be discussed in Section 

3.2 Twist-three effects in collinear pQCD 

The origin of this approach could be traced back to the early paper by Kane et al. [IS], where it was 
shown how in a collinear QCD parton model any single spin asymmetry should be proportional to the 
quark mass. Although this result cannot explain the large values of An observed, it was read as a 
signal of a higher twist effect. Another step towards this interpretation was developed by Efremov and 
Teryaev [186[ I187[ 1188] who pointed out that a nonvanishing SSA can be obtained in pQCD beyond the 
leading power expansion. However, the first consistent calculation in collinear pQCD of a sizable SSA 
in the forward region (large xp) was given by Qiu and Sterman in their paper on hadronic direct photon 
production [1891 119U] and their subsequent application to pion production [191] . The asymmetries are 
evaluated in the framework of pQCD generalized factorization theorems with the introduction of new 
twist-three quark-gluon correlator functions convoluted with ordinary twist-two parton distribution 
functions and a short-distance hard scattering part. The new functions, being expectation values 
between the hadronic states of three field operators, do not have a simple partonic interpretation. More 
precisely, for the numerator of in the process A^ B — > C + X (C unpolarized) they obtain, neglecting 
higher power corrections (we keep here the original notation of Ref. [191] ). 

dAa A i B _> c+x = ^I/at ( z i> z 2) ® <f>b/B(x') ® H a+b ^ c g> D c ^ c (z) 

abc 

+ E *<8)a{*) ® 4) ® K + b^ c ® A~c(s) (24) 

abc 
abc 
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where a, b, c represent parton flavours (q, q, g), and 06/_b(V) and D c ^c{z) are ordinary twist-two parton 

(2) 

distribution and fragmentation functions, respectively; 5q a j A {x) is the twist-two transversity distribution 
function. The new pieces, namely twist-three functions (see apexes), depend on two light-cone variables 
and can give three different contributions, from the initial polarized hadron, the initial unpolarized 
hadron, the final unpolarized hadron. As in ordinary pQCD factorization, only the hard scattering 
functions H are calculable and their forms depend on the definition of the twist-three distributions. 

In order to present this formalism we consider in some detail the first term (first line) in Eq. (124jl . 

Keeping in mind that the asymmetries are larger in the forward region, one can restrict to the 
valence flavour approximation for a in J2abci an d then consider only a valence twist-three distribution. 
Another important simplification comes from the fact that the H a+ b_> c part involves two contributions: 
one proportional to 5(x\ — x 2 ) and another one proportional to 5{xj). The first case sets the gluon 
momentum in the correlator to zero; in the second case, one of the quarks carries a vanishing momentum. 
These are referred to as soft gluon and soft fermion poles respectively. Moreover, while soft gluon poles 
enter with derivatives of the twist-three distribution, soft fermion ones do not. If, as for ordinary twist- 
two functions, the derivatives of twist-three functions are enhanced near the edges of phase space, then 
soft gluon poles could give the leading contribution to the large asymmetries for xp —* 1. 

To get the master formula in this approach one has to start with the factorized expression in a 
diagrammatic representation of the process 

dAa^j (2^4 t% Ta{K k2)Sa{ku **) ' (25) 

where ki s2 are valence quark momenta, T a is proportional to the non perturbative matrix element of 
the quark-quark-gluon operator between polarized initial hadron states (a twist-three matrix element) 
and S a refers to the rest of the process. Expanding S a in the collinear approximation enables to reduce 
the four-dimensional integrals to convolutions in the momentum fractions of partons with fcj ~ a^P. 
Notice that in this approach a nonzero spin dependence is found only from pole terms in the hard 
scattering. Without them the symmetries of the strong interactions would forbid any SSA. Keeping the 
first nonvanishing term in dAa, that is the leading pole structure of S a , one gets 

r - (99 

dAa « Y, dx 1 dx 2 ief' STnfl —^(x 1 ,x 2 )T ajF (x 1 ,x 2 ) , (26) 
where n, n are light-like unit vectors and 

T a Axi,x 2 ) = |^L^e^ p+ ^ + ^-^) p+ ^(P, ST |^(0) 7 + [e™ F a + {y^)] 1> a fa)\P, s T ) , (27) 

with F a + the gluon field strength. Since T at p is real only the poles of S a can contribute. 

The last step is to factorize the remaining function S a into a perturbatively calculable partonic term, 
H a+ b^ c , a corresponding target parton distribution, <fib/B, and a fragmentation function, D c ^c- 



/dx' f 
— 4>b/B{x') dz H a+h ^ c {k u k 2 ,x' \p c ) D c ^ c (z) , 

DC X J 

where J2bc runs over an parton flavours. 

Finally, the full factorized expression for dAa reads 

[ r r r - 9 

dAa = — ^ JdzD c ^ c (z) J — <f) b /B{x') Jdx 1 dx 2 T atF (x 1 ,x 2 ) i e pSTnn — pH a+b ^ c (x u x 2 , x', z) 

^ abc - 2 



f28l 



k p 2 =0 



(29) 
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where the integration over either X\ or 22 can be done by using the poles in H a+b ^ c . 

Notice that this results in a factorization formula with only a single momentum fraction for each 
hadron, similar to that for the spin-averaged cross section with <j> a /A(x) replaced by T a;F (x, x). 

A way to understand why T ajF is expected to be the dominant contribution is the following. From 
the sources of the ki (i = 1, 2) dependence in H a+b ^ c (ki, x',p c ) (on-shellness of the unobserved final 
parton d, on- and off-shellness of the propagators and the explicit dependence in the numerator) we get 
two kinds of terms: one proportional to (d/dx)T a)F (x, x) and one proportional to T atF (x, x). 

One has to remember that this approach was intended to address the description of the asymmetries 
in the forward region, where Xf is large. Here the dominant contribution comes from large net momen- 
tum fraction x for the polarized beam hadron, coupled with relatively small momentum fraction x' from 
the partons of the unpolarized target hadron. Since all distributions vanish for large x as (1 — x)@, with 
(3 > 0, (d/dx)T a>F (x,x) ^> T a)F (x,x) when x — > 1. Therefore, in the forward region, terms proportional 
to derivatives of the distributions T aF dominate. Moreover it is only the matrix element T UjF that 
inherits derivative terms, as a result of the collinear expansion involving soft gluon poles. Soft fermion 
poles have no such derivatives at LO and also do not correspond to the valence quark approximation. 
In this spirit only the derivative term was kept and the corresponding hard part was calculated. 

A recent development of this approach has been discussed in Ref. |192j . Still neglecting the soft 
fermion poles, the complete soft-gluon pole structure, including the non derivative term, has been 
computed. The main result of this work is that one can incorporate the two terms in a simple way since 
they have a common hard scattering part. This comes from a remarkable property of the hard scattering 
functions recently proven by Koike and Tanaka [193j . On the other hand the phenomenological role of 
this extra piece in the forward region is found to be negligible, as assumed in former papers. 

Without giving further details of the calculation (see [191} 1192] ). we only mention that by properly 
exploiting the pole structure of H a+b ^ c (ki, k 2 ) (initial and final-state poles), one gets for the numerator 
of the SSA, restoring all factors, 

d 3 Aa a 2 s r 1 dz ^ r 1 dx' 

d 3 Pc s ^ c Jz min z 2 J x ' min x' x's + t/ 



c- 



^r"Iz mi „ z 2 c ^°( ) j x i x > x's + t/z^ b ^ B ^ ^ ^ ^ 



X 



net; 




T a>F (x,x) -x (-^T a>F (x,x) 



H a+b _> c (s,t,u) 



where s, t, u (s, t, it) are the standard hadronic (partonic) Mandelstam variables, and x, x' m 
fixed by the usual elastic scattering condition. 

The H a+b ^ c are the final hard-scattering functions and read 



H a+b ^ c = H'(s, t, u) + HF(s, t, u) I 1 + j I , (31) 




where H^ +b ^ c (H^^) denotes the contributions due to initial (final) -state interactions. An important 
feature is that they are given in terms of the same Feynman diagrams needed to calculate the usual 
unpolarized partonic cross sections for the ab — > cd process but with different colour factors due to 
the extra initial (final) state interactions, see Appendix of Ref. |192j . Two aspects of the final result 
in Eq. fl30|) deserve a word: the denominator u gives explicitly the power suppression of the SSA and 
the extra factor \/a s w.r.t. the unpolarized cross section comes from the additional interaction with a 
gluon field. 

It is also worth to mention that both in the original and the more recent works the contributions 
involving the three-gluon twist-three correlation function, which would enter the gg partonic process 
and could probably be relevant in the mid-rapidity region, have been ignored. 

In order to perform a phenomenological study one has to rely on some parameterizations of the 
T a ,F twist-three functions. Since the operator defining it contains the same quark field operators as the 
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ordinary twist-two parton distribution functions, the most natural ansatz is: 

T a ,F(x, x) = N a (x)(j) a (x) . (32) 

In Refs. |189[ I19U[ I191j a simple functional form with N a (x) = n a \ was adopted, where A has the 
dimension of a mass and k u ^ = ±1 (for a proton). In the recent developments of this formalism 
|192j . with more data having become available, a more general form has been adopted, i.e. N a (x) = 
N a x aa (l — Xf a . In both cases the main features of the data for An in pion production in pp as well as 
in pp collisions can be described rather well. Notice, however, that all these papers use leading order, 
leading twist expressions for the unpolarized cross section in the denominator of An- 

This approach, applied to SSA's in prompt photon production, by adopting the same twist-three 
functions as extracted from the analysis of pion data (positive for up and negative for down quarks) 
gives a sizable negative asymmetry. This result, as we will discuss in section [51 could be of extreme 
relevance in discriminating between the GPM approach, based on TMD distribution and fragmentation 
functions, and the twist-three formalism. 

A complementary study has been also devoted to the role of the other two sources of SSA's in 
pp — > ir + X in the twist-three approximation, namely the chiral-odd contributions in Eq. (1241) (second 
and third lines). The twist-three distribution, (p^L^x'^x'^} and its relative hard scattering parts have 
been discussed in Refs. |194[ 1195] . This contribution, expected to be relevant in the negative xf region, 
comes out to be negligible over all ranges in xp due to the smallness of the hard scattering partonic 
asymmetry. In Ref. [196] a detailed analysis of the chiral-odd twist-three fragmentation function was 
performed, showing how, with some difficulties, a description of E704 data can be obtained. 

Remaining in the context of inclusive particle production in hadron-hadron collisions, this approach 
has been also applied to the study of the transverse A polarization in unpolarized pp collisions. In 
this case two twist-three effects could contribute: the twist-three chiral-odd distribution <f>a) p (xi,%2)i 
coupled with a twist-two unpolarized PDF and a twist-two chiral-odd fragmentation function (analogous 
to the twist-two transversity PDF); a twist-three chiral-even FF, that, following the same notation 
adopted in Eq. fl24l) . would read D c ^^ (z\ , z?) , coupled with two twist-two unpolarized PDF's. The 
first contribution has been studied in Ref. [197] . where by adopting a simple ansatz on </>^ p (xi, x 2 ), 
namely a proportionality relation with the chiral-odd transversity distribution, the main features of the 
data are reproduced, although some discrepancies still remain. 

The twist-three formalism has been also applied to SIDIS azimuthal asymmetries [198, 199J at large 
p T ~ Q (with Q the virtuality of the exchanged photon). In these papers the authors show that among 
the two possible mechanisms, the first coming from the twist-three chiral-even distribution and the 
second from the twist-three chiral-odd fragmentation (coupled with the transversity PDF), the last 
one is negligible. In particular they calculated the soft-gluon pole contributions for both cases. They 
also improved their calculation for the first and most relevant case, identifying all pole contributions 
(soft-fermion and hard-gluon poles included) and providing a systematic collinear expansion for this 
approach in SIDIS. 

Before concluding this section we mention some attempts to connect the TMD approach to the twist- 
three formalism. A first study has been presented in Ref. [2UU] . showing that the inclusion of transverse 
gauge links, which lead to T-odd effects, is connected to the Qiu-Sterman mechanism. More precisely, 
it was shown that the first k± moment of the Sivers function is equal to the twist-three chiral-even 
distribution, T a ^. This connection has been also addressed in Ref. [2 1 J . 

Recently a unified picture in terms of the twist-three approach and the TMD formalism |2U2} \'203\ 
1204^ I205j both in SIDIS and Drell-Yan processes has been presented. The main result is that although 
the two mechanisms have their own domain of validity, they describe the same physics in the kinematic 
region where they overlap. More precisely, focusing on SIDIS for the sake of clarity, when the transverse 
momentum of the observed hadron, pt, and the photon virtuality, Q, are much larger than Aqcd, the 
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spin-dependent cross section can be calculated in terms of a twist-three quark-gluon correlation. On the 
other hand, when px ~ Aqcd *C Q, single spin asymmetries can be generated from a spin- dependent 
TMD quark distribution. There is, however, a common kinematic region, Aqcd <S Pt <C Q, where both 
mechanisms should work. In this region px is large, so that the asymmetry is a twist-three effect, but 
at the same time px <C Q, so that the TMD factorization formalism also applies. The same arguments 
work for the Drell-Yan process by simply replacing the px of the observed hadron with the transverse 
momentum of the lepton pair. The resulting connection unifies the physical pictures for the underlying 
dynamics of transverse SSA's and imposes important constraints on phenomenological studies. 

We finally stress that both in the TMD and in the twist-three approach, on-shell massless partons 
are considered in the soft as well as in the hard parts. The role of off-shellness, still largely unknown, 
has been recently emphasized and investigated in Refs. [206[ 1207] I2U8] (see also references therein). 

3.3 A semiclassical model 

Among the nonperturbative models proposed to describe the SSA's observed in inclusive hadron pro- 
duction in hadron-hadron collisions we mention a picture based on a semiclassical geometrical approach: 
the orbiting valence quark model [9j 12091 12101 1211[ 12121 12131 1214] . The main ideas behind this approach 
can be summarized as follows: i) The constituents of a polarized hadron are assumed to perform an 
orbital motion about the polarization axis. Their relativistic ground state wave functions allow to cal- 
culate the polarization of valence quarks in a polarized proton. The results show that there is a flavour 
dependent asymmetry in the valence-quark polarization. The sea quarks come out unpolarized. ii) 
In the fragmentation region the main contribution to hadron production comes from quark- ant iquark 
annihilation via a direct-formation mechanism. Hi) Since hadrons are extended objects, a surface effect 
is expected in such production processes. 

Let us consider a proton-proton collision along the Z axis in the rest frame of the polarized proton. In 
the transversely upward {+X) polarized hadron, u quarks, being on average polarized along the hadron 
polarization, have an orbital motion such that they go preferably left (-Y) at the front surface and 
right (+Y) at the back surface of the hadron. When the polarized target interacts with the unpolarized 
projectile this preferred direction is preserved at the front surface (when the projectile and the target 
start overlapping) while it is partially spoiled at the back surface (when they start separating again) 
because of the colour interactions acting during the transition of the projectile through the target. 
The net result is that in the fusion mechanism of a u quark from the polarized proton and a suitable 
antiquark from the unpolarized proton, the formed meson goes preferably left, i.e. along —Y. More 
generally this means that, when a polarized proton collides, only colour singlet qq pairs formed at the 
front surface can get an extra momentum due to the orbiting motion of valence quarks. This determines 
the observed asymmetry in the distribution of the produced mesons. 

In this model one then expects that for p^p — > tt + X in the fragmentation region of the polarized 
proton, where mainly valence quarks contribute, 7r + and 7r° mesons go left, whereas 7r~'s go right. For 
the case of a polarized antiproton 7r + 's and 7r~'s reverse their direction, whereas 7r°'s keep going left. 
The left-right asymmetry is more significant at large Xp. According to this mechanism, also Drell-Yan 
lepton pairs should exhibit a left-right asymmetry. 

One of the main features of this model is that in order to have a nonzero SSA the produced hadron 
must have at least a common valence quark with the polarized hadron. In particular, for kaon production 
in p^p collisions one expects that for K + (us) An is similar (in size and sign) to An(tt + ), whereas for 
K~{us) it should vanish. 

This model has been also applied to the description of the transverse polarization of hyperons [Py(Y)] 
produced in unpolarized hadron-hadron collisions [215[ 1216] . 

For a detailed and comprehensive description of this approach and its phenomenological applications 
see Ref. [9] and references therein. 
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4 Spin asymmetries in the generalized parton model 



In this section we will present and discuss in detail the generalized, spin and k± dependent, parton model 
approach to polarization effects in high-energy hadronic collisions. As discussed in section [3j there are 
several approaches available in the literature. In what follows we will adopt the one developed in a series 
of papers, starting from the middle of 90s. This approach is an extension of early attempts to deal with 
parton intrinsic motion and unintegrated, transverse momentum dependent (TMD) parton distribution 
(PDF) and fragmentation functions (FF). It consistently includes spin and polarization effects, requiring 
the introduction of a new class of leading twist TMD distributions. Like any theoretical approach, the 
formalism we will use has both advantages and drawbacks. The main advantages are the following: a) 
It is a QCD-improved, generalized parton model approach, with direct inclusion, in the spirit of the 
parton model itself, of spin and transverse momentum (fej_) effects; b) It retains the simple partonic 
interpretation of all soft, leading-twist TMD distributions; c) Using the helicity formalism, it very clearly 
shows the connection among the polarization states of the particles involved and the role of spin effects 
in the soft and hard processes; d) It fully includes kinematical k± effects, without approximations. It 
has of course also some drawbacks: a) For some of the processes considered, factorization has not been 
proven explicitly yet. It is rather taken as a reasonable assumption for this generalized scheme; b) It 
only accounts for leading-twist soft distributions (in this sense its application to subleading-twist SSA's 
can be considered as an effective approach); c) It does not properly account for initial and final state 
interactions, which have been shown to be crucial for the nonvanishing of naively time-odd PDF's and 
FF's. Despite these problems, its phenomenological success is presently remarkable. In our opinion it 
is also the most convenient approach for a simple illustration of the pQCD-based formalisms and of 
their main phenomenological applications. In section we will discuss and compare the results and 
predictions of this approach with available experimental results and with those of other theoretical 
approaches. In particular, we will emphasize measurements which could help in disentangling among 
different formalisms and are at reach of present experiments. 



4.1 Polarized cross sections and helicity formalism 

Let us consider, to be definite, the doubly polarized hadronic process A(Sa) + B(Sb) — ► C + X at high 
energy and moderately large transverse momentum, px- Here px is the transverse momentum of the 
observed hadron with respect to the direction of the colliding beams. A, B, C are generic hadrons, 
respectively with spin state Sa, Sb, Sq (5/, = 0, 1/2, 1). Using the helicity formalism the invariant 
differential cross section for the production of an unpolarized particle in the polarized processes under 
consideration can be expressed through the following master formula: 

° d „ Vr = £ / ZtXX rfWV t rfW(*xc ■ Pc) J(h w ) (33) 

" PC a,b,c,d,{X} J iU7r X a X bZ b 

x Px a A xi A fa/A,s A (x a , k ±a ) p b ^f B f b/B ,s B (x b , k± b ) Mx c ,x d ,x a ,x b M^.y^ 5(s + t + u) D X x ^ x \ c (z, k ±c ) . 

Notice that with generic X c , X' c helicities (i.e. replacing J2x c D x c y C with D x c y c ) on the r.h.s. of the 
above equation, we get the (unnormahzed) helicity density matrix for the final particle C that allows 
to evaluate its polarization state. This is particularly relevant for the case of transverse hyperon 
polarization in unpolarized hadronic collisions. However, since here we want to focus on the basic ideas 
of the approach, for simplicity we will consider the production of unpolarized particles in the collisions 
of spin 1/2 hadrons. A detailed treatment of spin 1/2 particle production using the helicity formalism 
in a generalized parton model is under completion |217j . We will shortly discuss the formalism for the 
spin 1/2 case in section 14.3.21 and the phenomenology for transverse A polarization in section [5.1.31 
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Eq. (1331) gives the cross section for the polarized hadronic process under consideration as a factorized 
convolution of all possible hard elementary QCD processes, ab — > cd, with soft, spin and k± dependent 
PDF's and FF's. Let us briefly comment on the notation and physical content of this equation. Further 
details can be found in the original papers |135l 12181 1219j and in references therein. 

(i) A and B are initial, spin 1/2 hadrons in pure spin states denoted by Sa and Sb and corresponding 
polarization (pseudo) vectors P A and P B , respectively; 

(ii) Unless differently stated, we always consider the cm. frame of the initial colliding hadrons, the AB 
cm. frame or hadronic cm. frame; in this frame, hadron A, B, move respectively along the ±Z cm -axis 
direction and the final hadron C is produced in the (XZ) cm plane, with (pc)x cm > 0; 

(iii) The notation {A} indicates a sum over all helicity indices; x a , xj, and z are the usual light-cone 
momentum fractions of partons a, b inside the initial hadrons A, B and of hadron C produced in the 
fragmentation of parton c; k^ a (k±b) are the (two-dimensional) transverse momenta of parton a(b) with 
respect to hadron A(B); k±c is the generic three-momentum of hadron C in the hadronic cm. frame; 
notice that the delta function 5(k±c ' Pc) constrains k±c to be effectively two-dimensional and always 
orthogonal to parton c three-momentum, p c ; 

(iv) The factor (167r 2 x a Xfes) _1 , where s = (j>a + Pb) 2 is the total energy of the colliding beams, collects 
phase space factors related to the elementary cross sections and the corresponding flux factor for the 
noncollinear (when viewed in the hadronic cm. frame) parton scattering process; 

(v) The phase-space factor J(k±c) / z 2 is the invariant TMD Jacobian factor for massless particles con- 
necting the parton momentum p c with the hadron momentum pc, (see e.g. Appendix A of Ref. [218] ); 

(vi) p\ A y A is the helicity density matrix of parton a inside hadron A with spin state Sa, and analogously 
for p b / B jf B ; notice that the helicity density matrix always describes the polarization state of the particle 
in its helicity rest frame; for a massless particle no rest frame exists and the helicity frame is defined 
as the standard frame, reached from the hadronic cm. frame, in which the particle four-momentum is 
p" = (p,0,0,p) PI; 

(vii) f a /A,s A { x a, k± a ) is the spin and TMD distribution function of the unpolarized parton a inside the 
polarized hadron A, and similarly for parton b inside hadron B; these leading-twist TMD distributions 
generalize the usual partonic distributions in the collinear configuration; notice that in this paper, for 
the sake of clarity, we adopt the following convention: for all k^ dependent distribution and fragmenta- 
tion functions and soft amplitudes we will use a "hat" to indicate full, vectorial dependence on k±, while 
we will not use it for dependence on k± — \k±\, with factorized azimuthal dependence; notice also that 
this notation is partially at variance with those adopted in the papers where the original formulation 
of the GPM was developed, which by the way slightly differ from paper to paper; 

(viii) The Ma c ,a ;A a ,A ' s are the helicity amplitudes for the elementary process ab — > cd, opportunely 
normalized according to the kinematical phase-space factors explicitly shown in Eq. 

"A A' 

(ix) D x c \, c (z, k±c) is the product of nonperturbative fragmentation amplitudes for the polarized frag- 
mentation process c — > C + X . 

We will now consider in more detail the quantities which describe how the polarization of the initial 
hadrons is transferred, through the soft and hard elementary processes, to the final partons and to the 
produced hadron (if its polarization state is observed). 



4.2 Spin and TMD parton distribution functions at leading twist 

The quantity Px A y A fa/A,s A ( x a, k± a ) in Eq. (1551) encodes complete information on the polarization state 
of parton a with spin s a , as determined by the polarization state of the parent hadron A (which is fixed 
by the experimental conditions considered) and by the soft, nonperturbative processes related to the 
hadron structure (analogously for parton b inside hadron B). Indeed, introducing soft, nonperturbative 
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helicity amplitudes for the inclusive process A — > a + X, F\ a ,\ XA ;\ A (x a , k± a ), we can write: 



where X stands for a spin sum and phase space integration over the undetected system of hadron 

A remnants, p x ' x , is in turn the helicity density matrix of hadron A, 



a,s a If 1 + Pi P£-^P4\ 1 ( 1 + P£ P$ e-^-* \ 

P ^' A ~ 2\P£ + iP4 I -Pi ) ASA ~2\P£e^ \-P£ ) ASa > 

P A = (P<p cos <ps A , P-r srn 0.5U' Pl) * s its polarization vector and 4>s A the corresponding azimuthal 
angle, defined in the helicity reference frame of hadron A. We have also defined 

F\ A ,\ A ( X "> k±a "> = f< XA A £k,*x a *a( x "> k±a "> ^x' a ,x XA -,x' A ( x a, k± a ) ■ (36) 

"A A' 

The product of nonperturbative helicity amplitudes F x a ' x f has a simple physical interpretation and is 
directly related to the leading twist part of the well known hand-bag diagrams for the hadronic correlator 
in deeply inelastic scattering processes. The distribution amplitudes ^ r x a ,x XA -x A depend on the parton 
light-cone momentum fraction x a and on its intrinsic transverse momentum k± a , with modulus k± a and 
azimuthal angle <fi a (all defined in the parent hadron helicity frame), in the following way: 

F\ a ,\ XA ,\ A (x a , k ±a ) = J r x a ,x XA ,x A (x a , k ±a ) ex.p[i\ A (j) a ] . (37) 

Rotational and parity invariance imply that only eight functions as defined in Eq. (156"j) are independent. 
Notice that since some parity relations are different according to the parton type involved (quark or 
gluon), the two cases must be treated separately. For quarks, keeping in mind that upper (lower) 
helicity indices refer to the parton (hadron) respectively, one can easily associate the eight functions 

F++i F~, F+T, F;+, F-Z, F+Z , (38) 

to the eight independent, leading twist, polarized TMD quark distribution functions: 

1) The combinations F++±F+~Z are respectively related to the unpolarized and longitudinally polarized 
PDF's, f a / A and A L f a/A (denoted as ft and g lL in Ref. PMJ); 

2) The combinations F+Z ± F+Z are in turn related to the two possible TMD contributions to the 
transversity PDF, A T f a / A (corresponding to h 1T and h^ T in Ref. |108j ): 

3) The combinations Ftz ± F+Z are respectively related to the probability of finding an unpolarized 
(longitudinally polarized) quark inside a transversely polarized hadron, that is to the well-known Sivers 
function [TH [19] A N f a / A ^(or f^ [12]), and to the function Af a)Sz / At s T , (related to gir of Ref. [108] ): 

4) Finally, the combinations F+~+ ± FZZ are related to the probability of finding a transversely polarized 
quark respectively inside an unpolarized hadron (that is, to the Boer- Mulders function [108] . A N f a y A 
or hj~) or inside a longitudinally polarized hadron (that is, to the PDF Af a ^ T / A) s L or h^ L |108j ). 

One can also see that: the functions F++, are obviously purely real; the F+Z are purely real (purely 
imaginary) for quarks (gluons); the last four functions in Eq. (1551) are complex but not independent, 
and their sums and differences can be also expressed as the real and imaginary parts of F+Z and F+Z ■ 

Some additional comments are in order here to clarify the physical content of these functions and 
summarize their main features: among the eight functions, the only ones that survive in the collinear 
configuration (that is, under k± integration) are the three usual leading twist PDF's, f a /A, Ai,fa,/A, 
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At fa/ a- Those which are off-diagonal in the parton helicity indices, like Ay/ a /^, are chiral odd; there- 
fore, they cannot be directly measured in DIS processes with on-shell, massless partons. In any physical 
observable, they must combine with some other helicity-odd term (e.g., in SIDIS, with a chiral-odd FF, 
like the Collins function, or with higher twist terms due to parton masses and off-shellness); those func- 
tions proportional to the imaginary part of the F's are in turn naively T-odd, and require initial/final 
state interactions in order to be nonvanishing [1381 1139] (see also section [3]). 

The most general expression for the helicity density matrix of quark a inside a hadron, A, with 
polarization state Sa is 



a/A,s A = ( P a ++ P a +-\ U 1 + Pt K ~ iP% \ 1 ( 1 + PI P%e-»- \ 
x a ,x> -y p a_ + p a__ ) 2\P2 + iI» 1-P? ) ASA 2\P^e^ l-Pt ) A J A 



where P a = (P£,P£,P?) = (P$ cos (j) Sa , sin <j) Sa , P^) is the quark polarization vector and <p Sa its 
azimuthal angle, in the quark helicity frame (to which the x,y,z axes refer). In specific calculations, 
when performing the sum over the helicity indices A a , X' a , and A 6 , X' b in Eq. (I3"3"j) . one obtains terms of 
the form (j = x, y, z) 

Pj fa/A,S A = Isj/Sa ~ f-sj/s A = ^fs 3 /s A ' ( 40 ) 

where, as often in the sequel, we do not show the functional dependence of the /'s on x and k±. These 
terms refer to polarized quarks inside hadron A. Additionally, we can have also terms corresponding to 
unpolarized quarks inside polarized hadrons (the Sivers asymmetry) 

Af a /s A = fa/S A - fa/-S A ■ 

Parity and rotational invariance allows this function to be nonzero for transversely polarized hadrons only. 

By employing Eq. (|34l) . the relations among the functions in Eq. (138]) and those defined in Eqs. (1401) . 
(|4T|) can be more precisely cast in the following form: 



fa/A 


— fa/A,S L ~ 


- {ptt + *vt) 


fa/A,S T 


1 

= J a/ A + - 


Af a/ s T = + F--) + 2 ImF++ sin(0 5A - a ) 


Px fa/A,S L 


= Af Sx/ s L 


= 2 ReF++ 


Px fa/A,S T 


= Af Sx/ s T 


= (F+r + cos(0 5 . 4 - <f> a ) 


Py fa/A,S L 


= Py fa/ A z 


= Af Sy/SL =Af Sv/A = -2lmF+- 


Py fa/ A, S T 


= Af Sy /s T 


= Af Sy/A + A-f Sy/ST = -2 ImF+- + (f+ = - F~+) 


Pz fa/A,S L 


= ^fs z /S L 


= {FU - F^) 


Pz fa/A,S T 


= Af Sz/ s T 


= 2ReF++cos(0 5A -0a), 



(42) 

sin(0 SA - (p a ) 



where <fis A and cj) a are the azimuthal angles of the hadron polarization vector and of the quark transverse 
momentum k± a respectively, in the hadronic cm. frame. Notice that P£f a /A = ^fs x /A = 0, and that 
the function Af Sy /s T can be decomposed into two terms, one independent of the hadron transverse 
polarization, Af s /a, and a term that changes sign when the hadron polarization is reversed, A~ f Sy /s T = 
(1/2)[A/ S /5 T — Af s y /-s T ]- Analogous relations hold for quark b inside hadron B. 

As already said, apart from the well-known PDF's surviving in the collinear configuration, two of 
the remaining functions are of particular phenomenological relevance for the study of azimuthal and 
single spin asymmetries: the Sivers function, which may be responsible for several transverse single 
spin asymmetries; the Boer-Mulders function Af s /a that, besides contributing with additional terms 
to the transverse SSA's, may be responsible for azimuthal asymmetries in unpolarized processes, e.g. in 
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Drell-Yan. Appendix C of Ref. |135j gives a detailed comparison of the above notation with the one of 
the Amsterdam group |lU8j . often adopted in the literature (see also Eq. fTToT) ). 

For an on-shell massless particle with spin 1 the helicity density matrix is formally very similar to 
that of a spin 1/2 particle, considered above for the quark case. Therefore, most of the discussion and 
treatment presented there can be directly translated to the gluon case. There are two main differences: 
a) gluons cannot carry any transverse spin; there is however a formal analogy between transversely 
polarized quarks and linearly polarized gluons; b) the complete specification of the gluon polarization 
state requires the introduction of a polarization tensor Ty. Below, we briefly summarize the results for 
the gluon sector, emphasizing where they differ in some respect from the quark case. 

The helicity density matrix for an on-shell gluon inside hadron A in a spin state Sa can be cast as 

a/A,S A _ i ( i + ps r x 9 - i% 9 \ _ 1 / 1 + P? irc -pl \ 

px ^ ~2{v + i y i -pi ) a -2\-pL^ i-PLc ls A { } 

Pci rc corresponds to P%, the gluon longitudinal polarization. The off-diagonal elements of Eq. (1431) 
are related to the linear polarization of the gluons in the transverse (xy) plane at an angle <p and <p + 7i/2 
w.r.t. the x axis. The x, y and z axes refer to the standard gluon helicity frame, as reached from the 
hadronic cm. frame, in which its momentum is p^ = (p, 0,0,p). Pf in can be expressed in terms of the 
parameters Tf and T 9 , which are closely related to the Stokes parameters used in classical optics; they 
play a role formally analogous to that of the x and ^/-components of the quark polarization vector in 
the quark sector. The use of the parameters Tf and T 9 makes the gluon distribution functions formally 
similar to those for the quarks and considerably simplifies the formulae for the spin asymmetries. 

The nonperturbative helicity amplitudes J r x a ,x XA -\ A (x a , k± a ) and their products Fx A x A now refer 
to the process A — >• g + X and their parity properties are then slightly different. The rest of the 
calculation goes along the same way, i.e. adopting the eight functions in Eq. (I38p . with the difference 
that the functions F+I and F^ are now purely imaginary quantities. 

Without giving the detailed expressions of the gluon TMD distributions (see Ref. |135j ) we only 
mention that the Sivers function exists also for gluons and in this approach maintains the same physical 
interpretation as for the quark case. There is also an analogue of the quark Boer-Mulders function, re- 
ferred to in the sequel as the Boer-Mulders-like gluon function, T 9 f g /A = ^f^/A- However, the physical 
interpretation of this function is more involved: it represents the difference between the probabilities 
of finding, inside an unpolarized hadron, a gluon with linear polarization in the plane containing the 
gluon and hadron three-momenta and in the direction orthogonal to this plane. Similarly, there exists 
an analogue of the quark transversity distribution, involving again linearly polarized gluons. A more 
detailed comparison with the corresponding notation of the Amsterdam group |22Uj can be found in 
Appendix C of Ref. |135] (see also Ref. [35] ) . 

4.3 Spin and TMD fragmentation functions at leading twist 

Leading twist fragmentation functions can be described within the same formalism presented for PDF's 
in the previous section. Let us define the nonperturbative helicity fragmentation amplitudes, T>x c ,x x -,x c , 
for the process c — * C + X in which parton c, coming from the elementary hard scattering process 
ab — > cd, hadronizes into a jet containing the detected final hadron C, with light-cone momentum 
fraction z and transverse momentum k±c- We then define the product of fragmentation amplitudes 

D%${z, fcxo) = Vx c ,x x ,x c ^ fc±o) %a x; a'>, k xc ) ■ (44) 

Denoting by <j)^ the azimuthal angle of hadron C momentum in the parton c helicity frame, one has 

Vx c ,x x ,x c (z, k ±c ) = Vx c ,x x ,x c (z, k ±c ) e iX ^c . (45) 
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In the sequel, apart from a short discussion of transverse A polarization in sections I4.3.2[ 15.1.31 we 
will always refer to unpolarized final hadrons. In this case, the generalized FF D x c x ,°(z, kj_c) simplifies 
to 

D°% c (z, k ±c ) = E D X x ±(z, k ±c ) = D°&{z, k ±c ) e *CV-AU*g , ( 46 ) 

and its azimuthal independent term fulfills the simplified parity relations 

D^_ x ,(z,k ±c ) = (-1) 2 * (-1)^ D°%(z,k ±c ) , (47) 
which again imply some differences between the quark (s c = 1/2) and the gluon (s c = 1) case. 

4.3.1 Quark and gluon fragmentation into an unpolarized hadron 

From parity and rotational constraints we can easily see that for quark fragmentation into an unpolarized 
hadron there are only two independent fragmentation functions at leading twist: 

D%£{z, k ±c ) = D°l q (z, k ±c ) = D c/q (z, k ±c ) (48) 
D^(z, k xc ) = D£to(z, k ±c ) e^c=- D^(z, k ±c ) e ± ^ , (49) 

with 

[D?L q (z, k ±c ) }* = -D^ q (z, k ±c ) . (50) 
Dc/ q {z, kj_c) is the unpolarized TMD fragmentation function. When integrated over the intrinsic trans- 
verse momentum it gives the usual unpolarized FF in the collinear configuration, Dc/ q (z), 

Dc, q {z) = \ E/ d 2 k ±c D?X(z,k ±c ) = \ E / d 2 k ±c D X x ^(z,k LC ). (51) 

D+^ q (z, k± c ) is a purely imaginary quantity, related to the Collins FF, A N D c / q t (z, k± c ), 

- 2 % D^ 9 (z, k ±c ) = 2 lmD?L 9 (z, k ±c ) = A N D c/q i (z, k ±c ) , (52) 

which gives the difference between the number densities of unpolarized hadrons C produced in the 
fragmentation of a quark polarized transversely to the hadronization plane (containing the quark, p q , 
and hadron, pc, three-momenta) in the p q x Pc direction and in the opposite one. 
We can summarize these results in a form similar to Eq. ([421 : 



D c/q = D c/q ,s L = D^l q AD c/q , ST = -2 ImD°L q sin(0 S9 - <%) , (53) 

where <p Sq and, as already explained, 4>q are the azimuthal angles of the fragmenting quark polarization 
vector and of k±c m the quark helicity frame, as reached from the hadronic cm. frame. 

The gluon generalized fragmentation functions D±± (z, k±c) with equal helicity indices obey the 



same relations as for the quark case, see Eq. (1481) with q — > g, whereas, due to the different parton spins 
and parity properties, the gluon analogues of Eqs. (jl9]l . fl50|) for unequal helicity indices read 

D% 9 {z, k ±c ) = D^ 9 (z, k ±c )e^ = D%\z, k ±c ) e^" 

[D°L 9 {zMc)\ = D^ 9 (z,k ±c ). (54) 

Therefore, D+[_ 9 (z,kj_c) is now a purely real quantity. Similarly to what happens for the gluon 
parton distributions, D+^J?(z, k±c) is related to the fragmentation process of a linearly polarized gluon 
into an unpolarized hadron. In analogy to Eq. (|5"2"1) we define for gluons a Collins-like FF: 

2 D^( z , k ±c ) = 2 ReD?L 9 (z, k ±c ) = A n D c/t b(z, k ±c ) . (55) 

In the gluon helicity frame, in which the fragmentation process takes place in the (xz) plane, this 
quantity gives the difference between the number densities of unpolarized hadrons C resulting from the 
fragmentation of gluons linearly polarized along the x-direction and along the y-direction. 
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4.3.2 Quark fragmentation into a spin 1/2 hadron 

In this section we shortly describe the helicity formalism for the polarized TMD fragmentation functions 
for the production of spin 1/2 particles. We will limit to the case of quark partons. Clearly, in this 
case there is a strong analogy with the polarized TMD distribution functions for parton quarks inside 
a spin 1/2 hadron, discussed in section I4~2l The extension to the case of gluon fragmentation can be 
easily performed along the same lines, using the results presented in the previous sections. 

The polarization state of a spin 1/2 hadron C, produced in the fragmentation of an on-shell massless 
quark is described by its helicity density matrix 

Px c „x> c D c/ q , Sq (z, k ±c ) = £ D X X °J° (z, k ±c ) , (56) 

and depends on the polarization state of the fragmenting quark, through its own helicity density ma- 
trix p 9 \ S \i , and on the nonperturbative dynamics of the polarized hadronization process, through the 
generalized fragmentation functions introduced in the previous sections. The polarization state, and 
the helicity density matrix, of the fragmenting quark depend in turn on its production process. 

Analogously to the PDF sector, by using rotational and parity invariance we can easily see that 
at leading twist there are eight independent polarized TMD fragmentation functions, which can be 
connected to the generalized functions 

D++, D--, D+Z, D~t, D+t, D-Z, D%~ , DZZ , (57) 

through relations very similar to those of Eqs. f|4*2l) for the PDF case. 

Even if we will discuss spin 1/2 particle production only marginally, let us stress the following points: 

1) There are only three TMD fragmentation functions that do not vanish under k±c integration, that is 
the three usual FF's of the collinear approach: the unpolarized, longitudinally polarized and transversely 
polarized FF's (like in the PDF case, the last one has two TMD contributions; only the dominant one 
survives under k±c integration); 

2) There is a term corresponding to the Collins effect, describing the fragmentation of a transversely 
polarized quark into an unpolarized spin 1/2 hadron, as already discussed in section |4.3. 1( 

3) There is a very interesting new distribution, ADg Y / q = A N D c i/ q , related to the fragmentation of 
an unpolarized quark into a transversely polarized (in its helicity frame) spin 1/2 hadron. This TMD 
function, first introduced in Ref. [126j . needs not to vanish, and is in a sense the analogue of the Sivers 
PDF in the fragmentation sector. As it will be discussed in section [51 it can play a relevant role for the 
description of transverse A polarization in unpolarized hadronic collisions within a pQCD approach. Its 
phenomenological relevance was pointed out for the first time in Ref. |137j . where it was also named, 
given its role, as polarizing fragmentation function. 



4.4 Helicity amplitudes for the elementary process ab — ► cd 

Another essential ingredient of the master formula, Eq. (l3"3"|) . are the helicity amplitudes for the hard 
elementary process ab — > cd. These enter in the combination Mx c ,x d ;\ a ,x b M x , A . A , x , , which is diagonal 
only in the helicity indices of the unobserved parton d (fragmenting into a jet). For particles a, b, c, both 
diagonal and non diagonal helicity terms play a role, in different combinations corresponding to all the 
allowed contributions to a given polarized process. These contributions will be discussed in more detail 
in section 14. 5\ where the structure of Eq. f)33p will be explicitly given for some of the terms appearing 
in the case of doubly polarized processes. Helicity amplitudes are normalized so that the unpolarized 
cross section, for a collinear collision, is given by 

dd ab->cd ii 

~ir~ ~ t&tim A ft A i m w-*i • (58) 
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It is important to keep in mind that all hadronic polarization observables are measured in the 
hadronic cm. frame. On the other hand, in the helicity formalism the polarization state of each 
particle involved in the process (including partons) is given via its helicity density matrix, defined in 
the appropriate helicity frame, as reached by the chosen overall hadronic frame. Therefore, one must 
carefully connect all helicity frames with the hadronic cm. frame. A similar connection must be found 
between the helicity amplitudes M°, given in the canonical partonic cm. frame (that is, a frame in 
which partons a and b collide head-on moving along the ±z-axis direction respectively, creating particles 
c, d, laying in the xz plane), and their expressions in the hadronic cm. frame, M, which enter the master 
formula for cross sections. While this connection is irrelevant for any process involving only unpolarized 
elementary scatterings, it plays an essential role for the polarized case, which in principle introduces 
azimuthal factors that are not Lorentz invariants. As we will see in more detail in section [51 these 
factors, in addition to those coming from the polarized PDF's and FF's discussed previously, are crucial 
for a proper understanding of each term contributing to spin and azimuthal asymmetries. 

Let us now start from the hadronic cm. frame for the process AB —>■ C+X. Neglecting for simplicity 
all masses, the hadronic four-momenta are given by 

^ = ^(1,0,0,1), ^ = ^(1,0,0,-1), p£, = (E c ,PtAPl), (59) 

with E c = Jpj< + Pl, and s = (pa + Pb) 2 - As for partons, we introduce the usual light-cone momentum 
fractions x a = p+ /p\, Xb = p£ /p~b, z = Pclvt an d transverse momenta k± a , k±b and k±c- 

Due to the simplified kinematics, in the fcj_-integrated, collinear configuration the hadronic reaction 
plane (the (XZ) cm plane of our overall cm. frame) and the partonic scattering plane coincide. The 
hadronic and partonic cm. frames are connected by a simple Lorentz boost along the direction of the 
colliding beams. Clearly, this situation does not apply to a generalized parton model including all 
kinematical TMD effects. The elementary scattering, a{jp a ) + b{pb) —> c(p c ) + d(pd), is not planar 
anymore in the generic k± configuration and takes place out of the (XZ) cm plane. 

There are two different ways to find the elementary helicity amplitudes in the hadronic frame: 

1) One can directly compute the amplitudes M in this frame, with the given kinematical configuration, 
using well known effective calculational techniques. As a result, the amplitudes will be given in terms of 
azimuthal phase factors directly measured in the hadronic frame, having therefore a simple interpreta- 
tion. However, they will have complicated, not very illuminating, expressions. Moreover, their general 
properties, e.g. under parity transformations, will be not so evident. 

2) One can alternatively use the amplitudes M° in the canonical partonic frame, that are well known 
and have very simple expressions at leading order, together with their parity properties, which are also 
well known. As a second step, one can relate the amplitudes M° to the amplitudes M by considering 
the Lorentz transformation connecting the partonic and hadronic frames and its action on the helicity 
states of the particles involved. This method, which will be used in the sequel, leads to simple expres- 
sions for the amplitudes M, related by simple azimuthal phase factors to the amplitudes M°. However, 
the azimuthal phases will be given as very involved combinations of all angular variables coming from 
the kinematics and the Lorentz transformation. 

Following this procedure, the elementary scattering amplitudes computed in the hadronic cm. sys- 
tem are related to those computed in the partonic cm. system (in the XZ plane) by: 

Mi , , , = . , . e - i ( x aU+\£,b~Ktc-~\ d £, d ) e -i[{K-\)Z*-{K-\%\ e i(K-\)4>" (qq) 

where £j, £j and 0" are phases acquired as a result of the boost and rotations [I3J I135| . 
The parity properties of the canonical cm. amplitudes M° are the usual ones [T3j : 

M°A C ,-A d; -A a ,-A b = VaV b VM-lY a+Sb - Sc - Sd (-l) iK - Xb) - {K - Xd) K^\ > ( 61 ) 
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where rji is the intrinsic parity factor for particle i. From Lorentz and rotational invariance properties [13] 
one can also obtain useful expressions relating, up to an overall phase factor, irrelevant here, the 
canonical amplitudes for processes which only differ by the exchange of the two initial partons, a «-> b, 
or of the two final partons, c <-> d: 

<X?VJ ) = KtjUM - ) z- i<K ~ Xd) = Ktjth^ - °) e- lniK - X »\ (62) 

where the scattering angle 9 is defined in the canonical partonic cm. frame. 



4.5 Kernels for doubly polarized processes 

The computation of the invariant differential cross section corresponding to any polarized hadronic 
process A(S a) + B(Sb) — ► C + X, see Eq. (I33I) . requires the evaluation and integration, for each 
elementary process ab — > cd, of the general kernel 

£(S A ,S B )" = 

PxiX* f a / A > S A ( X a, k± a ) p b l B ^ B f b/B ,S B {Xb, k ±b ) Mx c ,X d ;X a ,X b ,X d ;X' a ,X' b D^K (*> k±C ^ " ( 63 ) 

{A} 

All the ingredients of the kernels have been separately discussed in detail in the previous sections. 
We have stressed the fact that, including k± effects, the non-collinearity of soft processes, together with 
the non-planarity of the elementary scattering processes, generate non trivial azimuthal phase factors 
that must be properly taken into account in the calculations. As it will be shown in section |5l these 
phase factors have direct consequences on the behaviour of spin and azimuthal asymmetries in term of 
the kinematical variables and can be experimentally tested. 

Before giving explicit kernels, Eq. (|6"3"|) . for some representative partonic channels contributing to 
the doubly polarized process A(Sa) +B(Sb) — > C + X, we note that for on-shell, massless partons there 
are only three nonvanishing, independent helicity amplitudes for the process ab — > cd: 

M++.++ = Ml + . ++ e itpi =M 1 V W 

M_ +; _+ = M°_ + ._ + e lip2 = M 2 °e^ 2 (64) 
M_ +;+ _ = M° +;+ _e^ 3 = M 3 °e^ 3 , 

where the phases ipi, i = 1,2,3, can be obtained from Eq. (|60l) by inserting the proper values for 
the helicities Xj, j = a,b,c,d and are such that they change sign by helicity inversion. The (+) and 
(— ) subscripts correspond respectively to (±1/2) helicities for quarks and to (±1) helicities for gluons. 
For shortness, in the kernels we will not show explicitly the PDF and FF dependences on the parton 
light-cone fractions and transverse momenta. Kernels for processes related to those presented by the 
exchange a <-» b, and/or c *-* d can be easily obtained by using Eq. fl62|) and exchanging properly some 
labels. The elementary amplitudes M°, i — 1,2, 3, and their products, can be found, e.g., in Ref. [135J. 

(1) QaQb -> QcQd processes 
2Z(S A ,S B ) qa " cqd = 

{ h/s A h/ SB (iMrT + ia^i 2 + |m 3 T) + ±ft,, SA Ml /SB (IMiT - |m 2 T - |m 3 T) 

~ { A t/s A ^fs y/ s B ~ ^f: y/ s A Ait/sJ sin(^ 3 - to) ] Ml M°} D c/c (65) 
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M° M 2 ° 

Ml M D c/C 



A ft/s A h/s B sin(<^i - ip 2 + <pc) ~ A fs y /s A h/s B cosfai - <p 2 + <j>%) 

+ fa/S A A fs x /S B sin (^l - V3 + <t>c) ~ fa/S A A fs y /S B COS (fl ~ <P3 + <Pc) 

This case includes all channels involving quarks and antiquaries like e.g. q a q\ — > q c qd, q a qb —> QcQd- 
(2) qg — > qg processes 

2Z(S A ,S B )" = 



1 2 



Af q A/ 9 /s 

Jsz/S A J s z /S B 



M\ 



|M 2 °| 2 



/ 9 /5./ s /5 fl (iMi°r+iM 

A /l/5 A /s/Sb Sin (^l - ^2 + <t>c) ~ A fs y /S A fb/S B C0S(V?1 - + 0c) 



C*/g 

M° M 2 ° A N D C/C r . 



(66) 



(3) 



Notice that in this case the amplitude M3 is zero because of helicity conservation. 
99 —* 99 processes 
2Z(S A ,S B )"° = 



fg/S A fg/S B 



+ |M 2 °| 2 + |M 3 °| 2 ) 



f 0|2 



|M 2 °| 2 



|M 3 °| 2N 



+ 2 



A h 



A fljs A A flis B (l^ 
_/s A A /li/5 B + A fr 2 /s A A fr 2 /s B ) cos (^3 - <p 2 ) 
( A fr l/SA A Kls B ~ A Kls A A Kls B ) sin (^3 " 
%s A fg/8 B cos (^i - + 20^) + A/f 2/5A / ff/5fl sin(^! - <p 2 + 20^) 
A/5. A f\/s B cos(^i - ^3 + 20g) + f g/ s A A f 9 r 2/SB Bm(vi ~ ¥?3 + 20^)] M° M° } 



M 2 ° M 3 ° S> D, 



(67) 



M° M 2 ° 

M° M 3 ° I A^D 



A few comments are in order here: 

1) The explicit expressions of the polarized TMD quark distributions, for specific S A and Sb combina- 
tions, can be read from Eq. (142 j) (for gluons see Ref. [135] ). Notice that in general these functions still 
depend on the azimuthal phases (given in the hadronic frame) of the hadron polarization vector and 
of the parton transverse momentum. These angles will therefore appear in the final azimuthal phase 
factors of a specific polarized cross section; 

2) Since for the fragmentation process into an unpolarized hadron there are only two independent TMD 
functions, the unpolarized FF and the Collins function, and only the last one depends on the azimuthal 
phase 4>q, this dependence has been explicitly factored out in the azimuthal phase factors; 

3) When considering processes obtained from those listed above by the exchange a b, and/or c <-> d, 
it is important to recall that this concerns also the phases cpi, i = 1, 2, 3 already factored out from the 
amplitudes M°, see Eq. ( 16"U|) . 

4) The unintegrated kernels still depend on the parton transverse momenta fcj_ a , k±b, kj_c- As a result, 
we have seen that the partonic process is not planar and does not lay on the hadronic production plane. 
Therefore, the kernels still contain terms that, under k± integration, must vanish for some specific 
asymmetry (as it has been checked both analytically and numerically). As an example, there are terms 
corresponding to longitudinal SSA's, that are forbidden, for the overall hadronic process, by parity 
and rotational invariance. This is somehow analogous to the case of single spin azimuthal asymmetries 
measured (in the virtual photon-target proton reference frame) in SIDIS off longitudinally polarized (in 
the laboratory reference frame) protons. Since the two ref. frames are not coplanar, the longitudinal 
proton polarization vector has a transverse component, of order M/Q, in the 7* — p cm. frame. 
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To conclude, given its phenomenological interest we show, as an example, the kernel for the q a qb 
q c qd channel for the transverse polarization of spin 1/2 particles in unpolarized AB collisions. 



JV„£(o,o)* 



(6* 



r> \ fa 



I a fb/B |M<f + |M 2 T+|M 3 T +2A/? M A/* /B cos(^ 



r 0|2 



,C/c 



+ / aM A/* y/B M° M 3 ° { A-L»g c /ST cos(^ - ^ 3 + 



Afsy/A fb/B M° 2 { A-L»g c /ST cos(^ 



C*/c 
Sy/s t 
C/c 
Sy/s t 



y?2 



sin(v?i - v? 3 + ^ 



COS i 



I cos AD^ c /st 
cos - ADg£ /sT sin(v9i - y?2 



sm 



sin 



where the polarized fragmentation functions for the process c — > C + X are the analogues of the quark 
distributions discussed in section H~2l Py cm is the transverse polarization of the spin 1/2 final particle 
w.r.t. the production plane, which is the (XZ) cm plane of the hadronic cm. frame. The angle <fi projects 
the transverse polarization of the particle, as given in its helicity frame, on the transverse direction in 
the hadronic frame just defined above. Notice that cos0 = Y cm ■ (p c x k±c)- 



5 Phenomenology of single spin asymmetries 

In this section we will consider in more detail, and from a phenomenological point of view, the theoretical 
computations aiming to a description of unpolarized and single-polarized cross section data for various 
high-energy processes of interest. We will consider in particular the generalized parton model approach 
in the helicity formalism (GPM), presented in section HJ making when possible a comparison with 
analogous results obtained within the other theoretical approaches discussed in section [3J 

Following an historical perspective, we will first consider single and, to a lesser extent, double 
inclusive particle production in pp collisions. We will give an overview on how spin and k± effects play 
a role in explaining some former and well-known observations, as well as the most recent data collected 
at RHIC on SSA's. Consideration of several final particles, like pions, kaons, heavy mesons, hyperons, 
and photons, in different kinematical configurations, will guide us in showing the most interesting 
aspects of SSA's and their description in the GPM approach. Let us recall however that for the above 
processes factorization in the TMD approach is still under investigation and, for the single inclusive 
case, SSA's appear at subleading-twist level. 

We will then move to a class of processes for which the k± factorization has been proved and 
SSA's show up already at leading twist. First we will discuss the Drell-Yan process, where both the 
unpolarized and single-polarized cases present interesting features to be analyzed in the context of the 
TMD approach. Finally, e + e~ and SIDIS processes will be also addressed, with emphasis on some recent 
data collected by the Belle, HERMES and COMPASS Collaborations. We will then give a discussion 
of the latest extractions of the Sivers function, the transversity distribution and the Collins function. 

We will show how rich this phenomenology can be, focusing also on the information one can gather 
from a combined careful study of all these cases. 

In the theoretical estimates at LO accuracy the MRST01 set |221j is adopted for the collinear 
unpolarized PDF, f a / P (x). Where it appears together with the helicity PDF, A L f a / p , the GRV98 [222J 
and the GRSV2000 [223] sets or, in some cases, the MRST01 and the LSS01 [221] sets are consistently 
used. For the FF (D h/c (z)), the KKP [225] and the Kretzer (K) [226] sets are adopted. 

5.1 Unpolarized cross sections and SSA's in pp —> C + X 

The inclusive production of large pt particles in high-energy hadronic collisions has been for a long time 
a crucial testing ground for perturbative QCD. We will focus here on the role of intrinsic transverse 
momentum effects, utilizing the helicity formalism for the computation of (un)polarized cross sections. 
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5.1.1 Pion, kaon and heavy meson production 

Pion production is by far the most representative case: a relatively large sample of data is already 
available, and more are going to come in the near future; moreover, pion SSA's may be considered as 
the original motivation behind the GPM approach. 

A crucial and distinguishing ingredient of this approach is the proper inclusion of full k^ kinematics 
and its interplay with the azimuthal phase factors appearing in the expressions of the polarized cross 
sections. We fix for the pp — > ir + X inclusive process the following configuration: we adopt the pp 
cm. frame (the hadronic frame), with the Z axis along the direction of the colliding protons, and the 
pion produced in the XZ plane with positive (p n )x- For the proton moving along +Z, the transverse 
polarization is defined along the Y axis with f (J.) = ±Y [4>s A = ± 7T /2 in Eq. (j42p ]. 

Let us start with the unpolarized case. The corresponding differential cross section for inclusive pion 
production is given in terms of the kernels in Eq. (jS3"|) . £(0,0) = 1/2[E(|,0) +£(|, 0) ], one for each 
active partonic channel, ab — > cd. As discussed in Ref. |135j . the main new result of the GPM approach 
is the appearance of additional terms besides the usual contribution from the unpolarized k± dependent 
PDF's and FF's, coupled with the unpolarized partonic cross sections. This term, once integrated over 
all fcj_'s, would be the only one in the factorization scheme in collinear configuration. It is also the only 
contribution considered in previous approaches using the parton model with TMD effects but neglecting 
the new spin and k± dependent, leading-twist, PDF's and FF's [223 [22E1 E2SJ [MI]. In fact, focusing 
here on the qq — > qq channel, the additional contributions come from: a) The convolution of two (chiral- 
odd) Boer-Mulders functions with the double transverse spin asymmetry for the qq interaction; 2) The 
convolution of one Boer-Mulders function and the Collins function (again, two chiral-odd quantities) 
coupled with the partonic spin-transfer factor. Analogously, all the other partonic channels contain 
some extra pieces as compared to the naive TMD extension of the collinear factorized approach. 

The study of so many additional contributions, given in terms of new unknown distributions, might 
look as hopeless from the phenomenological point of view. It is however the very structure behind this 
scheme, and its exact kinematical formulation, to offer a way out. 
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Figure 2: Invariant differential cross section for pp — * 7r° + X at y/s = 19.4 GeV and for 
different p? values (in GeV/c) vs. xp |231j . The corresponding results in the GPM approach 
with k± effects (a) and in the collinear partonic configuration (b) are also shown. All curves 
are rescaled by a fixed i^-factor, K = 2.5. Unpolarized PDF and FF sets: MRST01-KKP. 



Once the integrations over the nonobservable intrinsic motions have been performed, the many az- 
imuthal phases appearing in the elementary interactions, due to the noncollinear configuration, lead to 
large cancelations of most of the contributions involving the new unknown TMD functions. However, 
the inclusion of k± effects in the standard unpolarized contribution in collinear configuration could be 
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relevant for a good description of the unpolarized cross section data measured in inclusive particle pro- 
duction in different kinematical configurations. These cross sections are otherwise often underestimated 
by a factor of up to one order of magnitude [218] . As discussed in Ref. [117] threshold resummation 
effects might reduce to some extent this discrepancy. 

As an example, in Fig. [2] we show the invariant differential cross section for inclusive 7r° production 
in pp collisions, at y/s = 19.4 GeV and for different pt values, as a function of Xf- In the left panel we 
plot the estimated cross sections, obtained with inclusion of k± effects, while the corresponding results 
in the usual collinear partonic configuration are plotted in the right panel (all curves are rescaled by a 
fixed /^-factor, K = 2.5). With the aim of giving a consistent treatment of both SSA's and unpolarized 
cross sections at the same level of accuracy (namely leading order pQCD calculations with leading twist 
soft functions) a large collection of data on Drell-Yan, direct photon production and pion production 
in pp, pA collisions has been considered. NLO corrections to the unpolarized cross sections, usually 
parameterized as i^-factors, having values around 2-3, are assumed to cancel, at least partly, in the 
SSA. Some of the relevant parameters of this approach, namely the Gaussian widths of fc^-dependent 
PDF's and FF's, could then be reasonably fixed. For details see Ref. |218j . 

For higher-energy kinematical configurations, like e.g. that of the RHIC-BNL experiments, the 
dominance of the standard unpolarized term, f a / p ® f b / p <S> D n / C (but with inclusion of k± effects), 
w.r.t. all other additional contributions, comes out again in a neat way. In this case, however, the 
enhancing factor coming from the k^ effects is much less important. 

In this context, in Ref. [116] it has been claimed that experimental data on pion SSA's observed 
at fixed target experiments (moderate cm. energies) cannot be explained by pQCD, since the collinear 
approach (LO and NLO) fails to reproduce the corresponding unpolarized cross sections. On the 
contrary, at high energy a good description can be obtained at NLO accuracy. As shown in Fig. [2] the 
inclusion of k± effects reduces the gap between theoretical estimates in pQCD and data on unpolarized 
cross sections even for fixed target experiments. Therefore, a consistent treatment in the GPM approach 
of both unpolarized cross sections and SSA's for pion production can be given. 

Let us now consider the numerator of A^, focusing on two partonic channels: 
(1) q a qb -> q c qd processes 



[S(T, 0) - £(|, 0)]^™ = \ A/ a/pT (l tt) fcia) /, 
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Let us recall that the above equations have been obtained by properly exploiting the phases coming 
from the scattering helicity amplitudes and the fragmentation process, while the phases entering the spin 
and TMD initial distributions, if present, are still implicitly included in their vectorial kj_ dependence. 



40 



Again we briefly comment on the qq — > qq case. In the first term of Eq. fl69l) we find the contribution 
due to the Sivers effect; in the second term we have the Collins effect coupled with the two pieces of 
the TMD transversity distribution. Notice that these contributions, already familiar from previous 
work adopting a simplified k± kinematics, enter now with a complete phase structure without any 
simplification. Two extra terms appear: the third one contains the Boer-Mulders effect coupled again 
with the TMD transversity distribution; finally, the fourth term comes from a combination of the 
Sivers, Collins and Boer-Mulders effects. Notice that in the qg — > qg case only the Sivers effect and, 
separately, the transversity + Collins effect, could be active. For the gg — > gg process one gets an 
analogous structure as in the qq — > qq case, with a corresponding number of effects (involving linear 
gluon polarizations instead of transverse quark polarizations). 

The complex structure in terms of azimuthal phases plays an even more crucial role in the numerator 
of the SSA. To show that one can saturate all unknown polarized TMD distribution functions, using 
the corresponding positivity bounds. In some cases this is certainly an overestimate: e.g., for the 
transversity distribution it violates the Soffer bound [182] . 

In Fig. [3] (left panel) we plot all different maximized contributions to for the p^p — > ir + +X process 
in the E704 experimental configuration, for which very large values of An have been measured [4"71 l4"8t 
|4*9] . Notice that, whenever different terms could combine with different signs (e.g., Sivers or Collins 
functions for different quark flavours), their contributions have been summed always with the same 
sign, in order to avoid any kind of cancelations not resulting from azimuthal factors and phase-space 
integrations. One sees that potentially the Sivers mechanism is largely dominant, some role could be 
played by the Collins mechanism, while all other contributions are negligible [135] . 




Figure 3: Left panel: Maximized contributions to A^, plotted as a function of xp, for 
p^p — > 7r + + X processes and E704 kinematics. Solid line: quark Sivers mechanism alone; 
dashed line: gluon Sivers mechanism alone; dotted line: transversity <g> Collins effect. All 
other contributions are negligible. Unpolarized PDF and FF sets: MRST01-K. Right panel: 
Values oi An vs. xf (E704 kinematics), due to the Collins mechanism alone, obtained by 
saturating all proper bounds on the unknown soft functions. Data are from Refs. [J71 HH 09] . 
(Un)polarized PDF and FF sets: (MRSTOl)LSSOl-K. 

However, a clear-cut statement on the relevance of the Collins effect in A^(p*p — ► n + X) would 
require the combined study of the SSA for both charged pions. 

In Ref. |219j a detailed study on this issue was performed, by imposing the proper bounds for the 
TMD transversity distributions [181J with the LSS01 set for the helicity PDF's. This analysis showed 
how the phases involved, when the k± kinematics is treated carefully, are crucial, leading to a strong 
suppression of the asymmetry generated by the Collins mechanism, see Fig. [3] (right panel) . 
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The role of the azimuthal phases can be checked more explicitly by performing the simple exercise 
of setting all of them to zero (somewhat an analogue of the former simplified studies with planar 
kinematics): in such expected, one obtains considerably larger results. However, the competing 

contributions from qq — > qq and qg — > qg processes (due to the opposite sign in the partonic spin transfer 
factor), still prevent to get An values as large as those observed at large xp. This result is in contrast 
with former believes, that the remarkably large SSA measured e.g. by the E704 experiment could be 
generated by the Collins mechanism alone [131] . In fact, in these former studies, only the leading k± 
effects were taken into account in a simplified planar configuration for the hard scattering. 

It is worth to mention here that the Collins effect is not completely out of reach in hadronic collisions: 
indeed it has been recently shown that by measuring the azimuthal distribution of a pion inside a jet 
inclusively produced in p^p collisions, one can select this mechanism (see Refs. [232] and |233j ). 

From Fig. [3] (left panel) one also sees that the contribution of the gluon Sivers mechanism for E704 
kinematics at negative xp could be important. Indeed, very similar results can be obtained for the 
p^p — > 7r + X process in the kinematical configuration of the proposed PAX experiment at GSI [234J , 
where at large negative xf one could gain some information on the poorly known gluon Sivers function. 

The STAR experiment at RHIC-BNL has also measured non zero values of An in p^p — > tt° + X 
processes at forward rapidity [59j 1235] . Again, in the STAR kinematical regime the Sivers mechanism 
gives the largest contribution to An- There might remain some contribution from the Collins mechanism, 
while all other contributions are negligible. In this case, however, in contrast with the lower energy 
case, at negative xp all contributions are vanishingly small. 

Let us now try to explain some of the above results, namely, the different behaviour of the Sivers 
mechanism in different kinematical configurations. The azimuthal phase factor cos0 a entering the Sivers 
function, see section [U plays a crucial role and deserves some additional comments. The only other term 
depending on a in the contribution of the Sivers effect to the numerator of An is the partonic cross 
section, in particular via the corresponding Mandelstam variable t. At large positive xf and moderately 
large pt the (average) values of t are relatively small. Therefore, the (dominant) t-channel contributions, 
proportional to 1/t 2 , depend sizably on <f) a , so that An is not necessarily suppressed. This explains 
the potentially large effect at large xp. Instead, for negative values of xp, all partonic Mandelstam 
variables are much less dependent on <p a , so that one is roughly left with the d 2 k^ a cos0 a integration 
alone, which cancels the potentially large Sivers contribution. As a consequence, the possibility of 
gaining information on the gluon Sivers distribution from the recent STAR and BRAHMS data at 
negative values of xf is frustrated. Notice that, due to the much lower values of \fs involved, this 
suppression caused by the cos 4> a dependence is much less effective for the kinematical regimes of E704 
and of the proposed PAX [234] experiment. 

On the basis of the above results it is clear that in the GPM approach the Sivers effect is the 
only mechanism which could by itself describe the large An values observed by the E704, STAR and 
BRAHMS collaborations in the forward region. 

In Ref. [218] . in order to perform numerical estimates and carry out a comparison with available 
experimental data a simple factorized form for the Sivers function was introduced: 

A N f qM (x,k x ) =2Af q s (x)f q/p (x)H(k L )^e-^, (71) 

where (3 2 = l/(k 2 L ) is a flavour-independent width parameter defining the Gaussian k± shape of the 
unpolarized PDF (as extracted from the unpolarized cross section analysis: (k\) = 0.64 (GeV/c) 2 
[218] ). Notice that the natural positivity bounds are automatically fulfilled by imposing A/" 9 and Ti to 
be separately smaller than unity in magnitude in the full x and k± range respectively (see section [532]) • 
To start with, only valence u and d quark Sivers functions (QSF) have been considered, neglecting 
(for what concerns the Sivers effect only) possible contributions from antiquarks and gluons. The fact 
that sizable values of A N set up at large positive xp somehow supports this choice. 
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In order to describe the E704 data invoking the Sivers effect alone, one needs a positive (negative) 
u (d) quark Sivers function, A N f q / p r(x, k±) [recall the opposite sign between A N f q / p t and see 
Eq. Therefore, in the subsequent fragmentation process the unfavoured flavour contributions 

(e.g. d — > 7r + ) enter with opposite sign (all other factors in the kernel being positive) w.r.t. the favoured 
ones (e.g. u — > tt + ) in the charged pion SSA's. This implies that the extraction of the quark Sivers 
functions is sensitive to the relative weight of the leading and non-leading terms in a specific set of FF's. 
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X) vs. xp, at = 19.4 GeV and fixed = 1.5 GeV/c 



(Sivers effect); unpolarized PDF and FF sets: MRST01-KKP; data are from 03 S3 09] . 
Right: A N (p^p — > 7T° + X) at = 200 GeV; data are from [63]. A comparison with the 
Sivers effect in the GPM approach (dashed line) |218j and twist-three calculations (dotted 
line) [192j is also shown. 



Under these assumptions, a good description of E704 data was obtained, see Fig. [4] (left panel). 
Although somehow expected, the most striking feature of the Sivers functions so obtained was the 
marked valence-like behaviour (that is, f q / p t(x, k±) is strongly suppressed at small x). This analysis, 
which consistently treats the full kinematics in the calculation of both the numerator and the 
denominator of A n, also confirms all the main features of the Sivers contribution and of the Sivers 
function parameterizations found in previous papers based on a simplified planar kinematics [124[ 1125] . 

The so extracted Sivers function has been adopted to give predictions at much larger energies as 
those reached at RHIC [218] . The results have been found in good agreement with the subsequent 
published data [59J both in the positive and in the negative xf region. Latest STAR results, with 
enhanced statistics, have confirmed the peculiar xp behaviour of An- In Fig. H] (right panel) we show 
the theoretical estimates obtained by adopting the Sivers functions as extracted from E704 data (dashed 
line) and compare them with preliminary data from STAR collaboration [63J. As predicted theoretically, 
sizable pion SSA's, with features similar to those found in the E704 kinematical regime, are confirmed 
experimentally even at such large cm. energies. In Fig. H] (right panel) we also show theoretical results 
in the twist-three approach (dotted line). In order to compare the two phenomenological studies some 
words are mandatory: in Ref. [192J, adopting the twist-three approach a global fit on all available E704, 
STAR and BRAHMS data has been performed, leading to a good description of SSA results (see also 
Fig. [5]); LO cross sections in collinear configuration have been adopted to compute the denominator of 
An; in order to describe both low and high energy SSA data, calculations for E704 kinematics have 
been rescaled by a factor of 2. 
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Further data for charged pion SSA's have been collected by the BRAHMS collaboration [72], covering 
a different kinematical region (most of the data is at lower p?) and at two energy values. In this case, 
adopting the parameterization of the quark Sivers functions extracted by fitting E704 data, theoretical 
results underestimate the high-energy data (dotted lines in Fig. [5]). This indicates that an extension of 
the GPM fit procedure including all data is required to better fix the Sivers function. Indeed the global 
fit based on the twist-three approach gives a better description (Fig. [5], solid lines). 
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Figure 5: Left panel: A^(p^p ^ n ± + X) at \/s = 200 GeV for two scattering angles, 2.3° 
(left) and 4° (right). Right panel: A N at y/s = 62 GeV. Dotted line: GPM, Sivers effect; 
solid line: twist-three approach. Curves are calculated at px > 1 GeV/c. Data are from 



In Fig. [6] we show preliminary STAR data on A^{p^p — > 7r° + X), as a function of pt [63] in 
different Xf bins and compare them with theoretical predictions based on the Sivers effect alone. With 
the possible exception of the lowest pr points, at small xf, An is reasonably reproduced, taking into 
account that the GPM approach (based on pQCD), should be better suited for px bigger than 1-2 
GeV/c. Let us also remind that the theoretical curves are (genuine) predictions obtained from a simple 
parameterization of the quark Sivers functions, aimed to describe the main features of the large xf 
E704 data (at much lower energy). A similar pr-behaviour is also expected in the twist-three approach. 
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Figure 6: A N {p^p — > 7r° + X) at y/s = 200 GeV as a function of px and for different xp bins. 
Data are from [63J. Predictions for the Sivers effect in the GPM approach are also shown. 



We mention here that a preliminary study adopting the parameterizations of the Sivers function 
extracted from the azimuthal asymmetries observed in SIDIS (see section I5.5.2P leads to a much better 
agreement with BRAHMS data, still preserving a fair description of STAR data |233j . 
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Let us focus now on the mid-rapidity region. As already mentioned above, in this region only the 
Sivers effect could generate a sizable SSA in the GPM approach. In particular, at the large cm. energies 
reached at RHIC, the potentially dominant contribution should come from the gluon Sivers function 
(GSF). This region is indeed covered by the PHENIX experiment, whose recent data are consistent 
with a vanishing SSA [66J. This potential gluon dominance, together with the (almost) vanishing of 
all possible contributions to An other than the Sivers effect, allow to interpret the data - showing tiny 
values of An - in terms of useful constraints (upper bounds) on the magnitude of the GSF. As pt grows, 
however, x™ m , the minimum value of the parton light-cone momentum fraction in the polarized proton, 
increases and the dominance of the gluonic channels becomes less prominent. 

For the same reason, another set of data, for comparable rapidity and px ranges but at much 
lower energy, y/s ~ 20 GeV, cannot give significant constraints on the GSF. In this case, in fact, even 
at the smallest pt values x™ m remains large and a possible mixing with quark initiated contributions 
cannot be excluded, see Fig. [7J (left panel). Different scenarios have been considered in Ref. [236], where 
a significant constraint on the GSF has been obtained, see Fig. [7] (right panel). For instance, even 
allowing for the largest possible cancelations coming from maximized, sea-quark Sivers contributions 
(still totally unknown), in the small x region (below 0.05), where the gluon plays a relevant role, the GSF 
is bounded to be less then 30% in size of its positivity bound (twice the unpolarized gluon distribution). 
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Figure 7: A N (p 1 'p — > 7r° + X) at mid-rapidity vs. px- Left panel: E704 data [51] at 
y/s ~ 20 GeV. Maximized contributions, in magnitude, from the quark (dashed line) and 
gluon (dot-dashed line) Sivers functions. Right panel: PHENIX data [66] at y/s = 200 
GeV. Contributions of: maximized GSF (dot-dashed line); sea (maximized) + valence quark 
Sivers functions (dotted line); constrained GSF (dashed line). Overall total Sivers effect 
(thick, solid line). In both plots the thin, solid lines show the predictions using the u and d 
Sivers functions extracted in Ref. [218j . Unpolarized PDF and FF sets: MRST01-KKP. 



As an indirect tool to constrain the gluon Sivers function one could also consider the Burkardt sum 
rule of Eq. ( 1231) |184j . As discussed in section [3], this rule states that the net transverse momentum 
due to the Sivers mechanism must vanish. Strict use of the BSR would require integration over the full 
x range, including the poorly known low x region, of each single parton Sivers function, which might 
even result in divergences. For this reason at the present stage it can be only adopted as a useful 
cross-check. Similar conclusions on the smallness of a possible GSF have been reached in a different 
context in Ref. [237] . Let us also mention that in Ref. [238] the BSR, together with the large N c limit of 
QCD, was used to predict a suppression of the gluon Sivers function w.r.t. the non-singlet quark Sivers 
distribution at not too small x. 

Another interesting process to study is inclusive kaon production in p^p — > K + X, in particular for 
charged kaons, where hopefully one could learn more on T-odd mechanisms in the sea-quark sector. 
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To begin with, one can figure out what SSA's would result in the valence-like picture for the quark 
Sivers function discussed above. Consider for the sake of clarity the forward region (large positive xp). 
In this case, A N (K + ) should be quite similar to A N (n + ): in both cases, a valence u quark from the 
incoming proton fragments dominantly into the observed meson. This qualitative expectation of the 
TMD generalized parton model works, in principle, for the Sivers as well as the Collins effects, and 
holds true for the twist-three collinear formalism and for the orbiting valence quark model discussed in 
section [3l However, the K~ case is more subtle, since the valence u, d, quarks from the proton are both 
subleading in the fragmentation process. Therefore, the expected SSA is more sensitive to the details 
of the approach. We already said in section [3] that in the orbiting valence quark model one would get a 
zero SSA for K~ . This would also be the case for the Collins effect in the model of Artru [128j . because 
the specific Lund string model utilized gives a vanishing Collins function for subleading fragmentation. 
On the other hand the Sivers effect, involving unpolarized TMD fragmentation functions, would give a 
positive and, depending on the FF set adopted, sizable SSA |125] . This somehow "unexpected" result 
for An(K~) comes from a non negligible unfavoured u — > K~ fragmentation function and the fact 
that, at large x, the (positive) u quark Sivers function is bigger in size than that for a d quark. This 
expectation seems to be confirmed by recent preliminary BRAHMS results (see Fig. [8]), although the 
theoretical predictions still underestimate the data. This discrepancy could be reduced by including 
the sea-quark Sivers effect, or the non negligible contributions of the Collins effect from favoured and 
unfavoured FF's. 

Similar results have been obtained in the global fit based on the twist-three approach, adopting 
chiral-even twist-three distributions [192j . See also section [5.5.21 for further discussion on the kaon case. 
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Figure 8: An for inclusive charged kaon production in pp collisions, at ^/s = 200 GeV (left) 
and y^i = 62 GeV (right) as a function of xp. Dotted line: GPM, Sivers effect; solid line: 
twist-three approach. Data are from 



Another useful tool to access specific spin and TMD effects is the study of heavy meson production 
in pp collisions. D mesons originate predominantly from the fragmentation of c or c quarks, which 
at LO can be created either via qq annihilation, qq — > cc, or via the gluon fusion process, gg — > cc. 
In this case the heavy quark mass cannot be neglected. This somehow makes the structure of the 
kernels in the helicity formalism of the GPM more involved. The main difference is that now helicity 
conservation is broken in the hard scattering amplitudes. Once again, a detailed study including the 
exact k± kinematics shows that upon k± integration the azimuthal phase factors suppress all allowed 
terms but the Sivers contribution. 

In Ref. |239] it was shown how at RHIC energies, = 200 GeV, and px roughly in the range 1 -j- 3 
GeV/c, the gluon fusion mechanism gives the dominant contribution to the p^p — > D + X process, up 
to xf — 0.6, opening the possibility for a direct access to the gluon Sivers function. At lower energies, 
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like those reachable at the proposed PAX experiment at GSI, y/s = 14 GeV, the qq annihilation process 
becomes dominant, giving the opportunity to access directly the quark Sivers function. For polarized 
pp collisions at PAX the potential QSF dominance would be even more dramatic p40j. 

Let us finally mention here that the inclusion of a possible intrinsic charm component in the proton, 
as originally proposed by Brodsky et al. |241j (see e.g. Ref. [242] and references therein), could change 
drastically this picture. In this case, in fact, c-quark initiated processes and other TMD effects could 
play a role. This interesting topic would certainly deserve a more detailed study. 

5.1.2 Prompt photon production 

The inclusive production of direct photons in pp collisions is certainly another useful tool to access 
fc_i_-dependent PDF's. At LO, only two partonic subprocesses, q(q) g —> 7 q(q) and qq — > jg, contribute. 
The polarized differential cross section can be obtained from the general case, Eq. (1331 . by replacing 
the fragmentation function by a product of two Dirac delta functions, Dc/ C (z, k±.c) - > $(z) 8(kxc)- 

There are in principle three mechanisms that could contribute to the photon SSA: the quark Sivers 
effect for both subprocesses, the gluon Sivers effect via the Compton-like subprocess, and the Boer- 
Mulders effect (coupled with the transversity PDF) via qq annihilation. Notice also that electromagnetic 
couplings will enhance the it-flavour contributions. For xp > 0, and fixed pp, the dominant contribution 
to An should come from the quark Sivers effect. Adopting the QSF parameterizations extracted from 
the fit to A N (p^p — > 7r + X), in the GPM one gets a positive SSA, rising at larger x F , both at = 20 
GeV and = 200 GeV [218J. Notice that the collinear twist-three approach [TSDl HSU H22], while 
leading to a similar description of A N for pion production, predicts a photon SSA opposite in sign. 
The reason is just the opposite sign entering the hard scattering part for the dominant channel in the 
twist-three formalism. In this respect, this process can be considered as a good tool to discriminate 
between the two approaches. 

Concerning the backward rapidity region, in Ref. [243.] it was argued that at RHIC energies and large 
Pt values (around 20 GeV/c) An would be directly sensitive to the GSF. In fact, in a more quantitative 
approach, with proper treatment of the noncollinear kinematics and the relative azimuthal phases, it 
has been found that the best kinematical region to look for the GSF at RHIC energies is at negative 
xp and pt values around 5-8 GeV/c. In this region the contribution to An from the GSF could be as 
large as 10%, whereas the other mechanisms would give at most a 1% contribution [240]. 

An even more interesting case is the photon SSA at lower energies, like for instance at the PAX-GSI 
experiment. In this case, the QSF gives the main contribution in the forward as well as in the central 
rapidity region, whereas again the negative xp region is potentially dominated by the GSF. More- 
over, by colliding polarized protons off unpolarized antiprotons, another effect could become accessible: 
namely the Boer-Mulders function, A N f q t i p (or h±), coupled to the transversity distribution. In this 
configuration the Compton-like subprocess is suppressed because: i) the minimum value of x reached 
is still large, being of the order of Ippj \fs] ii) the q in p has a valence component. As a result, we have 
a clear dominance of the qq subprocess. Moreover, the integration over the fcj_-dependent phases does 
not wash out the partonic double spin asymmetry. By using the parameterizations so far extracted for 
the QSF one gets a photon SSA of the order 5-10% for pp around 2-4 GeV/c, while the (maximized) 
Boer- Mulders effect might give An values of the order up to 20%. This means that in this kinematical 
regime a large photon SSA could be a clear signal in favour of the Boer-Mulders effect, then giving 
another way to access the transversity distribution [240] . 

5.1.3 Transverse A polarization in unpolarized pp collisions 

Another longstanding problem of spin phenomena in particle physics is the transverse polarization of 
hyperons produced in unpolarized hadronic collisions, Pp(Y), that can also be considered as a SSA. A 
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huge amount of experimental results has been collected, in particular for the A particle. These data 
show that transverse hyperon polarizations can be quite sizable. Their behaviour as a function of xp and 
Pt is very peculiar and difficult to explain, as well as the differences in size and sign among different 
hyperons. Several nonperturbative models have been proposed, but none of them is able to give a 
comprehensive description of all data. We will not address in detail this subject. A quick survey of its 
experimental and theoretical aspects can be found in Refs. [6], [8] respectively, and a more extensive 
treatment in the references quoted therein. Most of the experimental data are at fixed target and in the 
very forward regime. For this reason, a pQCD approach to Pt(Y) was not seriously considered until 
recently. In fact, the description of the unpolarized cross section data [2] is difficult in pQCD, even 
when, as in the GPM approach, full k± effects are included (notice that also in this case the leading 
term in the unpolarized cross section is given by the usual convolution of TMD unpolarized PDF's 
and FF's). More precisely, data on unpolarized cross sections can be described only partially and in 
a limited kinematical range, corresponding to not too low pp and not too large xp. Limiting to the A 
case, for fixed target experiments the kinematical regions where data for the unpolarized cross section 
and for the A polarization are collected overlap only partially. A better agreement is found with the 
recent STAR data on unpolarized cross sections at larger energies [248J. 

In the GPM approach, focusing on the qq — > qq processes, there are essentially two possible sources 
for the A polarization (see also the end of section H]): the Boer- Mulders function from one of the incoming 
protons convoluted with a spin transfer partonic factor and the transverse fragmentation function (the 
analogue, in the fragmentation sector, of the transversity function); the so-called (chiral-even) polarizing 
fragmentation function (A D A t / q ), coupled with unpolarized PDF's and partonic cross sections. This 
last mechanism is somewhat the analogue of the Sivers effect in the fragmentation process. An ongoing 
preliminary analysis [217] seems to show that indeed only the polarizing effect can lead to sizable 
values of Pp{A). This study basically confirms the main findings of Ref. |137] where, adopting a more 
simplified planar kinematics and a simple parameterization of the polarizing FF, the main features of 
Pt(A) (including its size) were reproduced (see Fig. [9]). 
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Figure 9: Best fit |137j to -Pt(A) data from p-Be reactions vs. pp and for different xp 
bins. Data |244l I245[ I246] are collected at two different cm. energies, ^snn — 27 GeV and 
\J s nn — 39 GeV. Theoretical curves are evaluated at the mean xp value in the bins, and at 
^snn = 27 GeV; the results change very little with the energy. Unpolarized A FF set: [247] . 
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5.2 SSA's for double inclusive production in pp collisions 



Another important source of information on TMD distributions could come from the study of pp — > 
hi h 2 + X and pp — > jet jet + X processes. Indeed, by considering specific kinematical configurations 
they have been advocated as a tool to extract the gluon Sivers function. In Ref. [249] they considered 
the almost back-to-back asymmetric jet correlation in the mid-rapidity region in the GPM approach 
and showed that a sizable asymmetry could be observed at RHIC, mainly driven by the gluon Sivers 
effect. In Refs. [1541 1164] a more complete treatment was presented, including gauge-links effects and 
all possible spin and k± dependent distributions at leading twist for double inclusive hadron and jet 
production in pp collisions. Recent data on jet-jet correlations have been collected at RHIC showing 
an almost vanishing asymmetry [62]. A possible explanation of these results could come from the 
gauge-links effects properly incorporated in the TMD approach [164] . 

A recent phenomenological study for the SSA in pp — > 7 jet+X [163] (see also [25U] ) has addressed the 
comparison between the approach based on the generalized parton model and the colour-gauge-invariant 
QCD formalism, where the gauge links are believed to be responsible for initial/final-state interactions 
(in other words, to be the source of the phase shifts required to generate naively T-odd effects). For 
hadronic processes different from SIDIS and DY, for instance p^p — > hadrons, the Wilson line structure 
becomes more intricate [153] . It is therefore more challenging to derive clear-cut predictions for the sign 
of the SSA in these processes [164] . Again, in direct photon production the reduced number of channels 
and the colourless structure of the photon greatly simplify the analysis. 

Due to the intrinsic parton motions in the initial hadrons, in the process pp — > 7 jet + X the photon 
and the jet are only approximately back-to-back in the plane transverse to the beam collision direction. 
In other words, the azimuthal angles of the photon and jet transverse momenta, and Kj±, are 

related by <p- s — 7 = it + 50, where 5<p just originates from k± effects. By defining suitable weighted 
azimuthal moments of Ky±/M, sensitive to these effects, in Ref. [163] it has been shown how one could 
check the role of the gauge-links, that in this simple process should only appear as numerical prefactors 
for the standard LO partonic cross sections involved. 

In a simplified partonic kinematics, keeping only 0(k±/K^^) effects, and suitably choosing the 
kinematical configuration (that is, large positive values of the photon pseudorapidity, r/ 7 , small or 
negative values for the jet one, r/j) it has been shown that for RHIC energies the quark-Sivers mechanism 
dominates. Moreover, the weighted moment of interest could reach values of 4-5% in size. Notice that 
fixing the sign of the valence quark Sivers functions from SIDIS, the GPM approach would predict an 
opposite sign w.r.t. the colour-gauge-invariant QCD formalism. 

The experimental investigation of the sign of the weighted moment, discriminating among these 
predictions, seems to be accessible at RHIC and could have the same significance as the measurement 
of the relative sign of the Sivers effect in SIDIS and DY processes. However, a word of caution is needed 
here and for the double inclusive processes discussed above, being the validity of factorization for these 
processes still under debate, see Refs. [1601 11611 1162] and section 13.1.11 



5.3 Azimuthal asymmetries and SSA's in Drell-Yan processes 

The Drell-Yan process is definitely another valuable tool to access naively T-odd TMD effects in the 
PDF sector and discriminate among them. Starting already from the unpolarized cross section, we can 
learn on k^ effects in PDF's by studying the spectrum of the final lepton pair [218] . But there is much 
more to learn from the DY process. Indeed the unpolarized differential cross section, as measured in 
the so-called Collins-Soper frame [120] for the process N — > /i + fi~ + X, where N is either deuterium 
or tungsten, using a n~ beam with various energies in the range of 150-300 GeV [1021 11031 1106] . 

--^ = T~ T~ — ( l + A cos 2 # + u sin 26 1 cos </> + - sin 2 6» cos 20 V (72) 
a d\ I 47r A + 3 V 2 / 
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shows remarkably large values of v (see Fig. [10] for the proper definition of the angles). 




lepton plane (cm) 



Figure 10: Kinematics of the Drell-Yan process in the lepton center of mass frame. 

Leading and next-to-leading order perturbative QCD corrections cannot explain these results. Sev- 
eral explanations have been proposed, e.g. higher twist effects [2511 1252] and the role of QCD-vacuum 
structure in hadron-hadron scattering [2531 12541 1255j . 

In the context of a pQCD approach with TMD distribution functions, such a large cos 20 azimuthal 
dependence can arise at leading order only from a convolution of two Boer-Mulders functions, A N f Sy /A 
or h^. In Ref. |136] it was shown how by a suitable parametrization of this naively T-odd function, under 
some simple but reasonable assumptions a description of the cos 20 asymmetry could be obtained, see 
Fig. [TT] (left panel). More recently, the cos 20 asymmetry in the unpolarized DY process pp — > fi + fi~X at 
s = 50 GeV 2 , has been studied in detail in the same framework but using a more realistic model for the 
Boer-Mulders function in the quark-diquark spectator approach [256] . Notice that the negligible cos 20 
dependence observed recently for DY dimuons in pD collisions [107] could imply an almost vanishing 
sea-quark Boer-Mulders function. 

Considering the transverse SSA in the Drell-Yan process p^p — > £ + £~ + X, again two TMD mecha- 
nisms could play a role: the Sivers effect and the Boer- Mulders effect, which involves also the transver- 
sity distribution. The main difference with respect to the inclusive processes considered above, like 
pp — > h + X, is that in this case the measurement of the lepton-pair angular distribution automatically 
allows to select one specific effect. In Ref. |136j the full structure of the polarized cross section, at 
0(k±/Mh), has been presented. From there, one can easily see that integrating over the lepton-plane 
angular variables, only the Sivers effect survives in the kinematical regime where k± ~ qx <C M (M here 
is the invariant mass of the lepton pair). Notice that other mechanisms that could also generate SSA's 
in the Drell-Yan process |257l 1258] 259, 260], based on higher twist quark-gluon correlation functions in 
a generalized pQCD factorization theorem [191] . lead to expressions of which, being dependent on 
the angle between the proton polarization direction and the plane containing the final lepton pair [120J, 
vanish upon corresponding angular integrations. In Ref. [261] . considering in the hadronic cm. frame 
the polarized differential cross section da/dM 2 dyd 2 qT, where y is the rapidity of the lepton pair, a 
phenomenological analysis of the Sivers asymmetry has been performed in great detail. By adopting 
the Sivers functions as extracted from pion production data (with the reversed sign, see below) prelim- 
inary predictions for at RHIC were given, see Fig. [TT] (right panel). A factorized Gaussian ansatz 
for the fc^-dependent PDF's, with some reasonable approximations, allowed to perform analytically the 
integration over the transverse parton momenta, resulting in a simple factorized expression for A^. 

Let us comment on the sign of the SSA originating from the Sivers effect in the Drell-Yan process. 
As already discussed in section (3] if the Sivers effect is due to initial/final state interactions and colour- 
gauge invariance (Wilson lines) the corresponding asymmetry should have opposite sign in Drell-Yan 
and SIDIS processes, respectively because of past and future-pointing gauge-links, related, somehow, to 
s-channel and t-channel elementary reactions. On the other hand, the SSA in the p^p — »• n + X process 
at large positive xf and moderately large pr is dominated by t-channel, quark initiated processes. 
Therefore, one can tentatively assume that the Sivers function as extracted from pion SSA data should 
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Figure 11: Left panel: Data on the v coefficient in Eq. (E2) [1021 H03] at 194 GeV, vs. 
qr, compared with the fit of Ref. |136j . Right panel: Maximized An for the DY process, 
at RHIC energies, y^s = 200 GeV, as a function of the rapidity y and for two ranges of 
the invariant mass. The asymmetry is practically negligible in the range y < 0. PDF set: 
MRST01. 



contribute with opposite sign to DY processes. Notice that this assumption is not a prediction of 
the GPM approach, which would consider the Sivers function as universal, but rather an attempt to 
incorporate in it additional information coming from different approaches not yet directly applicable to 
the p^p —> it + X process. Given these considerations, even a simple comparison of the sign of these 
estimates with data might be highly significant. 

Notice also that, as we will see below, the valence-quark Sivers functions as extracted by SIDIS 
data have the same sign as those extracted from pp — > it + X: this somehow supports the sign change 
adopted for DY in this preliminary study. 

Predictions for SSA's in DY processes based on the extraction of the Sivers effect from SIDIS |262j 
(see also section I5.5.2P have been discussed in Refs. |238l 12501 12631 1264j . As pointed out in Ref. [264] . 
the Sivers asymmetry at RHIC should be strongly sensitive to the sea-quark Sivers distributions. 

We finally stress that the study of Ref. [261j included a combined description, in the GPM approach, 
of the unpolarized Drell-Yan cross sections. 

Another promising idea in the study of SSA's has been proposed in Ref. |265j . where they considered 
electroweak Drell-Yan processes. In this case by looking at a vg — i lepton pair one selects a VT-exchange 
channel and therefore only one parton helicity state contributes. The transversity function that could 
couple to the Boer-Mulders effect is therefore absent. This process would then give another way to 
directly access the Sivers effect. 

5.4 e + e annihilation in two nearly back-to-back hadrons 

Hadron production in e + e~ collisions would be in principle the most valuable and clean process for 
the study of TMD polarized fragmentation functions, like the Collins function, thanks to the lack of 
corresponding TMD effects in the initial state. Unfortunately, there is not in the process e + e~ — > qq 
any transverse polarization transfer to a single, on-shell massless quark. Therefore, it is not possible to 
observe the Collins effect individually by measuring, e.g., the asymmetry in the distribution around the 
jet trust axis (i.e. the fragmenting quark direction) of a hadron species (typically pions) produced in 
the quark fragmentation. Indeed, in hadron production in e + e~ — > qq — > 2 jets events, the Collins effect 
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can only be observed when the quark and the antiquark are considered simultaneously. As reported in 
section I2.3.3[ the Belle Collaboration at the KEK-B asymmetric-energy e + e~ storage rings has in fact 
performed a measurement of azimuthal asymmetries in hadron-hadron correlations for inclusive charged 
dihadron production, e + e~ — > ir + ir~+X |109j . This asymmetry has been interpreted as a direct measure 
of the Collins effect, involving however (and unavoidably) the convolution of two Collins functions. Two 
experimental methods have been used to extract the azimuthal asymmetry from raw data: 1) The first 
method, requiring the knowledge of the two-jet thrust axis, corresponds to the kinematical setup shown 
in the left plot of Fig. [12] and gives rise to a cos(0i + fa) modulation in the dihadron yields, the A\ 2 
asymmetry discussed below and shown in Fig. [13] [fa, <p2 are the azimuthal angles of the two hadron 
momenta); 2) The second method, leading to the asymmetry Aq (see Fig. [T3l) . does not rely on the 
knowledge of the thrust axis and corresponds to the setup shown in the right plot of Fig. [T2l 




Figure 12: Kinematics for the e + e — > h%h 2 + X process, according to two different config- 
urations considered by Belle and discussed in the text: setup 1 (left) and setup 2 (right). 



On the theoretical side this azimuthal asymmetries have been studied in Ref. [266], where their 
phenomenological interest was pointed out, and more recently in Refs. [2671 [268], where Belle results 
were used to extract additional information on the Collins function. Both Refs. [2661 126 7j have used 
in their analysis the setup 2. In Ref. |267j . utilizing the formalism of Ref. |266j and a quark soliton 
model for the transversity distribution, the Collins function was extracted by fitting the HERMES and 
COMPASS data on azimuthal spin asymmetries in the SIDIS process lp^ — > I'h + X (see next section). 
The parameterization extracted for the Collins function was shown to be compatible with Belle results 
on e + e~ — > 7r + 7r~ + X. In Ref. |268] both kinematical setups were considered, and SIDIS and e + e~ data 
were combined to get a simultaneous fit of the Collins and the transversity distribution functions. 

As an illustration, let us briefly summarize the case of setup 1, which is more simple to figure 
out. Within the GPM approach and the helicity formalism, the differential cross section for dihadron 
production can be written as (see Fig. [12] for notation): 

^l e l + COs2 9) D h!/q(Zl,P±l) D h2 /q(z 2 ,P± 2 ) 

AS q I 

+ - sin 2 6 A N D hl/q i(z 1 ,p ±1 ) A N D h2/ ^(z 2 ,Px2) cos0icos0 2 |- (73) 

The azimuthal asymmetry, as defined by the Belle collaboration (setup 1), can be obtained from here 
by the following steps: i) Performing a change of angular variables from {<pi^fa) to (0i,0i + fa) and 
then integrating over the moduli of the intrinsic transverse momenta, p±\, p± 2 , and over the azimuthal 
angle fa; ii) Normalizing the result to the azimuthally averaged cross section; Hi) Taking the ratio R 



dzi dz 2 d 2 p±i d 2 p± 2 d cos 9 
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Figure 13: Belle data on two azimuthal correlations in e + e~ — > nn + X processes |109j . Solid 
lines: results obtained by fitting the A\ 2 asymmetry |268j : dashed lines: results obtained by 
fitting the Aq asymmetry, according to the setup of Fig. [12] (right panel) [268] . The shaded 
area corresponds to the theoretical uncertainty on the parameters. FF set: K. 



of unlike-sign (U) to like-sign (L) pion-pair production, in order to eliminate false asymmetries: 

R ~ 1 + ccxifa + fa) A 12 (zi, z 2 ) , where A 12 {z u z 2 ) = \ ^ m %- (P v - P L ) , (74) 

4(1 + cos z U) 

the angle 9 is averaged over a range of values given by the detector acceptance, 

= E q e 2 q A N D n+/g ,(z 1 )A N D 7T - (+)/qr (z 2 ) 
U{L) ^qPi+iMDrMffa) ' 

and 

A N D h/q ^z) = f d 2 p L A N D h/g ,(z, P± ) = f d" P± ^- Ht q (z, Px ) = 4 H^ m {z) . (76) 

In this last equation, the relationship among the involved moments of the Collins function, as defined 
in Refs. [2671 1268] . has been also given. As an example, in Fig. [131 we show the quality of the fit 
to the asymmetry, as defined by the Belle Collaboration for the two experimental setups, according 
to Ref. [268] . This allowed to extract the first independent parameterization of the Collins function, 
that was then used to extract, in combination with HERMES and COMPASS results on azimuthal 
asymmetries, the u and d quark transversity distribution function (see also section 1531) . 

Let us finally remind that preliminary Belle data using the unlike-sign over the all charged (C) 
pion-pair ratio [110] have been found in agreement with the predictions of Ref. [267J. 

5.5 Semi- inclusive deeply inelastic scattering 

Semi-inclusive particle production in deeply inelastic lepton proton scattering, £p — > £'h + X, is a rich 
case of study over many points of view. By focusing once again on the role of intrinsic parton motion, 
both in the unpolarized and in the single polarized case, we will show how a consistent picture of various 
asymmetries observed experimentally can be given. 
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At variance with inclusive particle production in proton-proton collisions, the standard reference 
frame adopted to deal with SIDIS scattering is not the lepton-proton cm. frame, but rather the virtual 
boson-proton cm. one, for the process 7*p — ► h + X. Therefore, one cannot adapt in a straightforward 
way the kernel structure discussed in section H]to this case. Indeed, a complete treatment of full k± effects 
in SIDIS within the generalized parton model and the helicity formalism is still in progress |269j . What 
is currently adopted by many groups, using also different approaches, is a first order approximation in a 
(k±/Q) power expansion, where Q is the virtuality of the exchanged boson. In this context, a complete 
classification for SIDIS cross sections in the TMD approach has recently appeared [87] . 

As in Drell-Yan processes, in SIDIS the intrinsic parton motions result in azimuthally dependent 
terms in the unpolarized and single polarized cross sections for hadron production [270[ 12714 I272[ 1931 192] . 

Let us first of all set the kinematics that we will adopt in the rest of this section. As already said above 
we work in the 7*p cm. frame, shown in Fig.[TH We take the photon and the proton colliding along the z 
axis with momenta q and P respectively; following the so-called "Trento conventions" [75], the leptonic 
plane coincides with the x-z plane. We consider the kinematic regime in which Pt ~ Aqcd — k^, where 
Pt = \Pt\ is the final hadron transverse momentum. In this region QCD factorization holds |122] 
and leading order elementary processes, £q — > £q, are dominating: the soft Pt of the detected hadron 
is mainly originating from quark intrinsic motion rather than from higher order pQCD interactions, 
which, instead, would dominantly produce large Pt hadrons [2731 12741 1275] . 



5.5.1 Azimuthal dependence in unpolarized SIDIS 

In the QCD factorization scheme the SIDIS cross section for hadron production reads 

dx B dQ* dz h d 2 P T =^ e «J dkl - Afrfo *-0 4^ -qT- D ^P±) Y h Y{ 1+ xtQ-z) ' (77) 
where the elementary Mandelstam variables are given by (£ being the incoming lepton three-momentum) 



xs -2£-k ± -kl^[l-^£) u = - x (s-^)+2i-k L -kl^ • (7S 




The on-shell condition for the final quark implies x = x-q 1 + y 1 + 4kj_/Q 2 /2. 

Notice that p± is the transverse momentum of the hadron h w.r.t. the direction of motion of the 
fragmenting quark, and z is the light-cone momentum fraction carried by the produced hadron. After 
some algebra one can find an exact expression, at all orders in (k±/Q), for the light-cone fragmentation 
variables z and p±, in terms of the usual SIDIS and hadronic variables x&, Q 2 , Pt, %h — {P ' Ph) /(P'O)- 

Due to these relations among the kinematical variables, momentum conservation and the integration 
over intrinsic k± at fixed Pt = Pr{cos4>h, sin <frh, 0) result in a dependence on the azimuthal angle of 
the produced hadron, <ph, that is the angle between the leptonic and the hadronic planes, see Fig. [HJ 
This azimuthal dependence also appears in the SIDIS cross section. To see this, we will adopt for the 
k± dependent PDF's and FF's the usual factorized, flavour independent Gaussian forms (this helps in 
finding analytical expressions), that is: 

f q / P (x, k ± ) = f q/p (x) - JL. e~ k V^ , D h/q (z, P± ) = D h/q (z) ^ e^/ft) . (79) 

To have a closer look at the azimuthal dependence of the cross section it is instructive, and often quite 
accurate, to consider Eq. (177]) at the order 0(k±/Q). In such a case x ~ Xb, z ~ Zh and 



d 5 a i P ^£h+x 2ira 2 e 2 

m q fi( XB "> D h/ q ( z h) 



dxs dQ 2 dzhd 2 PT q Q 4 



,2 A (2-y)y/T^(kl)z h P T 

y) " 4 Wq cos0/l 



p|/H> ) (80 ) 
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Figure 14: Kinematics of the SIDIS process in the 7* — p cm. frame. 



where 

(Pt) = (Pi) + 4(kl) ■ (81) 

This approximate result illustrates very clearly the origin of the dependence of the unpolarized SIDIS 
cross section on the azimuthal angle <ph- As first observed by Cahn |270j . such a dependence is directly 
related to the parton intrinsic motion and vanishes when k± = 0. 




01234560123456 



i» h (radians) o h (radians) 

Figure 15: Fits to data [92] on the cos (fih dependence of the unpolarized SIDIS cross section 
for charged hadron production. Dashed line: exact kinematics; solid bold line: 0(k±/Q) 
calculation. The shadowed region corresponds to varying the parameters (k\) and (p\), 
given in Eq. flM} of the text, by 20%. Unpolarized PDF and FF sets: MRST01-K. 

Indeed, this effect (also known as Cahn effect) has been measured in different SIDIS experiments [921 
1931 [94] . In the TMD approach the available data can be used to extract information on the Gaussian 
widths entering the PDF's and FF's. A systematic analysis, using several sets of experimental data 
(see Fig. [15] as an example) has been carried out in Ref . [276] , leading to the following best values of 
the parameters: 

(hi) = 0.25 (GeV/c) 2 (p 2 ± ) = 0.20 (GeV/c) 2 . (82) 

This study leads to an excellent agreement with data for small values of the final hadron transverse 
momentum, P T . However, as somehow expected, it fails in reproducing data at higher P T , the turning 
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point being around Pt ~ 1 GeV/c. A similar conclusion was drawn in Ref. |273j . As discussed in 
Ref. |277j . by including also NLO corrections (like hard gluonic radiation and gluon-initiated elementary 
scatterings) the agreement with data can be extended at much larger Pt values. This analysis also 
confirms that the transition between the regime dominated by nonperturbative, intrinsic transverse 
momentum effects and the one mainly governed by NLO radiative effects sits around Pt — 1 GeV/c. 

We also mention that preliminary results |278j in lattice QCD for the first x-moment of the unpo- 
larized TMD distribution are compatible with a Gaussian k± dependence with (k^_) ~ 0.3 (GeV/c) 2 . 

In the context of azimuthal asymmetries in unpolarized SIDIS there could be also a cos 20^ depen- 
dence. Indeed a contribution to this term is already present in Eq. (1771) . due to kinematical (k±/Q) 2 
effects. Moreover, in the TMD approach including spin effects, an additional contribution appears, 
coming from the convolution of the Boer-Mulders and the Collins functions |108[ I168[ I279[ [87]. 



5.5.2 Transverse SSA's in SIDIS 

Let us now consider the single-transverse polarized cross section da[£p^ — > £'h + X), or daur- Retaining 
only leading terms in 0(k±/Q), in the TMD approach at leading twist the differential cross section can 
be formally expressed as 

^0 + o [^Sivers + ^Collins + ^Collins] ■ (83) 



dxs dQ 2 dzh d?PT d<ps 2 

The <ps dependence originates from the cross section dependence on the angle between the proton 
(transverse) polarization vector and the leptonic plane; in the configuration of Fig. [TH this angle is 
simply <frs — <f>£> = <fis, having chosen <j)g/ to be zero, with <ps varying event by event. 

The A term in Eq. (153"]) is, apart from a factor of 2tt coming from the <fts integration, the unpolarized 
cross section given in Eq. f)77p and, in the consistent 0(k±/Q) approximation, in Eq. flSD]) . The terms 
^sivers and A Collins in Eq. (|83|) correspond respectively to the Sivers effect and to the Collins effect 
(involving the k± dependent transversity distribution, ATq{x, k±)) and are related to the A s J^^ h± ^ s ^ 
asymmetries discussed in detail below. The fourth term, ^4'c iii ns , involving again the Collins function 
coupled with the h^ T distribution (in the notation of the Amsterdam group) provides eventually a 
sin(30h — 4>s) azimuthal modulation. Only recent preliminary results from COMPASS are available on 
this observable [86J. The measured asymmetry is compatible with zero within the statistical errors. We 
will not consider this contribution in the following fenomenological analysis. 

Exploiting the helicity formalism, the two interesting additional terms, Ag ivers and Acoiims, read 

^sivers = J2 e lJ d2k ^ AN fq/p ] ( x ^ k ^) s[n (^~^s)-^D h/q (z,p_ L ), (84) 
^Collins = Y, e lJ d2k ^ A T^,k ± )^^A N D h/q ,(z,p ± )sm{<f) S + if + <f) h q ). (85) 



Let us explain in more detail the above equations. <ps and if identify the directions of the proton 
spin, S, and of the quark intrinsic transverse momentum, k±, see Fig. [14j <p q is the azimuthal angle 
of the final hadron h, as defined in the helicity frame of the fragmenting quark. da/dQ 2 is the usual 
unpolarized elementary cross section for the £q — > £q process, see also Eq. (177j) . while (y = Q 2 /xbs) 

d(Aa) da^^ da^^ 1 Ana 2 . 

-lQ 2 - = ^Q 2 dQi- = -Q^ {1 - y) - (86) 

The Sivers effect comes directly from the possible azimuthal dependence of the number density for 
unpolarized quarks inside a transversely polarized proton (generically denoted by p\ with p^ denoting 
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the opposite polarization state), that can be written as (see also section [37Tj) : 

f q/p r(x, k ± ) = f q/p {x, k ± ) + \ A N f g/p r(x, k ± ) S-(Pxk L ), (87) 

where the triple product gives in our kinematical configuration the sin((^ — 0$) factor entering Eq. (I84p . 
Upon integration over k^, A SiveTS gets the well known sm((p h — (fi s ) azimuthal dependence. 

Because of the kinematics, the Collins case is a bit more involved. Here is the transverse spin 
polarization of the fragmenting quark (as inherited by the initial polarized quark) that originates the 
possible azimuthal asymmetry in the dependent fragmentation process. Namely, for a transversely 
polarized quark with spin polarization s and three- momentum p q , 

D h /q, s (z, p±) = D h/g (z, p ± ) + - A N D h/q1 (z, p ± ) s-(p q x p ± ) . (88) 

The sm((j)s + f J r<f> q ) azimuthal dependence in Eq. ( |85l) arises from the combination of the phase factors 
appearing in the transversity distribution function, in the non-planar Iq — > Iq elementary scattering 
amplitudes, and in the Collins fragmentation function. Neglecting 0{k 2 _/Q 2 ) terms, one finds 

cos <j)= — cos((f) h - if) - z — , sin = — sin(0/> - ip) . (89) 
q Pi p± P± 

The integration upon k± results, among other factors, in the well known sin(^/ l + <ps) dependence 
of the Collins effect in SIDIS. 

The different azimuthal dependence of the Sivers and the Collins contributions in SIDIS offers 
indeed a clear separation technique and the possibility to gain valuable information on both of them. 
This is particularly true for the Sivers mechanism, since the other terms in Eq. (18~4"|) are well known. 
On the other hand, the Collins contribution involves also the still unknown transversity distribution 
function. We have seen in section 15.41 how a combined analysis of SIDIS results with data on azimuthal 
asymmetries in e + e~ — > tttt + X processes is of great help in this respect. 

To separate the contributions in the polarized cross section coming from the two naively T-odd 
effects, one defines the moments of the azimuthal asymmetry A sl ^ h±( ^ s \ given in Eq. ([2]), where in this 

case da(4>s) corresponds to (d 6 a ep ^ e ' h+x )/(dxBdydz h d 2 PTd<f)s), see Eq. fl8"51) . 

In order to perform analytically all k± integrations one can choose suitable Gaussian parameteriza- 
tions for the TMD distributions appearing in the above equations, see Ref. |268j : 

A T q(x,k ± ) = -M q T (x) [f q/p (x) + A L q(x)] 6 X X (90) 

A N f q / p i{z, k±) = 2M q Hx)f q/p (x)V^e^ e -^ r (91) 

A N D h/q ,(z,p ± ) = 2M q c (z)D h/q (z)V^e^ , (92) 
where Aiq(x) is the usual collinear helicity distribution, and 

(k 2 ) - M ' 2{kl) (n 2 ) - M2{PI) W 
( ^ >S " M' 2 + W) ' {PL)C "MH (p\) ' (93) 

with |A/T' 5 ' C | < 1 to automatically fulfill proper bounds. With these choices one gets 



A M*H-*S) fa ^ Zj p r) „ fivers ^ A *M<t> h +<f> S ) fa y f Zj p r) „ Aocollins ^ (g4) 
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where, omitting a common overall factor, 2ira 2 /xy 2 s, [recall that x ~ xb and z ~ Zh at (9(/cj_/Q)], 



A(T: 



Sivers 



with 



i + (i - v? £ e ? 2A/ ?( a; ) /«/*(*) ^W*) 



(95) 



Aa Collins = ^ ^JfP e "^? >C (1 ~y) E^A^(x) + A L g(x)j Mq{z) D h / q (z) (96) 



cr = 2?T 



(pi) (PI) 



7r(P|) 



e -P|/(P|) 



(97) 



(^>ff=(*i>5 + ^(*i>, = (pl>c + * 2 (*i> , (98) 

and (P|) given in Eq. (l8"Ti). Eq. (l9"5|l shows that A^ h ~^ s) = when z = or P T = 0. Analogously, 
Eq. (EED shows that Af^ h+4>s) = for P T = 0. 




Figure 16: Left panel: HERMES data on A^ 4 *'*^ [78] for eV -> e + n ± + X. Right panel: 
COMPASS data [82] on A^ h "* s) for /iP T fiht + X. Curves are the results of the fit 
in Ref. [263]. Predictions for it asymmetries (left panel, upper plot) are also shown. The 
shaded areas span a region corresponding to one-sigma deviaton at 90% CL. Unpolarized 
PDF and FF sets: MRST01-K. 



Concerning the unknown functions, J\f q , in Refs. [2631 1268"] forms like 

Af q (x) = N q x a "(l - x) b "(a q + b q ) a « +b «/a a q «b b q « (99) 

have been adopted. In particular, for the Sivers function only u and d quarks were considered, with a 
total of 7 free parameters (including M'); for the u and d quark transversity distributions, a u = and 
b u = bd have been adopted, keeping the flavour distinction only in N q , besides the one coming from the 
flavour structure in Eq. (1901) (4 free parameters). Finally, for the Collins function, in the pion case both 
favoured and unfavoured fragmentation functions have been considered, with the same a, b but different 
iVf av and A^unf for a total of 5 parameters (including M). 

The curves obtained by fitting HERMES (COMPASS) data on Ag$* h ~* s) for pion (charged hadron) 
production are shown in Fig. [161 Preliminary data from HERMES [81] and COMPASS [HH [85] on kaon 
asymmetries are also available. The corresponding theoretical estimates, adopting the Sivers functions 
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Figure 17: Left panel: HERMES data [78] on the azimuthal asymmetry A s ™^ s+< ^ h " > for ^ 
production. Right panel: COMPASS data [53] on A^j^ s+ ^ h \ for charged hadron production 
in /j,D^ — > fih + X. The extra tt phase in the figure label keeps into account the different 
choice w.r.t. Trento conventions [73] • The curves are obtained from Eq. fl96l) with the pa- 
rameterizations of Ref . [268] . The shaded areas correspond to the theoretical uncertainty on 
the parameters. (Un)polarized PDF and FF sets: (GRV98)GRSV2000-K. 



so extracted, are in fair agreement with HERMES data, even if the K + results are underestimated in 
the low x region. Further work is needed, in particular concerning the study of the sea-quark Sivers 
functions and flavour decomposition in the fragmentation sector. 

For the Collins asymmetry a global fit on HERMES, COMPASS and Belle data in SIDIS and e + e~ 
annihilation has been performed (9 free parameters) allowing the first extraction of the transversity 
distribution. The results for the SIDIS case and a comparison with the HERMES and COMPASS data 
are presented in Fig. [T71 
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Figure 18: HERMES data [ST] on A^ 5 ^ for K ± , compared with the results of Ref. [268]. 
(Un)polarized PDF and FF sets: (GRV98)GRSV2000-K. 



Finally, we compare the azimuthal asymmetry A^^s+M f or fo e production of K mesons with ex- 
isting HERMES results [77J. Preliminary results are also available from COMPASS [HJ[S5]. These data 
have not been included in the best fit of Ref. [263], as they might involve the transversity distribution 
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of strange quarks in the nucleon, which have been neglected for SIDIS data on tx production. We show 
in Fig. [TS] the results obtained using the extracted u and d transversity distributions of Fig. [2UJ (a, b 
panels). For the Collins function the same and M as extracted for the pion case were tentatively 
used, with the appropriate unpolarized FF's, D K ±/ q . 

We notice that the above computations are in fair agreement with data concerning the K + produc- 
tion, which is presumably dominated by u quarks; instead, there seem to be discrepancies for the K~ 
asymmetry, for which the role of s quarks might be relevant. New data on the azimuthal asymmetry 
for K production, possible from COMPASS and JLab experiments, might be very helpful in sorting out 
the eventual importance of the sea quark transversity distributions in the nucleon. 

5.5.3 Phenomenological extractions of Sivers and Collins functions 

The phenomenological analysis of HERMES and COMPASS data on moments of the azimuthal asymme- 
tries in the SIDIS process, with the aim of extracting information on the Sivers and Collins functions, has 
been performed independently by three different theoretical groups |238[ 12801 12761 1263[ 12501 1267[ [268J . 
Although the three groups have used different theoretical approaches, at the level of accuracy consid- 
ered (calculations are performed within LO pQCD, using leading twist TMD parton distribution and 
fragmentation functions and keeping only the leading term in a (k±/Q) power expansion) the results are 
quite similar and easily comparable. Several other simplifying approximations were common between 
these groups, namely the neglect of the so-called soft-factors [1221 1146] and of Sudakov suppression 
factors [149] (see also section [3]). However, slightly different approaches were followed concerning the 
dependence of the distributions on the transverse parton momenta. As an example, different functional 
forms were used. Also, possible constraints coming from general bounds (e.g. positivity bounds, the 
Soffer bound) were imposed and used in different ways. Regarding the Collins function, in Ref. [267j a 
quark soliton model was used for the transversity distribution, while in Ref. [250] the Soffer bound was 
tentatively adopted. As we have already discussed, in Ref. [268] Belle e + e~ data on azimuthal asymme- 
tries in dihadron production were used to independently extract the Collins function. From there and 
HERMES and COMPASS data, the first extraction of the transversity distribution for u and d quarks 
was performed, see Fig. [2D] ((a), (b) panels). Also, different sets for the usual fc_|_-integrated PDF's and 
FF's, in different combinations, were adopted. A useful summary and more detailed comparisons may 
be found, in the case of the Sivers function, in Ref. [262J. 




0.2 0.4 0.6 0.8 x 0.2 0.4 0.6 0.8 x 



Figure 19: The first and 1 /2-transverse moments of the quark Sivers functions, see Eq. ( TTTj) . 
as extracted in Refs. [280[ 12631 1250] . The curves indicate the 1-a regions of the various 
parameterizations. 

In Fig. [19] the first two moments (in k±) of the Sivers function (multiplied by x), as extracted by 
the three theoretical groups, are shown for comparison (see also Eqs. (fT6|) . (TTTj) ). Notice that for each 
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Figure 20: The transversity distribution for u and d quarks and the favoured and unfavoured 
Collins FF's from Ref. [268] . (a): x A T u(x) (upper plot) and x A T d(x) (lower plot), as 
functions of x at Q 2 = 2.4 (GeV/c) 2 . The Soffer bound [182] is also shown (solid blue lines), 
(b): unintegrated transversity distributions, xAtu(x, k±) (upper plot) and x Ard(x, k±) 
(lower plot), as defined in Eq. fl90l . as functions of k± at x — 0.1. PDF sets: GRV98- 
GRSV2000. (c): the p± integrated Collins functions, Eq. (1761) . normalized to twice the 
corresponding unpolarized FF's, vs. z; results of Refs. |267] (dashed line) and |250] (dotted 
line) are also shown, (d): the Collins function from the fits of Ref. [268] vs. p± at z = 0.36 
and Q 2 = 2.4 (GeV/c) 2 (lower lines). The positivity bound is also shown (upper blue lines). 
FF set: K. The shaded areas correspond to the theoretical uncertainty on the parameters. 



parameterization the two curves delimiting the one-sigma deviation region are shown. The range of 
x values spanned by HERMES experiment is also depicted. Clearly, the parameterizations extracted 
suffer of larger uncertainties in the large x sector, which is not well covered by data. Roughly speaking, 
all parameterizations are in favour of a negative (positive) Sivers function (f^p) for u(d) quarks, with 
shapes and magnitudes comparable. This behaviour reasonably complies with the prediction [281] 1282] 
f£r = —fyr, also derived in the large Nc limit of QCD in Ref. [283] . See also the discussion on the 
theoretical models in section 13.1.31 

Fig. [201 ((c). (d) panels) shows the first moment of the favoured and unfavoured Collins fragmentation 
functions for pions, normalized to twice the corresponding unpolarized FF, as a function of z ((c) panel) 
and p± ((d) panel) at fixed Q 2 . The three theoretical extractions are compared: Ref. [263] (solid lines 
and the shaded, one-sigma deviation, area); Ref. [267] (dashed line); Ref. [250] (dotted line). For 
comparison, we also show the positivity bound (upper solid line). The p± dependence in panel (d) is 
not a prediction but has been imposed to the parameterization. Notice also that there is no curve from 
Ref. [250] here: this is because of some additional simplifying assumptions performed there. 



6 Conclusions and outlook 

In this paper we have presented a thorough review of our present theoretical understanding of azimuthal 
and transverse single spin asymmetries in high-energy hadronic processes in the framework of pertur- 
bative QCD. The basic features of the main theoretical approaches to SSA's and an overview of the 
available experimental information have been given. An in-depth treatment of the generalized parton 
model, including spin and transverse momentum effects within the helicity formalism, and a detailed 
discussion of the phenomenology concerning azimuthal and SSA's complete this analysis. Hopefully, 
this work will provide the reader with a useful overview of the present status of this field of research and 
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with a fresh feeling of its ever-growing development. In fact, in the past few years several important 
results have been reached, both from the experimental and theoretical points of view. These results 
have confirmed earlier indications that spin and intrinsic transverse momentum effects can be relevant 
even at high energies and large pt and led to a profusion of interesting physical consequences. In the 
last years, a wide dedicated experimental activity has definitely confirmed previous pioneering results, 
refining their quality and statistical significance and enlarging the kinematical regime explored and the 
number of experimental observables investigated. From the theoretical side, significant progress has 
been reached in the understanding of new spin and TMD functions, like the Sivers distribution and the 
Collins function. The crucial role of initial and final state interactions, of a proper account of gauge 
links and QCD colour gauge invariance and of their consequences has been clarified. Related to this, 
the universality of the new TMD functions both in the distribution and in the fragmentation sector has 
been investigated. QCD factorization theorems and collinear twist-three approaches have been extended 
to cover several processes and new kinematical regimes. The theoretical formulation of the generalized 
parton model has been refined and its phenomenological applications extended to include most of the 
new experimental data. The interplay between twist-three approaches and QCD-improved parton mod- 
els with inclusion of intrinsic parton motion has been extensively studied. The large amount of activity 
devoted to SSA's and its significant but partial successes have also naturally raised new theoretical 
problems and experimental challenges. They will hopefully set up the stage for new developments in 
the forthcoming years. 

From the theoretical side, some of the most relevant open points are: 1) The formal properties of spin 
and k± dependent functions, their universality and evolution with scale; 2) The systematic comparison 
among different theoretical approaches and a deeper understanding of their interplay and complemen- 
tarity; 3) The role of Sudakov and soft factors, and of parton off-shellness and doubly unintegrated 
parton distributions; 4) A better understanding of the physical mechanisms involved, according to their 
twist classification and their behaviour with the relevant scales; 5) A comprehensive phenomenologi- 
cal analysis of all experimental data available as a test of universality and of the effectiveness of the 
proposed parameterizations for all the new unknown soft functions involved; 6) An improvement of 
nonperturbative models and approaches, including lattice calculations, for the same functions; 7) An 
improved understanding of the connections among TMD functions, generalized parton distributions and 
the orbital angular momentum of partons and its role for hadronic structure. 

From the experimental side, interesting developments are awaited in the near future: 1) The 
BRAHMS, PHENIX and STAR collaborations at RHIC-BNL will provide us with new results at 62, 
200, and 500 GeV cm. energies, with larger statistics, an improved kinematical coverage and includ- 
ing new interesting physical observables; 2) Hopefully, the COMPASS Collaboration will soon measure 
azimuthal single spin asymmetries for several charged hadrons with a transversely polarized proton 
target; on a more extended time scale they will also study the polarized Drell-Yan process with a pion 
beam colliding on transversely polarized proton and deuteron targets; 3) The proposed JLAB 12 GeV 
upgrade, with polarized proton, neutron and deuteron targets and a large-acceptance detector will pro- 
vide a detailed kinematical coverage in both longitudinal and tranverse momentum, with a large xb 
experimental reach; 4) The foreseen fixed-target experiments at y/s ~ 10 GeV at the JPARC facility 
will help in testing spin asymmetries and TMD distributions in the valence region (xb > 0.1); 5) The 
proposed 70 GeV upgrade of the transversely polarized proton beam at the IHEP-Protvino accelerator 
will provide further information on the moderate energy regime; 6) On a longer time scale, the PAX 
experiment at GSI will improve our knowledge of the transversity distribution and the other spin and 
TMD functions, provided the capability of efficiently polarizing antiprotons will be confirmed. 

Hopefully, the new forthcoming experimental results and the prosecution of active theoretical in- 
vestigation will lead in the near future to significant improvements in our understanding of SSA's, the 
spin structure of hadrons and related issues. 



62 



References 



[i] 

[2] 
[3] 
[4] 

[5] 
[6] 
[7] 
[8] 
[9] 
[10] 
[11] 

[12] 
[13] 
[14] 
[15] 
[16] 
[17] 
[18] 
[19] 
[20] 
[21] 

[22] 

[23] 
[24] 
[25] 
[26] 
[27] 
[28] 
[29] 
[30] 
[31] 
[32] 
[33] 
[34] 
[351 



C. Bourrely, J. Soffer, and E. Leader, Phys. Rep. 59 (1980) 95 
L. G. Pondrom, Phys. Rep. 122 (1985) 57 

M. Anselmino, A. Efremov, and E. Leader, Phys. Rep. 261 (1995) 1, |hep-ph/9501369 



S. M. Troshin and N. E. Tyurin, Spin phenomena in particle interactions, World Scientific, 
Singapore (1994). 

U. Stiegler, Phys. Rep. 277 (1996) 1 

K. J. Heller, Spin96 Proceedings, World Scientific, Singapore (1997), p. 23. 
S. B. Nurushev, Int. J. Mod. Phys. A 12 (1997) 3433 
J. Felix, Mod. Phys. Lett. A 14 (1999) 827 



Z.-t. Liang and C. Boros, Int. J. Mod. Phys. A 15 (2000) 927, hep-ph/0001330 



B. Lampe and E. Reya, Phys. Rep. 332 (2000) 1, |hep-ph/98 10270 



G. Bunce, N. Saito, J. Soffer, and W. Vogelsang, Ann. Rev. Nucl. Part. Set. 50 (2000) 525, 
|hep-ph/000"72i~8 

V. Barone, A. Drago, and P. G. Ratcliffe, Phys. Rep. 359 (2002) 1, |hep-ph/0104283| 
E. Leader, Spin in particle physics, Cambridge University Press, Cambdridge (2001). 



B. W. Filippone and X.-D. Ji, Adv. Nucl. Phys. 26 (2001) 1, |hep-ph/Oloi224 



V. Barone and P. G. Ratcliffe, Transverse spin physics, World Scientific, River Edge (2003). 



S. D. Bass, Rev. Mod. Phys. 77 (2005) 1257, |hep-ph/0411005 
W. Vogelsang, J. Phys. G 34 (2007) S149 
D. W. Sivers, Phys. Rev. D 41 (1990) 83 
D. W. Sivers, Phys. Rev. D 43 (1991) 261 
J. C. Collins, Nucl. Phys. B 396 (1993) 161, |hep-ph/9208213| 

J. C. Collins, S. F. Heppelmann, and G. A. Ladinsky, Nucl. Phys. B 420 (1994) 565, 
|hep-ph/9305309| 

D. Mueller, D. Robaschik, B. Geyer, F. M. Dittes, and J. Horejsi, Fortschr. Phys. 42 (1994) 



101, |hep-ph/9812448| 



X.-D. Ji, Phys. Rev. Lett. 78 (1997) 610, |hep~-ph/9603249 



X.-D. Ji, Phys. Rev. D 55 (1997) 7114, hep-ph/9609381 



A. V. Radyushkin, Phys. Rev. D 56 (1997) 5524, hep-ph/9704207 



X.-D. Ji, J. Phys. G 24 (1998) 1181, hep-ph/9807358 



M. Diehl, Phys. Rep. 388 (2003) 41, |hep~-ph/0307382 



A. V. Belitsky and A. V. Radyushkin, Phys. Rep. 418 (2005) 1, |hep-ph/0504030| 



M. Diehl, A. Manashov, and A. Schafer, Eur. Phys. J. A 31 (2007) 335, |hep-ph/06iTT01 



M. Burkardt, Int. J. Mod. Phys. A 18 (2003) 173, hep-ph/0207047 



M. Burkardt, Phys. Rev. D 66 (2002) 1 14005, |hip~-ph/0209179 



M. Burkardt and D. S. Hwang, Phys. Rev. D 69 (2004) 074032, [hip-ph/0309072"1 
M. Burkardt, Phys. Rev. D 72 (2005) 094020, [hep-ph/0505189 



M. Diehl and P. Hagler, Eur. Phys. J. C 44 (2005) 87, |hep-ph/0504175 



S. Meissner, A. Metz, and K. Goeke, Phys. Rev. D 76 (2007) 034002, [hep-ph/0703176 



63 



[36] HERMES Collaboration, P. B. van der Nat, Transversity 2005, World Scientific, Singapore 
(2006), hep-ph/0512019 

COMPASS Collaboration, C. Schill, DIS 2007, arXiv: 0706. 1459 [hep-ex] 
A. V. Efremov, L. Mankiewicz, and N. A. Tornqvist, Phys. Lett. B 284 (1992) 394 
J. C. Collins and G. A. Ladinsky, |hep-ph/9411444| 
X. Artru and J. C. Collins, Z. Phys. C 69 (1996) 277, |hep-ph/9504220| 
R. L. JafTe, X.-m. Jin, and J. Tang, Phys. Rev. Lett. 80 (1998) 1166, |hep-ph/9709322| 
A. Bianconi, S. Boffi, R. Jakob, and M. Radici, Phys. Rev. D 62 (2000) 034008, 



hep-ph/9907475 



A. Bacchetta and M. Radici, Phys. Rev. D 69 (2004) 074026, hep-ph/0311173 



F. A. Ceccopieri, M. Radici, and A. Bacchetta, Phys. Lett. B 650 (2007) 81, |hep-ph/070 3265 

G. Bunce et al, Phys. Rev. Lett. 36 (1976) 1113 

G. L. Kane, J. Pumplin, and W. Repko, Phys. Rev. Lett. 41 (1978) 1689 
E581 Collaboration, D. L. Adams et al, Phys. Lett. B 261 (1991) 201 
E704 Collaboration, D. L. Adams et al, Phys. Lett. B 264 (1991) 462 
E581 Collaboration, D. L. Adams et al, Z. Phys. C 56 (1992) 181 
E704 Collaboration, A. Bravar et al, Phys. Rev. Lett. 77 (1996) 2626 
E704 Collaboration, D. L. Adams et al, Phys. Rev. D 53 (1996) 4747 
E704 Collaboration, D. L. Adams et al, Phys. Lett. B 345 (1995) 569 
E704 Collaboration, A. Bravar et al, Phys. Rev. Lett. 75 (1995) 3073 
E704 Collaboration, D. L. Adams et al, Nucl. Phys. B 510 (1998) 3 
V. D. Apokin et al, Phys. Lett. B 243 (1990) 461 



V. V. Abramov et al, Nucl. Phys. B 492 (1997) 3, |hep-ex/0110011 
K. Krueger et al, Phys. Lett. B 459 (1999) 412 
C. E. Allgower et al, Phys. Rev. D 65 (2002) 092008 

STAR Collaboration, J. Adams et al, Phys. Rev. Lett. 92 (2004) 171801, [hep-ex/0310058 



STAR Collaboration, J. Adams et al, Phys. Rev. Lett. 97 (2006) 152302, |nucl-ex/06020TTl 



D. de Florian, R. Sassot, and M. Stratmann, Phys. Rev. D 75 (2007) 114010, [hep-ph/0703242 
STAR Collaboration, B. I. Abelev et al, Phys. Rev. Lett. 99 (2007) 142003, 



arXiv: 0705. 4629 [hep-ex] 



STAR Collaboration, L. Nogach, AIP Conf. Proc. 915 (2007) 543, |hep-ex/0612030 
l, F. Simon, AIP Conf. Proc. 870 (2006) 428, |hep-ex/0608050 



PHENIX Collaboration, S. S. Adler et al, Phys. Rev. Lett. 91 (2003) 241803, hep-ex/0304038 



PHENIX Collaboration, S. S. Adler et al, Phys. Rev. Lett. 95 (2005) 202001, [hep-ex/0507073 
PHENIX Collaboration, A. Adare et al, |arXiv: 0704 .3599 [hep-ex]~| 



PHENIX Collaboration, S. S. Adler et al, Phys. Rev. D 71 (2005) 071102, hep-ex/0502006 



PHENIX Collaboration, S. S. Adler et al, Phys. Rev. Lett. 98 (2007) 012002, [hepjjex/0609031 
PHENIX Collaboration, C. Aidala, DNP Spin Workshop 2006, unpublished 



BRAHMS Collaboration, I. Arsene et al, Phys. Rev. Lett. 98 (2007) 252001, |hep-ex/070104T 



BRAHMS Collaboration, J. H. Lee and F. Videbaek, AIP Conf. Proc. 915 (2007) 533 

A. Bacchetta, U. D'Alesio, M. Diehl, and C. A. Miller, Phys. Rev. D 70 (2004) 117504, 
|hep-ph/0410050| 



61 



[74] 
[75] 
[76] 
[77] 

[78] 
[79] 
[80] 
[81] 
[82] 

[83] 
[84] 
[85] 
[86] 
[87] 
[88] 
[89] 

[90] 

[91] 
[92] 
[93] 
[94] 
[95] 
[96] 
[97] 
[98] 
[99] 
[100] 
[101] 
[102] 
[103] 
[104] 
[105] 
[106] 
[107] 

[108] 
[109] 
[110] 



HERMES Collaboration, A. Airapetian et al, Phys. Rev. Lett. 84 (2000) 4047, |hep-ex/9910062 



HERMES Collaboration, A. Airapetian et al, Phys. Rev. D 64 (2001) 097101, [hep-ex/010 4005 
HERMES Collaboration, A. Airapetian et al, Phys. Lett. B 562 (2003) 182, |hep-ex/0212039 
HERMES Collaboration, A. Airapetian et al, Phys. Rev. Lett. 94 (2005) 012002, 
|hep-ex/0408013| 



HERMES Collaboration, M. Diefenthaler, AIP Conf. Proc. 792 (2005) 933, hep-ex/0507013 



HERMES Collaboration, A. Airapetian et al, Phys. Lett. B 622 (2005) 14, hep-ex/0505042 



HERMES Collaboration, A. Airapetian et al, Phys. Lett. B 648 (2007) 164, hep-ex/0612059 



HERMES Collaboration, M. Diefenthaler, AIP Conf. Proc. 915 (2007) 509, |hep-ex/06 12010 



COMPASS Collaboration, V. Y. Alexakhin et al, Phys. Rev. Lett. 94 (2005) 202002, 



hep-ex/0503002 



COMPASS Collaboration, E. S. Ageev et al, Nucl. Phys. B 765 (2007) 31, |hep-ex/0610068 



COMPASS Collaboration, A. Martin, Czech. J. Phys. 56 (2006) F33, |hep-ex/070 2002 



COMPASS Collaboration, F. Bradamante, AIP Conf. Proc. 915 (2007) 513, |hep-ex/0702007 
COMPASS Collaboration, A. Kotzinian, DIS 2007, arXiv: 0705. 2402 [hep-ex] 



A. Bacchetta et al, JEEP 02 (2007) 093, hep-ph/06 11265 



CLAS Collaboration, H. Avakian et al, Phys. Rev. D 69 (2004) 112004, [hep-ex/030100"5| 

CLAS Collaboration, H. Avakian, P. Bosted, V. Burkert, and L. Elouadrhiri, Transversity 2005, 
World Scientific, Singapore (2006), p. 88. 

CLAS Collaboration, H. Avakian, P. E. Bosted, V. Burkert, and L. Elouadrhiri, AIP Conf. 
Proc. 792 (2005) 945, |nucl-ex/0509032| 



M. Ahmed and T. Gehrmann, Phys. Lett. B 465 (1999) 297, |hep-ph/9906503| 
European Muon Collaboration, M. Arneodo et al, Z. Phys. C 34 (1987) 277 
European Muon Collaboration, J. J. Aubert et al, Phys. Lett. B 130 (1983) 118 
E665 Collaboration, M. R. Adams et al, Phys. Rev. D 48 (1993) 5057 



ZEUS Collaboration, J. Breitweg et al, Phys. Lett. B 481 (2000) 199, hep-ex/0003017 



ZEUS Collaboration, S. Chekanov et al, Phys. Lett. B 551 (2003) 226, |hep-ex70210064 



ZEUS Collaboration, S. Chekanov et al, Eur. Phys. J. C 51 (2007) 289, |hep-ex/0608053 

C. S. Lam and W.-K. Tung, Phys. Rev. D 18 (1978) 2447 

C. S. Lam and W.-K. Tung, Phys. Lett. B 80 (1979) 228 

C. S. Lam and W.-K. Tung, Phys. Rev. D 21 (1980) 2712 

NA3 Collaboration, J. Badier et al, Z. Phys. C 11 (1981) 195 

NA10 Collaboration, S. Falciano et al, Z. Phys. C 31 (1986) 513 

NA10 Collaboration, M. Guanziroli et al, Z. Phys. C 37 (1988) 545 

E. L. Berger and S. J. Brodsky, Phys. Rev. Lett. 42 (1979) 940 

E. L. Berger, Z. Phys. C 4 (1980) 289 

J. S. Conway et al, Phys. Rev. D 39 (1989) 92 

E866/NuSea Collaboration, L. Y. Zhu et al, Phys. Rev. Lett. 99 (2007) 082301, 
|hep-ex/0609005| 



D. Boer and P. J. Mulders, Phys. Rev. D 57 (1998) 5780, |hep-ph/971 1485 



Belle Collaboration, K. Abe et al, Phys. Rev. Lett. 96 (2006) 232002, |hep :: ex70507063 



A. Ogawa, M. Grosse-Perdekamp, R. Seidl, and K. Hasuko, DIS 2006, hep-ex/0607014 



65 



Ill] Belle Collaboration, A. Ogawa, M. Grosse Perdekamp, R.-C. Seidl, and K. Hasuko, AIP Conf. 
Proc. 915 (2007) 575 

N. I. Kochelev, JETP Lett. 72 (2000) 481, |hep-ph/9905497 



112 
113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 
131 
132 
133 
134 
135 
136 
137 

138 
139 
140 
141 
142 
143 

144 
145 
146 
147 
148 



D. Ostrovsky and E. Shuryak, Phys. Rev. D 71 (2005) 014037, [hep-ph/0409253 



P. Hoyer and M. Jarvinen, JEEP 10 (2005) 080, hep-ph/0509058 



P. Hoyer and M. Jarvinen, JEEP 02 (2007) 039, |hep-ph/06 11293 



C. Bourrely and J. Soffer, Eur. Phys. J. C 36 (2004) 371, |hep-ph/0311110| 

D. de Florian and W. Vogelsang, Phys. Rev. D 71 (2005) 114004, [hep-ph/050 1258] 
R. D. Field and R. P. Feynman, Phys. Rev. D 15 (1977) 2590 

R. P. Feynman, R. D. Field, and G. C. Fox, Phys. Rev. D 18 (1978) 3320 
J. C. Collins and D. E. Soper, Phys. Rev. D 16 (1977) 2219 
J. C. Collins and D. E. Soper, Nucl. Phys. B 193 (1981) 381 

X.-d. Ji, J.-p. Ma, and F. Yuan, Phys. Rev. D 71 (2005) 034005, [hep-ph/0404183"1 



X.-d. Ji, J.-P. Ma, and F. Yuan, Phys. Lett. B 597 (2004) 299, [hep 1 ph/0405085 



M. Anselmino, M. Boglione, and F. Murgia, Phys. Lett. B 362 (1995) 164, [hep 1 ph/9503290 



M. Anselmino and F. Murgia, Phys. Lett. B 442 (1998) 470, hep-ph/9808426 



P. Mulders and R. Tangerman, Nucl. Phys. B 461 (1996) 197, |hep-ph/9510301 
J. P. Ralston and D. E. Soper, Nucl. Phys. B 152 (1979) 109 



X. Artru, J. Czyzewski, and H. Yabuki, Z. Phys. C 73 (1997) 527, |hep-ph/ 9508239 
B. Andersson, G. Gustafson, G. Ingelman, and T. Sjostrand, Phys. Rep. 97 (1983) 31 
J. Czyzewski, Acta Phys. Polon. 27 (1996) 1759, |hep-ph/9606390| 



M. Anselmino, M. Boglione, and F. Murgia, Phys. Rev. D 60 (1999) 054027, hep-ph/9901442 



A. Metz and M. Schlegel, Eur. Phys. J. A 22 (2004) 489, |hep-ph70403182 



A. Bacchetta, P. J. Mulders, and F. Pijlman, Phys. Lett. B 595 (2004) 309, hep-ph/0405154 



K. Goeke, A. Metz, and M. Schlegel, Phys. Lett. B 618 (2005) 90, hep-ph/0504130 



M. Anselmino et al, Phys. Rev. D 73 (2006) 014020, hep-ph/0509035 



D. Boer, Phys. Rev. D 60 (1999) 014012, hep-ph/9902255 



M. Anselmino, D. Boer, U. D'Alesio, and F. Murgia, Phys. Rev. D 63 (2001) 054029, 
|hep-ph/0008l86l 



S. J. Brodsky, D. S. Hwang, and I. Schmidt, Phys. Lett. B 530 (2002) 99, [hep 1 ph/0201296 



J. C. Collins, Phys. Lett. B 536 (2002) 43, hep-ph/0204004 



X.-d. Ji and F. Yuan, Phys. Lett. B 543 (2002) 66, hep-ph/0206057 



A. V. Belitsky, X. Ji, and F. Yuan, Nucl. Phys. B 656 (2003) 165, hep-ph/0208038 



J. C. Collins, Acta Phys. Polon. B 34 (2003) 3103, [hep 1 ph/0304122 



J. C. Collins, D. E. Soper, and G. Sterman, in Perturbative Quantum Chromodynamics, World 
Scientific, Singapore (1989), 1-92. 



S. J. Brodsky, D. S. Hwang, and I. Schmidt, Nucl. Phys. B 642 (2002) 344, |hep-ph/0206259 
A. Metz, Phys. Lett. B 549 (2002) 139 



J. C. Collins and A. Metz, Phys. Rev. Lett. 93 (2004) 252001, [hep-ph/0408249 



A. Idilbi, X.-d. Ji, J.-P. Ma, and F. Yuan, Phys. Rev. D 70 (2004) 074021, [hep-ph/0406302"1 



A. Idilbi, X.-d. Ji, and F. Yuan, Phys. Lett. B 625 (2005) 253, [hipph/0507196 



66 



[149] 
[150] 
[151] 
[152] 
[153] 
[154] 

[155] 
[156] 
[157] 
[158] 

[159] 

[160] 
[161] 
[162] 
[163] 

[164] 

[165] 
[166] 
[167] 
[168] 

[169] 
[170] 
[171] 
[172] 

[173] 
[174] 
[175] 
[176] 
[177] 

[178] 

[179] 

[180] 



D. Boer, Nucl. Phys. B 603 (2001) 195, |hep-ph/0102071 



A. A. Henneman, D. Boer, and P. J. Mulders, Nucl. Phys. B 620 (2002) 331, hep-ph/ 01 04271 



F. A. Ceccopieri and L. Trentadue, Phys. Lett. B 636 (2006) 310, |hep-ph/05 12372 
F. A. Ceccopieri and L. Trentadue, |arXiv: 0706 .4242 [hep-ph] | 



C. J. Bomhof, P. J. Mulders, and F. Pijlman, Phys. Lett. B 596 (2004) 277, |hep-ph/0406099l 

A. Bacchetta, C. J. Bomhof, P. J. Mulders, and F. Pijlman, Phys. Rev. D 72 (2005) 034030, 
|hep-ph/0505268| 



C. J. Bomhof, P. J. Mulders, and F. Pijlman, Eur. Phys. J. C 47 (2006) 147, hep-ph/0601171 



C. J. Bomhof and P. J. Mulders, JEEP 02 (2007) 029, |hep-ph/ 0609206 



J. Collins and J.-W. Qiu, Phys. Rev. D 75 (2007) 114014, arXiv : 0705 . 2141 [hep-ph] 



J.-W. Qiu, W. Vogelsang, and F. Yuan, Phys. Lett. B 650 (2007) 373, 
|arXiv: 0704 .1153 [hep^phj 



J.-W. Qiu, W. Vogelsang, and F. Yuan, Phys. Rev. D 76 (2007) 074029, 



arXiv: 0706 .1196 [hep-ph] 



J. Collins, arXiv: 0708. 4410 [hep-ph] 



W. Vogelsang and F. Yuan, Phys. Rev. D 76 (2007) 094013, arXiv: 0708. 4398 [hep-ph] 



C. J. Bomhof and P. J. Mulders, |arXiv: 0709. 1390 [hep-ph] 



A. Bacchetta, C. Bomhof, U. D'Alesio, P. J. Mulders, and F. Murgia, Phys. Rev. Lett. 99 



(2007) 212002, hep-ph/0703153 



C. J. Bomhof, P. J. Mulders, W. Vogelsang, and F. Yuan, Phys. Rev. D 75 (2007) 074019, 
|hep-ph/0701277| 



M. Burkardt, Phys. Rev. D 69 (2004) 057501, hep-ph/0311013 



D. Boer, S. J. Brodsky, and D. S. Hwang, Phys. Rev. D 67 (2003) 054003, [hep-ph/0211110 
G. R. Goldstein and L. Gamberg, |hep-ph/0209085| 

L. P. Gamberg, G. R. Goldstein, and K. A. Oganessyan, Phys. Rev. D 67 (2003) 071504, 
|hep-ph/0301018| 



F. Yuan, Phys. Lett. B 575 (2003) 45, hep-ph/0308157 



A. Bacchetta, A. Schaefer, and J.-J. Yang, Phys. Lett. B 578 (2004) 109, |hep-ph/0309246| 



Z. Lu and B.-Q. Ma, Nucl. Phys. A 741 (2004) 200, hep-ph/0406171 



I. O. Cherednikov, U. D'Alesio, N. I. Kochelev, and F. Murgia, Phys. Lett. B 642 (2006) 39, 



hep-ph/0606238 



S. J. Brodsky and F. Yuan, Phys. Rev. D 74 (2006) 094018, [hip-ph/0610236 



D. Amrath, A. Bacchetta, and A. Metz, Phys. Rev. D 71 (2005) 114018, [56^^0504124 



L. Gamberg, G. Goldstein, and M. Schlegel, arXiv: 0708. 0324 [hep-ph] 



M. Burkardt and B. Hannafious, arXiv: 0705. 1573 [hep-ph] 



A. Bacchetta, R. Kundu, A. Metz, and P. J. Mulders, Phys. Lett. B 506 (2001) 155, 



|hep-ph/0 102278 



A. Bacchetta, R. Kundu, A. Metz, and P. J. Mulders, Phys. Rev. D 65 (2002) 094021, 
hep-ph/020109l] 



L. P. Gamberg, G. R. Goldstein, and K. A. Oganessyan, Phys. Rev. D 68 (2003) 051501, 
|hep-ph/0307139| 



A. Bacchetta, A. Metz, and J.-J. Yang, Phys. Lett. B 574 (2003) 225, |hep-ph/0307282| 



67 



[181] 

[182] 
[183] 
[184] 
[185] 

[186] 
[187] 
[188] 
[189] 
[190] 
[191] 
[192] 

[193] 
[194] 
[195] 
[196] 
[197] 
[198] 
[199] 
[200] 
[201] 
[202] 

[203] 
[204] 
[205] 
[206] 
[207] 
[208] 
[209] 
[210] 
[211] 
[212] 
[213] 
[214] 
[215] 
[216] 
[217] 
[218] 



A. Bacchetta, M. Boglione, A. Henneman, and P. J. Mulders, Phys. Rev. Lett. 85 (2000) 712, 
|hep-ph/9912490| 

J. Soffer, Phys. Rev. Lett. 74 (1995) 1292, |hep-ph/9409254 



A. Schafer and O. V. Teryaev, Phys. Rev. D 61 (2000) 077903, hep-ph/9908412 



M. Burkardt, Phys. Rev. D 69 (2004) 091501, hep-ph/0402014 



K. Goeke, S. Meissner, A. Metz, and M. Schlegel, Phys. Lett B 637 (2006) 241, 
|hep-ph/0601133 



A. V. Efremov and O. V. Teryaev, Sov. J. Nucl. Phys. 36 (1982) 140 
A. V. Efremov and O. V. Teryaev, Sov. J. Nucl. Phys. 39 (1984) 962 
A. V. Efremov and O. V. Teryaev, Phys. Lett. B 150 (1985) 383 
J.-w. Qiu and G. Sterman, Phys. Rev. Lett. 67 (1991) 2264 
J.-w. Qiu and G. Sterman, Nucl. Phys. B 378 (1992) 52 



J.-w. Qiu and G. Sterman, Phys. Rev. D 59 (1998) 014004, hep-ph/9806356 



C. Kouvaris, J.-W. Qiu, W. Vogelsang, and F. Yuan, Phys. Rev. D 74 (2006) 114013, 



hep-ph/0609238 



Y. Koike and K. Tanaka, Phys. Rev. D 76 (2007) 011502, hep-ph/0703169 



Y. Kanazawa and Y. Koike, Phys. Lett. B 478 (2000) 121, hep-ph/0001021 



Y. Kanazawa and Y. Koike, Phys. Lett. B 490 (2000) 99, |hep-ph/0007272 
Y. Koike, Nucl. Phys. A 721 (2003) 364, |hep-ph/02 11400 



Y. Kanazawa and Y. Koike, Phys. Rev. D 64 (2001) 034019, hep-ph/00 12225 



H. Eguchi, Y. Koike, and K. Tanaka, Nucl. Phys. B 752 (2006) 1, hep-ph/0604003 



H. Eguchi, Y. Koike, and K. Tanaka, Nucl. Phys. B 763 (2007) 198, hep-ph/0610314 



D. Boer, P. J. Mulders, and F. Pijlman, Nucl. Phys. B 667 (2003) 201, hep-ph/0303034 



J. P. Ma and Q. Wang, Eur. Phys. J. C 37 (2004) 293, hep-ph/0310245 



X. Ji, J.-W. Qiu, W. Vogelsang, and F. Yuan, Phys. Rev. Lett. 97 (2006) 082002, 



hep-ph/0602239 



X. Ji, J.-w. Qiu, W. Vogelsang, and F. Yuan, Phys. Rev. D 73 (2006) 094017, hep-ph/0604023 



X. Ji, J.-W. Qiu, W. Vogelsang, and F. Yuan, Phys. Lett. B 638 (2006) 178, hep-ph/0604128 



Y. Koike, W. Vogelsang, and F. Yuan, arXiv: 0711 .0636 [hep-ph] 



G. Watt, A. D. Martin, and M. G. Ryskin, Phys. Rev. D 70 (2004) 014012, hep-ph/0309096 



O. Linnyk, S. Leupold, and U. Mosel, Phys. Rev. D 75 (2007) 014016, [hep-ph/0607305 
J. C. Collins, T. C. Rogers, and A. M. Stasto, |arXiv: 0708 . 2833 [hep^ph] 



Z.-T. Liang and T.-C. Meng, Z. Phys. A 344 (1992) 171 

C. Boros, Z. T. Liang, and T. C. Meng, Phys. Rev. Lett. 70 (1993) 1751 

Z.-t. Liang and T.-c. Meng, Phys. Rev. D 49 (1994) 3759 

C. Boros, Z. T. Liang, and T. C. Meng, Phys. Rev. D 51 (1995) 4867 



C. Boros, Z.-t. Liang, and T.-c. Meng, Phys. Rev. D 54 (1996) 4680, hep-ph/9603232 



H. Dong, F.-z. Li, and Z.-t. Liang, Phys. Rev. D 69 (2004) 017501, hep-ph/0311121 



C. Boros and Z.-t. Liang, Phys. Rev. D 53 (1996) 2279, hep-ph/9510319 



Z.-t. Liang and C. Boros, Phys. Rev. Lett. 79 (1997) 3608, |hep-ph7970 8488 
U. D'Alesio, S. Melis, and F. Murgia, unpublished, 2007[ 



U. D'Alesio and F. Murgia, Phys. Rev. D 70 (2004) 074009, hep-ph/0408092 



6n 



[219] M. Anselmino, M. Boglione, U. D'Alesio, E. Leader, and F. Murgia, Phys. Rev. D 71 (2005) 
014002, |hep-ph/0408356] 



P. J. Mulders and J. Rodrigues, Phys. Rev. D 63 (2001) 094021, hep-ph/0009343 



A. D. Martin, R. G. Roberts, W. J. Stirling, and R. S. Thorne, Phys. Lett. B 531 (2002) 216, 



hep-ph/0201127 



M. Gliick, E. Reya, and A. Vogt, Eur. Phys. J. C 5 (1998) 461, |hep-ph/ 9806404 

M. Gliick, E. Reya, M. Stratmann, and W. Vogelsang, Phys. Rev. D 63 (2001) 094005, 



hep-ph/0011215 



E. Leader, A. V. Sidorov, and D. B. Stamenov, Eur. Phys. J. C 23 (2002) 479, |hep-ph/01 11267 



B. A. Kniehl, G. Kramer, and B. Potter, Nucl. Phys. B 582 (2000) 514, hep -ph/00 10289 
S. Kretzer, Phys. Rev. D 62 (2000) 054001, |hep-ph/0003177 



A. P. Contogouris, R. Gaskell, and S. Papadopoulos, Phys. Rev. D 17 (1978) 2314 
C.-Y. Wong and H. Wang, Phys. Rev. C 58 (1998) 376, hep-ph/9802378| 
X.-N. Wang, Phys. Rev. C 61 (2000) 064910, |rmcl-th/981202Tl 

Y. Zhang, G. I. Fai, G. Papp, G. G. Barnafoldi, and P. Levai, Phys. Rev. C 65 (2002) 034903, 



hep-ph/0109233 



G. Donaldson et al, Phys. Lett. B 73 (1978) 375 
F. Yuan, |arXiv: 0709 .3272 [hep-ph] | 

U. D'Alesio, Talk delivered at the ECT* workshop "Transverse momentum, spin, and position 
distributions of partons in hadrons" Trento, June 11-15, 2007. 



PAX Collaboration, V. Barone et ai, hep-ex/0505054 



STAR Collaboration, C. A. Gagliardi, DIS 2006, hep-ex/0607003 

M. Anselmino, U. D'Alesio, S. Melis, and F. Murgia, Phys. Rev. D 74 (2006) 094011, 



hep-ph/0608211 



S. J. Brodsky and S. Gardner, Phys. Lett. B 643 (2006) 22, |hep-ph/060 8219 



A. V. Efremov, K. Goeke, S. Menzel, A. Metz, and P. Schweitzer, Phys. Lett. B 612 (2005) 233, 



hep-ph/0412353 



M. Anselmino, M. Boglione, U. D'Alesio, E. Leader, and F. Murgia, Phys. Rev. D 70 (2004) 



074025, hep-ph/0407100 



U. D'Alesio and F. Murgia, MP Conf. Proc. 915 (2007) 559, |hep-ph/0612208| 
S. J. Brodsky, P. Hoyer, C. Peterson, and N. Sakai, Phys. Lett. B 93 (1980) 451 
J. Pumplin, Phys. Rev. D 73 (2006) 114015, hep-ph/05081841 

I. Schmidt, J. Soffer, and J.-J. Yang, Phys. Lett. B 612 (2005) 258, |hep-ph/0503127 
K. J. Heller et al, Phys. Rev. Lett. 51 (1983) 2025 
B. Lundberg et al, Phys. Rev. D 40 (1989) 3557 
E. J. Ramberg et al, Phys. Lett. B 338 (1994) 403 



D. Indumathi, H. S. Mani, and A. Rastogi, Phys. Rev. D 58 (1998) 094014, hep-ph/9802324 



STAR Collaboration, B. I. Abelev et al, Phys. Rev. C 75 (2007) 064901, lnucl-ex/06070331 



D. Boer and W. Vogelsang, Phys. Rev. D 69 (2004) 094025, hep-ph/03 12320 



W. Vogelsang and F. Yuan, Phys. Rev. D 72 (2005) 054028, hep-ph/0507266 



A. Brandenburg, S. J. Brodsky, V. V. Khoze, and D. Mueller, Phys. Rev. Lett. 73 (1994) 939, 
|hep-ph/9403361| 



69 



[252] K. J. Eskola, P. Hoyer, M. Vanttinen, and R. Vogt, Phys. Lett. B333 (1994) 526, 



hep-ph/9404322 



O. Nachtmann and A. Reiter, Z. Phys. C 24 (1984) 283 

A. Brandenburg, O. Nachtmann, and E. Mirkes, Z. Phys. C 60 (1993) 697 

D. Boer, A. Brandenburg, O. Nachtmann, and A. Utermann, Eur. Phys. J. C 40 (2005) 55, 



hep-ph/0411068 



L. P. Gamberg and G. R. Goldstein, Phys. Lett. B 650 (2007) 362, |hep-ph/0506127 



N. Hammon, O. Teryaev, and A. Schafer, Phys. Lett. B 390 (1997) 409, hep-ph/96 11359 



D. Boer, P. J. Mulders, and O. V. Teryaev, Phys. Rev. D 57 (1998) 3057, hep-ph/9710223 



D. Boer and P. J. Mulders, Nucl. Phys. B 569 (2000) 505, |hep-ph79 906223 
D. Boer and J.-w. Qiu, Phys. Rev. D 65 (2002) 034008, [hep-ph/0108179 



M. Anselmino, U. D'Alesio, and F. Murgia, Phys. Rev. D 67 (2003) 074010, hep-ph/0210371 



M. Anselmino et al, Transversity 2005, World Scientific, Singapore (2006), hep-ph/0511017 



M. Anselmino et al, Phys. Rev. D 72 (2005) 094007, hep-ph/0507181 



J. C. Collins et al, Phys. Rev. D 73 (2006) 094023, hep-ph/05 11272 



S. J. Brodsky, D. S. Hwang, and I. Schmidt, Phys. Lett. B 553 (2003) 223, |hep-ph/02lT212 



D. Boer, R. Jakob, and P. J. Mulders, Nucl. Phys. B 504 (1997) 345, |hep z ph/9702281 



A. V. Efremov, K. Goeke, and P. Schweitzer, Phys. Rev. D 73 (2006) 094025, [hip~-ph/0603054l 



M. Anselmino et al, Phys. Rev. D 75 (2007) 054032, hep-ph/0701006 



M. Anselmino, M. Boglione, U. D'Alesio, F. Murgia, and A. Prokudin, unpublished, 2007 

R. N. Cahn, Phys. Lett. B 78 (1978) 269 

R. N. Cahn, Phys. Rev. D 40 (1989) 3107 

A. Konig and P. Kroll, Z. Phys. C 16 (1982) 89 

J.-g. Chay, S. D. Ellis, and W. J. Stirling, Phys. Rev. D 45 (1992) 46 

H. Georgi and H. D. Politzer, Phys. Rev. Lett. 40 (1978) 3 



A. Daleo, D. de Florian, and R. Sassot, Phys. Rev. D 71 (2005) 034013, hep -ph/0411212 
M. Anselmino et al, Phys. Rev. D 71 (2005) 074006, |hep-ph/050119"6 



M. Anselmino, M. Boglione, A. Prokudin, and C. Turk, Eur. Phys. J. A 31 (2007) 373, 



hep-ph/0606286 



LHPC Collaboration, B. U. Musch et al, [arXiv: 0710 .4423 [hep-lat] 



V. Barone, Z. Lu, and B.-Q. Ma, Phys. Lett. B 632 (2006) 277, hep-ph/0512145 



J. C. Collins et al, Phys. Rev. D 73 (2006) 014021, hep-ph/0509076 



M. Anselmino, V. Barone, A. Drago, and F. Murgia, hep-ph/0209073 



A. Drago, Phys. Rev. D 71 (2005) 057501, hep-ph/0501282 



P. V. Pobylitsa, hep-ph/0301236 



70 



